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Abstract 

A new, global inversion is used to estimate the large scale oceanic circulation based on 
the World Ocean Circulation Experiment and Java Australia Dynamic Experiment hydro- 
graphic data. A linear inverse “box” model is used to combine consistently the transoceanic 
sections. The circulation is geostrophic with an Ekman layer at the surface and oceanic 
layers defined by neutral surfaces. Near-conservation of mass, salt and top-to-bottom silica 
is required and, in addition, heat and the phosphate-oxygen combination (170[PO 4 ]+[O 2 ]) 
are conserved in layers that are not in contact with the surface. A globally-consistent solu¬ 
tion is obtained for a depth-independent adjustment to the thermal wind field, freshwater 
flux divergences, the Ekman transport, and the advective and diffusive dianeutral fluxes 
between layers. A detailed error budget permits calculation of statistical uncertainties, 
taking into account both the non-resolved part of the solution and the systematic errors 
due to the temporal oceanic variability. The estimated water mass transports during the 
WOCE period (1985-1996) are generally similar to previous published estimates. How¬ 
ever, important differences are found. In particular, the inflow of bottom waters into the 
Pacific Ocean is smaller than in most previous estimates. Utilization of property anomaly 
conservation constraints allows the estimation of significant dianeutral diffusivities in deep 
layers, with a global average of 3 ± lcm 2 s _1 north of 30°S. Dianeutral transfers indicate 
that about 20 Sv of bottom water is formed in the Southern Ocean. Significant ocean- 
atmosphere heat fluxes are found, with a global heating of 2.3 ± 0.4 PW in the tropical 
band and a corresponding cooling at high latitudes. The signature of a large-scale average 
export production is found for nutrients in several temperate regions. Despite the large 
uncertainties, the production magnitudes are consistent with independent measurements 
from sediment traps and isotopic data. Net nutrient sources or sinks are found in sev¬ 
eral regions, suggesting either transport of dissolved organic matter or a seasonal alias. 
Oxygen indicates large exchanges with the atmosphere, with intake at high latitudes and 
outgassing/remineralization at low latitudes. 


Thesis Supervisor: Carl Wunsch 

Title: Cecil and Ida Green Professor of Physical Oceanography 
Massachusetts Institute of Technology 
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Overview 


By its ability to transport important quantities of energy and dissolved carbon dioxide, 
the ocean plays a major role in climate regulation. One of the main objectives of the 
World Ocean Circulation Experiment (WOCE) was the determination of such transports. 

Oceanic surface waters are heated in low latitudes and transported by ocean currents 
towards higher latitudes. As a result, the ocean carries about the same amount of ther¬ 
modynamic energy 1 as does the atmosphere from the equator to the poles. Estimates of 
oceanic energy transports can be obtained by integration of the residuals from atmospheric 
models, but the uncertainties associated with this method are relatively large (e.g., Keith, 
1995). Alternatively, the energy transport can be derived from oceanic currents across 
transoceanic hydrographic sections. This later method is believed to be more accurate, 
and allows uncertainty calculations (e.g., Wijffels et al., 1999; Robbins and Toole, 1997; 
Macdonald and Wunsch, 1996; Bryden et al., 1991; Rintoul and Wunsch, 1991; Hall and 
Bryden, 1982, and many others). 

Similarly, carbon dioxide (CO 2 ) and oxygen (O 2 ) are absorbed and released at the 
ocean surface. Atmospheric C0 2 concentration is one of the key variables that affect 
the Earth s temperature and climate. The ocean is believed to be an important sink of 
anthropogenic C0 2 , in particular, the photosynthetic activity in the surface layer that 
consumes dissolved C0 2 fixes carbon and exports it towards the deep layers, through this 
“biological pump” effect (Sarmiento et al., 1990). The resulting depletion of C0 2 in the 
surface layer increases the amount of atmospheric C0 2 transferred to the ocean relative 
to a hypothetical ocean with no biological activity. Simulations suggest that this process 
may be lowering the atmospheric C0 2 concentration by a factor of 2 (e.g., Sarmiento and 
Toggweiler, 1984). 

The C0 2 export is thus intimately related to nutrient transport in the ocean, which can 
be calculated from hydrographic data. Nutrients are consumed by planktonic organisms 
in the surface layers of the ocean. They in turn are converted to particulate and organic 
matter, go through the food chain and end up in sinking fecal pellets and skeletons. 

1 Often mislabeled “heat”. 
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Some of the particles are recycled in the euphotic zone (regenerated production) while 
the remainder sink towards the deep oceanic layers (export production) where most are 
remineralized. These dissolved nutrients then flow slowly with the currents until they 
reach the surface to again support production. 

Existing estimates of new and export production are very uncertain. Because of large 
spatial and temporal variability in biological activity, direct measurements of export pro¬ 
duction and remineralization are difficult to interpret. Most estimates are local (e.g., 
Sarmiento et al., 1990; Spitzer and Jenkins, 1989) and rely upon uncertain assumptions, 
such as a one-dimensional balance over the water column. On the other hand, it is possible, 
in principle, to estimate new and export production from an accurate dissolved nutrient 
budget in the ocean. An estimate of oceanic nutrient transports should reveal, in the 
presence of export production, a dissolved nutrient depletion in the euphotic zone and a 
dissolved nutrient supply at deeper levels. 

In this thesis, I seek a global circulation scheme that is consistent with oceanic dynam¬ 
ics and simple biology and geochemistry. The benefits are two-fold: on the one hand, the 
dynamics provide information about the nutrient fluxes; on the other hand, the nutrient 
fluxes are able to provide constraints on the circulation. This is a first attempt to use the 
WOCE hydrographic data in such a global and consistent framework. 

Inverse methods provide convenient tools to quantify an ocean circulation that obeys 
simple conservation statements. Since the pioneering works of Riley (1951) and later Wun- 
sch(1978), steady geostrophic box inversions have been extensively used in most oceanic 
regions (see Wunsch, 1996 for a review). The method used here follows the earlier global 
calculations of Macdonald and Wunsch (1996) and Macdonald (1998) with a number of 
significant improvements, including the following: 

1. The present data consist almost entirely of recent, high-quality hydrographic sections 
from the WOCE survey. 

2. A more quantitative estimate of the uncertainties in property transports is obtained 
through an analysis of 1 /4° nominal lateral resolution General Circulation Model run 
at the US Naval Post Graduate School (see Semtner and Chervin, 1992; Stammer et 
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al., 1996). Temporal variability dominates the errors of inversions done using synop¬ 
tic section data. An important example is the large estimated potential deviation of 
the silica flux from its long-term average value when computed from any particular 
section. 

3. The use of conservation equations in terms of property anomalies (e.g., McDougall, 
1991a; Wijffels, 1993) rather than full property values themselves-but based upon 
making explicit the necessary assumptions concerning the error covariance between 
mass constraints and property constraints (Wunsch, 1996, p.273). 

4. Use of neutral surfaces (McDougall, 1987; Jackett and McDougall, 1997) to define 
layer boundaries with calculation, as part of the inversion, of the dianeutral transfers, 
dianeutral diffusion of tracers, freshwater fluxes along with the reference velocities 
and ageostrophic transport adjustment in the upper-most layer (that I ascribe to 
the Ekman transport uncertainty). 

o. The Ekman transport in the surface layers is derived from the mean wind fields from 
the NCEP reanalysis (Kalnay et al. 1996). 

6 . The nutrient and oxygen fluxes and divergences are analyzed in relation to biogeo¬ 
chemical cycles. 

The major part of the thesis concerns the estimation of a globally consistent large-scale 
oceanic circulation. In the Atlantic Ocean, my estimated water mass transports during the 
WOCE period (1985-1996 globally, although Atlantic sections span 1991-1993) are similar 
to previously published estimates. However, differences are found, noticeably in the Pacific 
Ocean, where the inflow of bottom waters is smaller than most previous estimates. 

By imposing stringent heat and salt conservation constraints, property anomalies al¬ 
low the estimation of consistent and significant dianeutral transfers and diffusivities, in 
contrast with many previous box inversions. Deep dianeutral advection and diffusivi¬ 
ties that are required by the tracer fluxes (mainly heat and salt) average respectively 
0.1±0.04x10 4 cms~ 1 and 3±1 cmV 1 . Larger diffusivities of about 20 cm 2 s _1 are found 
in bottom layers. 
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Energy transports evidence significant ocean-atmosphere heat fluxes, with a global 
heating of 2.3 ± 0.4 PW in the tropical band. Most of the cooling occurs in the northern 
hemisphere (1.7 ±0.2 PW) while in the southern hemisphere only 0.7 ±0.3PW is restored 
to the atmosphere. Nutrient residuals are in general uncertain. There is an interesting, 
but weak indication in several regions of nutrient consumption in the surface layers and 
deeper remineralization. The large uncertainties are due to time variability in the measured 
temperature and salinities and to the limited vertical resolution of the model. Nevertheless, 
the results provide an upper bound on the average remineralization rate in several depth 
ranges. Oxygen residuals indicate large and significant exchange with the atmosphere, 
with oxygen intake at high latitudes and outgassing at low latitudes. Significant oxygen 
utilization rate estimates indicate remineralization in subsurface layers in several regions. 

The remaining uncertainty in estimating the global oceanic circulation and fluxes is 
dominated by oceanic variability in the density and other properties. In particular, the 
large variability of the Pacific-Indonesian Throughflow degrades the property budgets in 
both the Indian and Pacific Oceans. Similarly, the variable North Brazil Current impacts 
the Atlantic’s budget. 

The thesis is organized as follows. Firstly, the method used to estimate the circulation 
and transports is described (Chapter 1). A detailed analysis of the uncertainties due 
to oceanic variability and measurement noise follows in Chapter 2. My estimate of the 
physical oceanic circulation for each oceanic basin is described in Chapter 3. The data 
treatment and model constraints that are specific to each ocean are presented there. In 
Chapter 4, I analyze the large scale fluxes of freshwater, heat, nutrients and oxygen in a 
global context. Chapter 5 gives a general summary of the thesis and a discussion with 
conclusive remarks and perspectives for the future. 
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Chapter 1 


Data and Model 

1.1 The data 

1.1.1 Hydrographic sections 

The circulation estimate is mostly based upon World Ocean Circulation Experiment 
(WOCE) hydrographic sections with high-quality water sample data. The sections that 
define the model are displayed on Figure 1-1 and listed on Table 1.1. In the Atlantic at 
36°N, we attempted to use the older AT109 data (Roemmich and Wunsch, 1985) from 
1981 because the 1993 data made available to us (Sokov et al., 1998) showed a large bias 
m the nutrients and oxygen. However, the 1981 data turned out to be incompatible with 
the WOCE lines A5 and A2, producing large and significant residuals in the upper layers, 
and it was decided to use neither section at this latitude. The Pacific section P4 (KPN) 
was also removed from the model because it caused large mass residuals and downwelling 
at intermediate levels, probably because of baroclinic variability in the Mindanao Current 
that was aliased in this section (S. Wijffels, personal communication, 1999). 

Given the data, the model solution, at best, represents an average of the circulation for 
the 1989-1996 period in most oceans, except in the North Pacific where the sections were 
taken earlier, in 1985. Many sections are not public or were not public when I started this 
work. They were used with the consent of the principal investigators and I thank them for 
letting me use their data, either by making them publicly available or by providing them 
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Section 

Date 

Reference 

A2 Gauss 226 

Jul. 93 

Koltermann et al. (1999) 

AT109-36 

Jun. 81 

Roemmich and Wunsch (1985) 

A5 

Jul. 92 

Parilla et al. (1994) 

A6 

Feb. 93 

Arhan et al. (1998) 

Oudot et al. (1998) 

A7 

Jan. 93 

A8 

Mar. 94 

Siedler et al. (1996) 

A9 

Feb. 91 

A10 

Dec. 93 

All 

Dec. 93 

Saunders and King (1995) 

A21 

Jan. 90 

Roether et al. (1993) 

A12 

May 92 

Lemke (1994) 

16 

Feb. 96 

Park et al. (1999) 

15 

Nov. 87 

Toole and Warren (1993) 

14 

Jun. 95 


13 

Jun. 95 


12 

Dec. 95 

Johnson et al. (1998) 

110 

Nov. 95 

Bray et al. (1997) 

JADE 1989 

Aug. 89 

Fieux et ah, 1994, 1996 

I9S 

Dec. 95 

Hufford et al. (1997) 

P12 

Jan. 95 

Rosenberg et al. (1997) 

P14S 

Jan. 96 

McTaggart and G. Johnson (1997) 

P6 

May 92 

Tsimplis et al. (1998); Wijffels et al. (1999) 

P21 

May 94 

Banks et al. (1995) 

P4 

Mar. 89 

Wijffels et al. (1996) 

P3 

May 85 

Roemmich et al. (1991) 

PI 

Apr. 85 

Talley et al. (1991) 


Table 1.1: Hydrographic sections 


directly to me. Their respective names are listed in the physical description of individual 
oceans in Chapter 3 and in the acknowledgments. More details on the data are provided 
in Chapter 3. 


1.2 Model formulation 

Inverse methods applied to the determination of oceanic circulation are now widely doc¬ 
umented in the literature. I used here the so called “inverse box model” described by 
Wunsch (1978) and subsequently used in various studies (e.g., Wunsch, 1982; Rintoul 
and Wunsch, 1991; Wijffels, 1993; Macdonald and Wunsch, 1996; Tsimplis et al., 1998, 


16 























































Top-to-Bottom 

within individual layers 

Mass 

yes 

yes 

Thermodynamic energy (heat) 

no 

below surface 

Salinity 

yes 

yes 

Silica 

yes 

no 

“NO” ( 9 .INO 3 +O 2 ) 

no 

no 

“PO” (I 7 OPO 4 +O 2 ) 

no 

below surface 


Table 1.2: Conservation constraints for the standard solution. Heat, salinity and “PO” 
use the anomaly formulation. 

etc.). Wunsch (1996) gives a complete description of the various methods along with their 
statistical background. 

The geostrophic circulation with a surface Ekman layer is required to conserve mass 
and other tracers within individual layers bounded by neutral surfaces (McDougall, 1987). 
Table 1.2 gives a summary of the different constraints and how they are applied. Con¬ 
servation constraints take into account advective and diffusive exchanges between layers, 
hereafter referred to as dianeutral fluxes. 

The thermodynamic energy (usually mislabeled “heat”, e.g., Warren (1999)) and the 
phosphate-oxygen “PO” (I 7 OPO 4 +O 2 , Anderson and Sarmiento, 1994) combination are 
approximately conserved in layers that are not in contact with the atmosphere. (The sur¬ 
face layer is non-conservative because of heat and oxygen exchanges with the atmosphere.) 
The value 170 for the PO Redfield ratio is controversial (e.g., Minster and Boulahdid, 1987 
; Anderson and Sarmiento, 1994). However, an experiment with a lower ratio of 140 showed 
that the circulation was not changed significantly. Conservation equations for the tracer 
“NO” ( 9 .IXNO 3 +O 2 ; Broecker (1974)) are, in general, indistinguishable from the “PO” 
equations and I chose to conserve “PO” alone because the nitrogen cycle is the more com¬ 
plex one. The total (top-to-bottom) silica is approximately conserved, but silica residuals 
are allowed in individual layers to accommodate consumption and dissolution through the 
biological cycle (Section4.3.1). 

Bulk advective (w*) and diffusive (n*) dianeutral transfers are allowed between the 
layers for each property and (Section 1.2.2). The same w* and k* are used for all properties. 
By definition there is no diffusion in the mass equation (e.g., Pedlosky, 1987, page 10). (It 
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is also possible to derive a density equation from the energy conservation equation, but 
this equation is redundant with heat and salt conservation). The bulk average property 
and property gradients on a neutral surface are calculated from the bounding sections. 
Freshwater flux divergences and adjustments to the Ekman transports in the surface layers 
are solved for as well (Section 1.2.3). 

To better condition the system, property anomaly equations are used for salt, heat and 
PO (Section 1.2.4) but not for silica, i.e. Section 1.2.4. Solutions are obtained using the 
Gauss-Markov estimator so that the uncertainties include the null-space (the part of the 
solution that is not resolved by the constraints), in contrast with standard least-squares 
and singular value decomposition methods (Section 1.3). All conservation equations were 
weighted according to their a priori uncertainty that will be determined in the next chapter 
(Chapter 2). 

1.2.1 Flux through a single hydrographic section 

Consider a hydrographic section running in the zonal (x) direction, composed of a set of 
stations, where temperature (T), salinity ( S ), nutrients and oxygen were measured over 
the water column. As the circulation is in geostrophic balance to a good approximation, 
the velocity field is obtained by integrating the thermal wind equation, 

v abs(^, z) = V re l(x, z) + b(x) = —f ~^-dz + b(x), (1.1) 

PofJ Zo ( x )dx v ’ 

where z is the vertical coordinate; x is the horizontal coordinate along the section; zq[x) is 

an arbitrary reference depth at each station pair; v a i, s (x, z) is the absolute velocity , the sum 

of the velocity relative to the reference depth, v re i{x,z) computed from the data, plus the 

reference velocity, b(x), which is unknown 1 , g is the local gravitational acceleration and 

/ is the Coriolis parameter. The surface formed by zq(x ) is called the initial zero-velocity 

surface, or reference surface. Depending upon the section, this surface may be a straight 

line (reference level), a neutral surface, or may be specified for each pair according to 

the surrounding water properties. Both b and v are normal to the section and defined as 

the horizontal average between two consecutive stations (a station pair). Classically, the 

: I have attempted to keep my notations consistent with Wunsch (1996). 
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circulation was inferred by elaborate guesses of the reference surface, zo(x), at which b(x) 
would vanish. The inverse box model is cast as providing an estimate of the adjustment, 
{6} = b, to this classical a priori solution, along with an error covariance P, which 
satisfies a number of requirements on the circulation and the a priori variance of b. These 
requirements may include, for instance, mass, salt and silica conservation. 

The net transport of a tracer across the section is, by definition: 

/■East rTop{x) 

T= / dx dz [b(x) + v re i(x,z)] p(x,z) C (x,z) (1.2) 

JWest J Bot(x) 

= T b + r (i.3) 

where C is the tracer concentration per unit mass ((7 = 1 for mass), T b is the contribution 
from the reference velocities and T is the contribution from the relative velocities. In the 
discrete form, T b is a linear combination of the b's, T b = a T b, where the elements of a 
contain the integrals of pC over the water column at each station pair and b is the column 
vector containing b(x). T may be the transport through an individual layer, Top(x) and 
Bot(x ) then defining the upper and lower boundaries of the layer. 

1.2.2 Conservation equations 

In steady-state, mass conservation (continuity) may be written as 

V• (pv) = 0 

while a tracer conservation equation may be written as 

V - (pCv) = V • k V (pC) , (1.4) 

where C is the tracer concentration; v is the three-dimensional velocity vector, and k is the 
eddy diffusion tensor. In the present model, it is assumed that the spatial sampling along 
a section is fine enough to resolve most of the horizontal tracer transport due to eddies at 
a given time (the ergodicity assumption is discussed in the next chapter). The smaller- 
scale lateral processes (including horizontal diffusion) are contained in the noise. One can 
integrate (1.4) over a layer bounded horizontally by two zonal hydrographic sections and 
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land, and vertically between selected neutral surfaces and rearrange the terms to obtain 


/ Bot(x) 


dz pCv 


-[II 


dxdy I k 


d( P C) 


— wpC 


where [ X= X North — x South , w is the velocity perpendicular to the top or bottom 
interface and k the dianeutral diffusivity. The left-hand-side can be calculated from the 
velocities across the sections, but the right-hand-side is not directly computable. The 
right-hand-side is decomposed as follows, assuming that the second order Reynolds terms 
are negligible: 


JJ dxdy 


d(pc) 


— wpC ) ~ A x k 


<d( P C) 


w*pC\, 


where the horizontal overbar denotes a horizontal average; A is the horizontal area of a 
given neutral surface bounded by the sections, and w*, k* represent the equivalent average 
dianeutral velocity and diffusivity across the interface. Such a formalism assumes for 
«* > 0 down-gradient turbulent fluxes in the dianeutral direction. Both pC and ^§p- 
are calculated from the section data, that is, at the boundaries of each ocean sector. 
The horizontal area between sections, A, was calculated according to the neutral surface 
positions as derived from the Levitus (1994) climatology. 

Using (1.3) and (1.6) and accounting for possible sources and sinks of each tracer, 
Eq. (1.5) may be written as 


l T b] Smith + A X W*pC - K 


■ d(pC) 
dz 


= [-r ]S£5 + Q, 


where Q is a possible source or sink of tracer C. (Q includes the freshwater flux in the 
mass budget and the Ekman flux divergence for all tracers including mass.) Eq. (1.7) is 
written in the discrete form as 

a T b + n = y, (1.8) 


} i=l,Npair 

with b = 

k=l,Nlay—2 
{ K k) k=l,Nlay-2 

A FW 
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b is the column vector containing the unknowns of the system, respectively, the reference 
velocities at each station pair {&*}; the adjustment to the Ekman transport (see below); 
the dianeutral fluxes across each interface {w^} and {«£}; and the freshwater flux in the 
mass equation, A FW. The term y contains the relative transport divergence (T) and 
the source term Q. A noise n has been introduced in Eq.1.8 as the equation cannot be 
satisfied exactly with real data. A quantitative discussion of the noise follows in Chapter 2. 
A similar equation is added for each layer, for the total flow, and for several tracers 
including mass (C = 1). All the equations of the form (1.8) may be grouped to form a 
single matrix equation 

Ab + n = y, (1.9) 

where A is now a matrix of size M x N, where M is the number of equations and N is 
the number of unknowns, (n and y are M x 1.) Constraints on the net flux through a 
particular set of station pairs and layers (e.g., Florida Strait flux, topographic constraints) 
are added by choosing the appropriate elements of a and y. 

1.2.3 Ekman transport 

The Ekman transport is included in (1.9). Annual variations in the Ekman transport can 
be large (e.g., Levitus, 1988). However, the relevant quantity to our time-average circu¬ 
lation model is the time-average Ekman transport. The Ekman transport was obtained 
from an average wind stress of the NCEP/NCAR reanalysis (Kalnay et al., 1996) over the 
period 1992 to 1997 2 . In the model formulation, the Ekman transport is decomposed into 
two contributions: 


T & e = T™ t + AT Ek , (1.10) 

where Tjjjj?* is the initial Ekman transport, calculated from the climatology, and AT#*, is 
the correction required by the conservation equations. A Tsk is included in b and solved 
for in the inversion. The a priori uncertainty on the Ekman transport, AT^, is 50% of the 

2 This period does not exactly overlap with the hydrographic data. The model-derived correction to the 
Ekman transport is assumed to accommodate the temporal differences. 
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initial value, with a lower bound of 1 Sv (see the discussion in Chapter 2, Section 2.3.4). 
The value for initial Ekman transports T 1 ^ 1 and for the correction after inversion A T Ek 
are given for each ocean in Chapter 3 (Tables 3.2, page 78; 3.6, page 95); 3.11, and page 112. 

The Ekman transport of tracers is estimated with the product ( T Ek + AT Ek ) xC where 
C is the average property concentration in the upper layer(s) on each full section. In the 
places where a layer outcrops along a section, the Ekman transport was calculated over 
the interval of outcrop. But C was taken as the average over the whole layer. This simple 
calculation excludes the effects of correlations between horizontal property variations and 
Ekman transport within the surface layer. It is however a good approximation for the net 
temperature transport which is mainly controlled by the vertical temperature difference 
between the surface and deep ocean, horizontal gradients playing a minor role (S. Jayne 
and J.Marotzke, personal communication, 1998). Equivalently, the nutrient fluxes are 
not much affected by the approximation because nutrients are generally depleted in the 
upper layers at low latitudes, where the Ekman transport is strong. The effect of the 
approximation on oxygen was not estimated; it is assumed to be included in the noise. 


1.2.4 Anomaly Equations 


The noise in the mass conservation equations tends to dominate most tracer conservation 
equations. As a result, the original full property equations produce an ill-conditioned 
matrix (e.g., McDougall, 1989; Wijffels, 1993). To better condition the system, one can 
subtract the mass divergence times the average property concentration for each layer from 
each tracer conservation equation. As discussed by Wunsch(1996) , the main issue, when 
using anomaly equations, is the determination of the degree of cancellation in the residual 
noise after subtraction. 

We now estimate this degree of cancellation, or equivalently the size of the noise term 
in anomaly equations. The salt equation is taken for illustration. For a single zonal 
hydrographic section, and within a layer, the net transports of mass and salt are written 


Tm = JJ dx dz (pv ( x , z) + n v ( x, z)) 
and Ts = JJ dx dz (pv ( x , z) S ( x , z) + n vs (x, z )), 


( 1 . 11 ) 

( 1 . 12 ) 
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where v(x,z) is the meridional velocity; p is density; the terms n v (x,z ) and n vs (x,z) 
refer to noise in the advective transports of mass and salt at individual station pairs; x 
is longitude and z is the vertical coordinate. Let S be the mean salinity of the layer and 
S' = S~S the salinity anomaly. Then by subtraction (1.12) — S x (1.11) one obtains an 
expression for the salinity anomaly transport: 


T s ’ = j j pv ( x , z) S (x, z) dxdz + n' 5 , 

(1.13) 

with n' 5 = JJ ( n vs ( x,z) — Sn v {x,z))dxdz. 

(1.14) 


A mass transport error (the subscript m stands for mass) is anticipated to produce 
a salt transport ns of roughly S x riM such that should be much smaller than ns- 
We now formalize this statement. Assuming that the noise in the salt flux, n vs (x, z ), 
is mainly due to noise in the mass flux n v ( x , z) (as opposed to being mainly caused by 
variations in the salinity itself), n vs (x,z) may be written as 


n vs (x, z) = n v (x, z) x S (x, z ), 

(1.15) 

or n vs (x, z) = n v (x, z) x (5 + S' (x, z)) , 

(1.16) 

so that n' s = JJ n v ( x , z) S' (x, z) dxdz. 

(1.17) 


may be divided into a component owing to the zonally averaged mass transport over 
the layer and a deviation that we call the horizontal eddy component: 


n's = n M // S' (x, z ) dxdz + JJ ( n v (x , z) — n m) S'dxdz, 
with n m = J J n v (x,z) dxdz. 


(1.18) 

(1.19) 


nM is the noise in the zonally integrated mass flux over a layer. S is defined over a 
whole layer, possibly bounded by several sections, so that the first term of (1.18) (the 
zonal average) is in general non-zero. As an approximate estimate for the variance of 
the zonally averaged component, we chose to use (n|^) (A5 2 ), where AS is the a priori 
size of salinity variations within the layer. The size of the horizontal eddy component 
is unknown a priori, and is assumed to be of the same magnitude or smaller than the 
zonally averaged component with possible correlation between the two. This speculative 
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assumption is supported by studies of the heat (energy) flux variability (salt was used 
as an example, but similar anomaly equations are written for other tracers including 
internal energy). Jayne (1999) found that at all latitudes, variations in the baroclinic heat 
transport-which appears as a noise, ng, in the heat equations-are about twice as large 
as variations in the horizontal eddy contribution to the heat transports. Instantaneously, 
salinity anomaly transports display a similar behavior as heat anomaly transport (e.g., 
Wunsch, 1996, Fig. 4-15) so that the noise in the salinity anomaly may be reasonably 
assumed to be similar. 

Thus, the ad hoc, a priori, noise in the tracer anomaly equations for any tracer C is 

( n c) = 4 x (n 2 M ) x (AC 2 ) , (1.20) 

where AC is the standard deviation of tracer variations within the layer and the factor 
of 4 accounts for possible correlations between the zonally averaged and horizontal eddy 
components. The average tracer concentration ( C ) in each layer was used in defining the 
anomaly equations. Expression (1.20) relies on the assumptions listed above. It is no more 
than a best guess, which turns out to produce consistent results (in terms of error bars) 
when used in the inversion. The main test of consistency comes from the fact that the 
residuals of all conservation equation are indistinguishable from zero within a posteriori 
error bars. 

Wijffels (1993) weighted her anomaly equations in a similar way, but assumed that the 
salinity variations were dominated by vertical differences. Because we have no information 
on the noise in vertical exchanges, it is not clear whether this approach is more appropriate 
than ours. Nevertheless, vertical salinity differences are on average of the same order as 
the horizontal variations within a layer, but display more irregularities, which make them 
difficult to use as a reference. 

Silica and anomaly equations 

Anomaly equations are inappropriate for top-to-bottom silica conservation. Generally, 
most of the spurious mass residuals in layers occur close to the surface in part because of 
the larger variability in the thermal wind shear there. Because the average concentration 
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over the water column ( C ) is much higher than the concentration in the surface layers 
where the mass divergence actually occurs, the term Cum (or S for salinity) of (1.14) is 
artificially large, creating an erroneous source or sink of silica in the top-to-bottom anomaly 
equation. Hence, using the surface concentration for creating the top-to-bottom anomaly 
equations or, for simplicity and consistency, using non-anomaly (original) equations for 
total silica conservation is more appropriate. 

Diagnostic of tracer divergences 

Tracer divergences were diagnosed following the anomaly formulation, effectively reducing 
the impact of noisy mass divergences from the tracer divergences. Doing so in each indi¬ 
vidual layer i by using the average property Ci yields a corrected tracer divergence XJ.C l . 
But if Ci varies substantially with depth (as for temperature or nutrient anomaly), the 
corrected tracer divergence in the top-to-bottom integral (V-Ct 2 s) ' s n °t the same as the 
sum of the anomaly layer divergences: 

V-Ct 2B 7^ 53 V-Cj- (1-21) 

layers 

Intuitively, using \J.C 1 in individual layers is more reliable than the top-to-bottom cor¬ 
rected divergence because the mass divergences are removed more efficiently (see the dis¬ 
cussion about silica just above). Still, if the difference (1.21) is significant within error bars, 
it will be interpreted as a warning that the a priori uncertainties (1.20) are inconsistent. 

1.3 The Gauss-Markov estimator 

The solution estimate b that satisfies (1.9) is obtained using the statistical Gauss-Markov 
method (i.e. Wunsch, 1996, page 184). We explicitly use the a priori knowledge of the size 
of the solution b, R-b,b = (t>b r ), and of the priori noise level variance, R n , n = (nn T ) (the 
() operator stands for the ensemble average). In the absence of a priori knowledge of the 
correlation between different equations and solution elements, Rb,b and R n ,n are taken 
to be diagonal. Determination of the a priori noise in the conservation equations is the 
focus of Chapter 2. The a priori size of the solution b is given in Section 2.4 of Chapter 2 
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and for individual basins in Chapter 3. 

Given Rb,b and R nn , the estimate b that minimizes the deviation from the hypothet¬ 
ical true value is given by: 

b = Ky (1.22a) 

with K = R b ,bA T (AR b , b A T + Rn,n) _1 (1.22b) 

If there is no bias (^b - bj} = 0), that is, both (y) = 0 and (n) = 0, (1.22a) is the Best 
Linear Unbiased Estimator (BLUE). 

The error covariance is 

p = (b - b Tru^j (b - b True ) ^ 

= R b ,b — K A R b , b (1.23) 

The error estimate includes the part of the solution which is unresolved by the equations 
(the null-space). The residuals and uncertainties are computed from (1.9), 

h = y - A b 
P„ n = APA t . 

For consistency, it is checked that the residuals, after inversion, are indistinguishable from 
zero within error bars, which is a necessary condition for obtaining realistic error bars. (It 
follows that the residuals are within the a priori range R n n ). 

Estimates for property transports (t) and uncertainties (sf^j are derived from (1.8) 
and (1.22a,b-1.23): 

f = r re/ + a T b (1.24) 

5T = V a T Pa, (1.25) 

where T re i is again the property transport due to the relative velocity and a T is a row 
vector containing the integrated property concentration at each station pair. 
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Chapter 2 


The Error Budget of Inverse Box 
Models 


In this chapter I estimate the uncertainty in mass and silica transport estimates due to 
the non-synopticity of hydrographic sections (i.e., the departure of nearly-synoptic fields 
from true time-means), numerical errors in transport calculations and the measurement 
noise. In particular, a an eddy-permitting general circulation model (GCM) is used to 
estimate the extent to which a hydrographic cruise is representative of the true average 
density field and associated transports. This study is relatively extensive, and the reader 
may chose to go directly to the summary of a priori uncertainties that will be used to 
weight the conservation equations, Section2.5 (page62) and Section2.6 (page66). 

2.1 Introduction 

As we have seen, the linear inverse box model can be condensed into one equation: 

Ab + n = y, (2.1) 

where A is a matrix, b is a vector containing all the unknowns of the system, y is the 
observation vector and n is the noise. Each row of (2.1) is an equation (or constraint) 
on b and, given some reasonable assumptions about the variance of b and n, the in¬ 
verse machinery provides an estimate for b and its covariance. In the present analysis, b 
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contains the unknown reference velocities and the product Ab represents the transport 
and transport divergence of different properties (mass, heat, nutrients...) due to the ref¬ 
erence level velocities, y contains the transport from thermal wind, Ekman transport and 
other possible sources and sinks in the conservation equations. The expected values of the 
solution and of the noise are assumed to be identically null, that is, < b >=< n >= 0. 

The uncertainty in the solution b and in the associated fluxes depends upon: 1) the 
weight of the constraints (i.e., the a priori variance of n, denoted (nn T )) and, 2) the 
a priori variance of b (denoted (bb T )). In most previous studies, these two a priori 
quantities were associated with the error in the net mass flux that could enter or leave an 
oceanic area (freshwater flux, mass storage). For instance, Macdonald (1995) and Rintoul 
and Wunsch(1991) required mass conservation using an a priori error of ±1 to ±2Sv-an 
optimistic value, as will be suggested in this chapter. 

In a related work, McIntosh and Rintoul (1997) tested the ability of an inverse box 
model to reproduce an instantaneous snapshot of the southern ocean circulation as com¬ 
puted with a quasi-eddy resolving general circulation model (GCM). They found the box 
model successful as long as the prior statistics on the solution variance (i.e., (bb T ) in our 
notations) were chosen to correspond to the true (GCM) solution. The object of this 
chapter is to review and quantify the different sources of error entering the conservation 
equations, and to provide estimates for the prior statistics (nn r ) and (bb r ). 

Mass Flux Estimation 

At a given time t, and assuming for the moment that geostrophy is a good approximation 
away from the surface layer, the mass transport through a transoceanic hydrographic 
section, T(t), is the sum 


T(t) = T R (t) + T b (t) + T Ek (f), 

(2.2) 

with T r (t) = - JJ dxdzx 

(2.3) 

Tb ( t ) — JJ dxdz x pb(x ); 

(2.4) 

r x E T (f \ x 

and T Ek (£) = -/ dx x p . 

Jx w J 

(2.5) 
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Tfi ( t ) is the thermal wind transport relative to a chosen “initial zero-velocity-surface” 
(IZVS), or reference surface. T b ( t ) is the transport contribution from the reference velocity 
(the unknown) and T Ek ( t ) is the near-surface Ekman transport, p is density, g is gravity, 
r ( t) x the zonal wind stress; / is the Coriolis parameter and t is time. 

The data I use span long time intervals: it takes one to three months to make a 
transoceanic section; and two different sections used in a box model may be separated by 
several years (Table 1.1, page 16). The object of the inverse model is to estimate the true 
time mean transports across all the sections involved. Write this average 

f = f R + T b + T Ek , (2.6) 

where the overbar denotes a time average. Because Tr is estimated from the measurements 
at a single time to, the existence of a solution that is mass conserving relies on the as¬ 
sumption that the large scale thermal wind transport at time to, Tr (to) , is representative 
of the time-average Tr to a good approximation. 

It is useful to separate the measurement noise, n measurement , and the model error, 
n modei_ r pp c measurement noise comes from uncertainties in the measurement of Tr (to) 
(internal wave noise, eddy aliasing, bottom triangle, etc.) while the (inverse) model error 
represents the limitation to a barotropic adjustment, permitting no change in the vertical 
shear-thus in Tr. It follows that the model error corresponds to the difference between 
Tr (to) at the time of each section, and the time-average Tr. This relationship may be 
written by 

rj-i true _ /^measured [Iq'j _f- ^measurement 

Tr = Tjf ue (to) + n model . (2.7) 

As we will see, the error from the oceanic baroclinic variability, n model , dominates the total 
error. [In a more elaborate inverse model like the ones of Mercier (1986) or Wunsch (1994), 
change is allowed in the density field, thus reducing the model error. Because no such 
change is allowed here, any noise in Tr may lead to a net mass imbalance in the solution.] 
The sizes of n model and n m easurement are estimated at three latitudes in the Atlantic: 
24°N, 36°N and 48°N; n model is estimated with the help of a near-eddy-resolving general 
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circulation model (GCM), while the measurement errors are estimated from data and 
experiments. The model error is investigated in Section 2.2. There it is shown that the 
oceanic response to Ekman transport fluctuations is essentially a balance between Teh (t) 
and T(, ( t ) (Section 2.2.2) so that the variability in the relative transport is relatively small 
(but not negligible, Section 2.2.3). Along with the error budget in the mass flux, the model 
error in the silica equations is explored in Section 2.2.4. Measurement errors are quantified 
in Section 2.3; those include internal wave noise (Section 2.3.1); geostrophy (Section 2.3.2); 
bottom triangle approximation (Section 2.3.3) and mean Ekman transport (Section 2.3.4). 
The reference velocity variance in the GCM is quantified in Section 2.4. In Section 2.5, 
all components of the uncertainty are summarized and discussed. The application to our 
global hydrographic inverse model is described in Section 2.6. 

2.2 Model Error 

2.2.1 The general circulation model (GCM) 

The GCM that I used to test the representativity of one-time hydrographic sections is a |° 
lateral resolution global circulation model run at the US Naval Post Graduate School (see 
Semtner and Chervin, 1992; Stammer et ah, 1996). The model has 20 vertical levels and 
was forced with the daily wind stress fields from the European Center for Medium range 
Weather Forecast (ECMWF). The surface temperature and salinity were restored to the 
monthly fields of Levitus et al. (1994) with a 30-day time scale [see Stammer et al. (1996) 
for more details]. Stammer et ah (1996) compared the model output with both altimetric 
and hydrographic observations and found that the GCM reproduced successfully the major 
features of the large scale oceanic circulation and the structure of the variability. However, 
the GCM underestimated the amplitude of the variability by at least factor of 2 to 4 on 
most scales. The discrepancy was highest in the mesoscale range. The model run, labeled 
POCM4-B, was extracted by S. Jayne (1999) and made available to me for the purpose of 
this study. Temperature, salinity, velocity and sea surface height data were subsampled 
every third day for a period of 8 years (1987 to 1994) and at three latitudes corresponding 
to actual hydrographic data across the Atlantic: 24°N, 36°N, and 48°N. 
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2.2.2 Oceanic response to Ekman fluctuations 

Areal integration of the GCM sea surface height between 36°N and 48°N revealed that 
in the GCM, the rate of change in mass storage did not exceed ±0.2 Sv between the two 
latitudes. Similarly, altimetric observations of the ocean do not show large variations in 
the mean sea level (Wunsch, personal communication, 1997), suggesting that the Ekman 
mass inflow is almost instantaneously compensated by an interior flow in the opposite 
direction (e.g., Jayne, 1999). (On seasonal time scales there are variations from steric 
height change, though.) 

Boning and Hermann (1994) studied the structure of the compensating interior flow 
and found it depth-independent at mid-latitudes in a near eddy-resolving model of the 
North and Equatorial Atlantic. An integration of the present GCM transport in different 
layers confirmed that the temporal variation in the shear was smaller and slower than 
the depth-independent part. This barotropic response suggests that the thermal wind is 
stable at mid-latitudes, as will be confirmed in Section 2.2.3: most of the fluctuations in 
the Ekman transport are represented by (t) (Eq.2.2). At low latitudes, baroclinic waves 
can propagate faster than at mid and high latitudes, and the oceanic response is different 
(Boning and Herrmann, 1994). The study of the low latitude response to variable wind 
forcing is postponed to a possible extension of this work. 

The Ekman transport was computed across each section for one year in the Atlantic 
Ocean from the ECMWF wind stress fields (Fig.2-1, thick line) using the classical Ekman 
formula (2.5). The Ekman transports change quickly with variations of up to 30 Sv in a few 
days at all latitudes. (The accuracy of the ECMWF wind stress themselves are discussed 
later in Section 2.3.4.) The same winds were used to compute the Ekman transports as 
were used to force the GCM. 

On Figure 2-1 is overlain the Ekman transport and the total flow beneath 100 m in 
the GCM. The balance between the two curves is not exact because of energetic inertial 
oscillations. But if high frequencies are filtered, the imbalances do not exceed ±0.8 Sv rms 
(not shown), showing that Eq.(2.5) is appropriate to use in the GCM. (Indeed, the numer¬ 
ical parameterization of momentum transfers in the GCM upper layers is close to Eq.2.5.) 
The immediate consequence of the oscillations is a large fluctuation in the energy flux due 
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to the temperature contrast between the Ekman layer and the average temperature over 
the water column (Jayne, 1999). 

Because the Ekman transport fluctuations are almost balanced by a barotropic return 
flow, different hydrographic sections can be consistently combined in an inverse box model 
as long as the time-mean Ekman transport is used for all sections. Then, the calculated 
energy flux is representative of the time-average energy flux. 

2.2.3 Baroclinic variability 

I now turn to the error caused by baroclinic variability, that is, the deviations between the 
instantaneous transport Tr (t) and the average Tr. Locally, mesoscale eddies introduce 
important deviations from the mean transport, so that estimating the time-average ther¬ 
mal wind current would require an enormous data set (would be inappropriate to represent 
small-scale property transports). However, eddy effects tend to cancel out when integrated 
on large scales so that the zonally integrated thermal wind transport is relatively stable. 

The baroclinic variability in the ocean is a broad subject. In their theoretical study of 
the seasonal variability, Gill and Niller (1973) found that the ocean response (baroclinic 
and “barotropic”) to seasonal forcing was relatively small away from the Equator. Air-sea 
exchanges of heat are stored locally and changes in the horizontal velocities are weak. 
Even on seasonal time scales, the authors found that most of the response is barotropic 
(i.e., independent of the density structure) and moreover, that horizontal advection of 
temperature and salinity by barotropic currents does not produce significant baroclinic 
variations. But although barotropic current variations are locally small (±3mms _1 ), 
they produce transports of ±30 Sv when integrated on large scales. 

Variability in the upper North Atlantic on longer time scales has been investigated 
by several authors, including Reverdin et al. (1997) and Curry et al. (1998) using the 
few available long time series of temperature and salinity. Large variations occur in the 
temperature and salinity field on 10 year time-scales. Salinity anomalies propagate in the 
North Atlantic and influence the convection regimes and in turn, the production of deep 
water masses (Reverdin et al., 1997). For instance, the Labrador Sea Water outflow is, for 
a similar density, alternatively cold and fresh or warm and salty. In turn, changes in the 
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surface conditions are imprinted in the deep waters, and the corresponding signals can be 
observed in the deep Bermuda waters (Curry et ah, 1998). The variability relates to the 
North Atlantic Oscillation index and sea ice coverage (i.e., Reverdin et al., 1997; Curry et 
al., 1998). 

Based on a simple, baroclinic, linear Rossby-wave model, Sturges et al. (1998) report 
decadal variations in the zonal differences of sea surface height of ±4 cm in the Atlantic at 
34°N and ±10 cm at 18°N (away from the Gulf Stream). Assuming these variations were 
confined within the surface 1000 m, the baroclinic variation in the gyre strength would 
be respectively ±5 Sv at 34°N and ±20 Sv at 18°N. These variations are entirely due to 
changes in the wind stress curl, and their model produced a reasonable agreement with 
hydrographic data. However, the variations in the net, coast-to-coast flux, are probably 
much smaller than the ones they reported because the western boundary current, which 
they did not include, probably balances the variations in the gyre strength. 

Greatbach et al. (1991) obtained variations of similar amplitude -but with the opposite 
sign-using a linear diagnostic model of the North Atlantic based on climatological wind 
and density data. Comparison of the mean circulation for the period 1955-1959 with that 
of 1970-1974 showed a 30 Sv decrease in the Gulf Stream gyre. Differences in the sampling 
between the two pentads are a large source of uncertainty, however, half of the variations 
were explained by changes in the density field above 1500 m following a weakening of the 
subtropical gyre. The other half is uncertain because of the large time and space sampling 
differences in the deep waters between the two time periods considered. In other oceanic 
regions, similar variations may be associated with ENSO and other large scale atmospheric 
patterns. 

Repeat hydrographic sections are appropriate to study the long term interannual vari¬ 
ability in the net fluxes. At 24°N and 36°N, sections were done during IGY (1957-1959), 
then in 1981 and recently in 1992-1993. Evidence of time shifts in the water mass charac¬ 
teristics have been found (e.g., Bryden et al., 1996, Arbic and Owens, 1999). Roemmich 
and Wunsch(1985) found the geostrophic transports (thermal-wind±barotropic adjust¬ 
ment) from the two first occupations to be qualitatively similar, but the vertical profiles 
of meridional transport exhibited noticeable differences in their structure and water mass 
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fluxes. Their differences are compared with my results at the end of this section. 

In this section, I use the GCM output to quantify the thermal-wind variability. To 
remove the large variability in the GCM associated with inertial oscillations (to first order 
unrelated to the thermal wind, Section2.3.2), all calculations were based on the absolute 
geostrophic velocities from the pressure field. That is, from the thermal wind referenced 
at the surface and the sea surface elevation (r?): 


Vga{x,z) 


9 drjjx) g f z dp{x,z) 

f dx pof Jo dx 


(2.8) 


where x is the zonal distance, z the depth, g the gravity, / the Coriolis parameter and 
p the density. The Rossby number in the model did not exceed e = 0.2. Near the 
Gulf-Stream the effect of horizontal friction in the GCM is relatively small so that, away 
from the Ekman layer and except for inertial oscillations, geostrophy is expected to first 
order. Still, a large ageostrophic transport was found in the GCM output. It occurred 
on the grid scale near the Gulf Stream, but neither the non-linear terms nor the friction 
was able to support the 5Sv deviation observed (50% of the local transport). Several 
experiments led to the conclusion that the ageostrophy was due to a numerical artifact, 
possibly the bottom drag or an implicit diffusion acting on small scales (Appendix A). As 
the ageostrophic flow is depth-independent the present analysis of baroclinic variability 
was not directly affected. 

Because the GCM was integrated for 8 years (1987 to 1994) but forced seasonally 
with climatological temperature, salinity and heat fluxes (Stammer et al., 1996), the only 
interannual patterns come from the variations in the wind stress. Because of the repetitive 
seasonal forcing, the variability in the GCM thermohaline circulation is underestimated. 
Because the GCM underpredicts the oceanic variability (page 30), the variability over 
those 8 years provides a lower bound on the actual variability. (On the other hand, model 
drift, which was not removed, introduces additional spurious variability.) 

I start by investigating the seasonal and interannual variability in Tr assuming all 


hydrographic sections are measured instantaneously. Then, a one-month-long section is 
simulated in the GCM to assess the effect of the finite time required to sample the ocean. 
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Instantaneous section 


Following the dynamic method, the velocity was split into a reference velocity (denoted 
hereafter RV or b(x)) plus a relative velocity vr(x,z) so that 


Vga(x,z) = V R (x,z) + b(x). (2.9) 

The velocity is referenced to a deep isopycnal surface, as in the hydrographic model of 
Macdonald (1995): o 4 = 45.81 at 25°N and 36°N, and cr 2 = 36.87 at 48°N. At 25°N and 
36°N the reference isopycnals are about 3000 m deep, while at 48°N the reference surface 
lies between 1500 and 2000 m. 

The absolute and relative geostrophic transports were integrated in individual layers 
similar to the ones used by previous authors (Macdonald , 1995; Rintoul and Wunsch ,1991; 
Roemmich and Wunsch, 1985) and listed in Table 2.1. However, for simplicity, interface 
depths were computed without allowing changes in the reference pressure-in contrast to 
Macdonald’s definition 1 . Neutral surfaces (Jackett and McDougall, 1997, McDougall, 
1987) are more consistent with the Macdonald definition, but the code for computing 
neutral surfaces is too slow for computation at each time step. A snapshot of the interface 
positions is shown on Figure 2-2 at 36°N. 

The integration was done every third day for eight years. Figure 2-3 shows the total 
transport (Tr + T&) across latitudes 25°N, 36°N and 48°N along with the contribution 

'in the Macdonald definition a reference pressure is assigned to an isopycnal at one end of a hydrographic 
section. The isopycnal position is computed along the section. If the isopycnal deviates too much from the 
original reference, the local reference pressure is changed to continue the calculation along the section.. 


00 

02 

04 

26.4 

36.87 

45.81 

26.8 

36.94 

45.85 

27.10 

36.98 

45.87 

27.3 

37.02 

45.895 

27.5 


45.91 

27.7 


45.925 


Table 2.1: Potential density interfaces. Reference is taken at the surface, 2000 db and 
4000 db. 
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Figure 2-2: Isopycnal interface position along section A36°N 


from the relative velocity (Tr). The layer thicknesses and depths plotted correspond to an 
average across the section. Note the northward transports are displayed in each individual 
layer in contrast to the transport streamfunction that is used in the remainder of this 
thesis. 

In the deep layers, most of the variability in the total transport ( Tr ± T&, Fig. 2-3, 
shaded areas) may be ascribed to the depth-independent contribution (T&, not shown) 
while a small fraction (10 to 30% across 25°N and 36°N) is reflected in the relative trans¬ 
port Tr. On the other hand, the relative transport variability dominates in the surface 
layers. The variability in the thermal wind transport, Tr, changes with latitude but gener¬ 
ally it does not exceed ±0.5 Sv below 2000 m while it increases from ±1 Sv to ±2 Sv above 
2000m (all numbers are the one standard deviation variability). Variability in the deep 
layers is higher at 48°N in the 2000-4000 m range, possibly because the deep western bound¬ 
ary current variations are stronger there. The variability in the total relative transport 
(numbers below each plot) is relatively large: ±6.2 Sv rms(25°N), ±5 Svrms (36°N) and 
±2.5 Svrms (48°N). On short time scales, such deviations are in balance with a barotropic 
(or Ekman) return flow. However, those variations must be taken into account for com¬ 
bining consistently different sections, as discussed later (page 64). 
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T b +T R A25N T r A25N 




T+T A36N 

b R 



total: 1.2 ± 6.5 Sv 


T b +T R A48N 



total: 8.4 ± 3.7 Sv 


T r A36N 



total: 25 ± 4.9 Sv 


T r A48N 



total: 41 ± 2.5 Sv 


Figure 2-3: Individual layer mass transports across 25°N (upper), 36°N (middle) and 48°N 
(lower) in the GCM. The plots on the left give the total transport while the plots on the 
right give the thermal wind contribution. The eight year average is indicated by the lines 
while the variability is given by the shaded area (one standard deviation). 
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In principle, in the real ocean, the net transport across each section (numbers below 
the graphs on the left) is equal and opposite to the average Ekman transport; but because 
of the ageostrophy in the GCM, this balance is not achieved by the time-mean geostrophic 
flux, producing an error of order 10 Sv at 36°N. 

For comparison, similar calculations were done over a single year (not shown). The 
variability in Tr is smaller by a factor of 2 in the deep layers while it is similar in the 
surface layers and for the total flux. 

As mentioned above, the variability may be underestimated by a factor two or four, 
according to the GCM evaluation of Stammer et al. (1996). Roemmich and Wunsch (1985) 
compared the geostrophic flux through repeat sections at 24° N and 36° N taken about 25 
years apart (1957-1981). They found deviations of order 1 Sv in individual isopycnal 
layers (the same layers as I used) and up to 2Sv near the surface. Those errors were 
found to be systematic over several layers, leading to relatively large deviations in water 
mass transports: up to 7Sv. The authors pointed out that part of the deviation may 
come from actual oceanic variability while the other part may be an enhancement of the 
variability by the simple model used, i.e., a uniform correction to the reference velocity. 
Direct comparison with the GCM results is somewhat difficult because different reference 
surfaces were used and because the Roemmich and Wunsch transports are already adjusted 
for mass conservation over the two sections so that their reported variability includes a 
barotropic adjustment. Still, the magnitude of the variations is consistent with those 
derived from the GCM. 

At 24°N, Lavin (1993) did a similar comparison including the more recent (1992) 
section at 24°N. The top-to-bottom geostrophic transports on the three occupations from 
Bahamas to Africa referenced to 1100 m were -41.0 Sv (1957), -31.7 Sv (1981) and -36.9 Sv 
(1992). Assuming a constant transport in the Florida Strait (e.g., Schott et al., 1988b) the 
geostrophic variations range between +9.3 Sv (1957—->1981) and —5.2Sv (1981—>4992). 
Both deviations are consistent with our ±6Sv GCM variability in the total flux at 24°N. 

At 24°N, 36°N, and 48°N, Koltermann et al. (1999) compared the transports from 
repeat sections occupied since the IGY through recent intensive WOCE surveys. Based 
on the assumption of Sverdrup balance, and after adjustment to the initial velocity field, 
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the meridional overturning in the interval 1957 to 1997 varied between 8 and 21 Sv at 36°N 
and between 10 and 20 Sv at 48° N. This variability is larger than the one inferred from the 
GCM. However, not only does it include part of the barotropic variability after correction, 
but it also relies upon the Sverdrup relation, which is inappropriate in the sub-polar North 
Atlantic Ocean (Luyten et al., 1985). As a result, it is not possible to distinguish the part 
of this variability which is due to changes in the density field alone, which I explored 
in the GCM. Also, at 48°N, measurement errors are enhanced because the deep western 
boundary current flows against the continental shelf (see below, Section 2.3.3). 

At 48°N in Koltermann et al. (1999) and the GCM, a shallow reference surface was 
used. In the presence of oceanic eddies, such a choice creates large velocities at the bottom, 
reflecting the large eddy velocities at the thermocline depth. Above rough topography, 
these eddy-related bottom transports do not cancel out, creating an artificially large vari¬ 
ability in the transport. Hence, starting with a deep reference surface is more likely to 
minimize this error, as long as enough range of adjustment is allowed in the reference 
velocities. 

On the whole, the GCM variability at 24°N, 36°N and 48°N is in agreement, within a 
factor of two or smaller, with the inferences from repeat sections. 

Eddy aliasing 

Another problem that arises while sampling a constantly varying ocean is aliasing. Two 
sorts of alias take place: l)eddy “Doppler” aliasing and 2) spatial aliasing. The eddy 
“Doppler” appears because it takes one month for a ship to cross the Atlantic, so that 
the field “seen” by the ship is different from any instantaneous field. An extreme example 
would be a ship moving at the same speed as an eddy: the T, S measurements would show 
flat isopycnals over the section. In a more realistic situation, eddies moving with the ship 
exhibit smoother gradients, while eddies moving against the ship exhibit steeper gradients, 
with possible effects on property transports. 

To explore the effect of this possible alias on the mass flux, cruises were simulated along 
36°N in the GCM so that the outputs were sampled in one month across the Atlantic. The 
mass flux was then compared with the annual average flux obtained from synoptic model- 
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data. One cruise was simulated for each month of 1993. The layer transports did not 
show departures from the 1993 average that would exceed the natural variability reported 
in the preceding section. Still, the effect of aliasing is visible in the top-to-bottom thermal 
wind transport: Tr varied by ±4.8 Sv rms between the 12 cruises while the variability 
from instantaneous sections at 36°N from 3-day snapshots and for one year is ±3.5 Svrms. 
The difference, 40% of the natural variability, does not exceed the ±5 Sv rms interannual 
variability there and for this “order of magnitude” error budget, the ship aliasing effect will 
be neglected with respect to the natural variability. Because at mid-latitudes a research 
ship crosses the Atlantic about 100 times faster than the first baroclinic Rossby wave 
(assuming it takes 10years for a baroclinic Rossby wave to cross the ocean), aliasing is 
unlikely to introduce a large error. 

The second possible error, spatial aliasing, is caused by finite station spacing, which 
is generally close to the first Rossby radius. Because much of the eddy energy occurs 
on scales comparable to the first Rossby radius, the data may alias the scales that are 
below the station spacing into large scales. Unfortunately, the GCM output is inadequate 
to quantify this error because subsampling the GCM fields would incur errors from the 
representation of the bottom topography. Roemmich and Wunsch (1985) compared the 
transports between two coarse-resolution and two fine-resolution sections taken at different 
times in the North Atlantic. They concluded that temporal variability dominated the 
changes in geostrophic transports. Still, it is not certain that modern sections have enough 
resolution. A possible test would be to sub-sample a high resolution XBT line to study 
the effect on property transports. To my knowledge, no such study has been done and I 
have to assume that the spatial aliasing error is not large compared to the other errors. 

2.2.4 Silica Variability 

Silica is a useful constraint in inverse box models because to a good approximation it 
is believed to be conserved in the ocean (e.g. Section 4.3.1; Robbins and Bryden, 1993; 
Robbins and Toole, 1997; Rintoul and Wunsch, 1991; Wunsch et al., 1983). Like most 
nutrients, the dissolved silica concentration is depleted in shallow waters. Consequently, 
the silica flux does not behave like the mass or energy flux. As for the mass conservation 
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budget, the silica variability introduces a noise n™° del in the silica conservation equations. 
In this section, I provide an estimate for the silica flux variability for each layer in the 
Atlantic at 36° N. The variability in the meridional flux of silica is impossible to estimate 
from existing observations and the GCM does not contain any nutrient information. To 
remedy this situation, an artificial silica variable was created in the GCM based on a 
simple relation between potential temperature, longitude and silica. I first constructed 
this relation, based on a hydrographic section at 36°N, then, I implemented it in the 
GCM output fields to obtain the variability in the advective silica flux. 

Figure 2-4 shows the silica concentration versus temperature in the Atlantic at 36°N. 
An analytical function was sought in the form of Chebyshev polynomials tapered with an 
exponential to remove the dominant vertical structure: 

6 

Si(t,l) = ^2 a i x Ti(6) x e - 1Ox (0-i) _|_ a7 x i (2.10) 

2 = 0 

where S; is the silica concentration, l, 0 are the longitude and the temperature scaled to the 
interval [-1,1]. Ti(t) is the ith Chebyshev polynomial (type 1) and the coefficients a* are 
obtained by regression. The regression was done independently for each side of the mid- 
Atlantic ridge as the silica pattern varies significantly. The regression curves are overlain 
on Figure 2-4. Polynomials of degree 6 were sufficient to represent approximately the 
behavior of silica with temperature. Residuals are shown on Figure 2-5. The rms residuals 
do not exceed 2/rmol/kg, which is much lower than the average silica concentration. 

Using the analytical representation, the silica field was computed from the GCM out¬ 
put every third day for eight years. Figure 2-7 compares the observed silica from the 
hydrographic section (upper panel) with a snapshot of the “GCM silica” (lower panel). 
The GCM silica is smoother than the observations because of the smoother temperature 
field. The step above the mid-Atlantic ridge arises from the change in the regression 
function. This unrealistic feature did not introduce significant spurious variability. 

The geostrophic flux of silica was obtained, as for the volume flux, every third day 
through each of the isopycnal layers defined in Section 2.2.3. Figure 2-6 shows the total 
mean flux and variability (Tr + 7),, left plot) and the thermal wind contribution (Tr, 
right plot). Because of the increase in the silica concentration with depth, the variability 
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Silica fit residuals, a36n 



Figure 2-5: Silica residual after regression, at each depth. For comparison, the mean silica 
concentration is overlain. Silica is in /xmol kg -1 while depth is in decibar. (The accuracy 
of very best routine measurements of silica is 0.5-1.0/imolkg _1 ). 
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Figure 2-4: Silica versus potential temperature at 36°N in the Atlantic. Data east of 
the mid-Atlantic ridge are represented by pui-ple open circles while on the west they are 
represented by the blue dots. The two curves (respectively green and blue) correspond to 
regressions obtained on either side of the mid-Atlantic ridge. 

in the total flux now increases with depth, in contrast with the mass flux variability. 
However, the variability due to the thermal wind, ^(«™„ pdel ) 2 )> is relatively constant 
with ±10kmolSis _1 in each layer. The variability in the net flux from thermal wind 
is about ±50kmoISis -1 . While I believe that my results give a lower bound on the actual 
oceanic variability, the mean GCM silica fluxes (Fig. 2-6, lines) are probably unrealistic 
because 1) of the ageostrophic issue discussed above and 2) the re-created silica field is 
not necessarily compatible with the average thermohaline circulation. Silica was chosen 
among other nutrients because of its relatively simple biogeochemistry, but the results are 
adaptable to other nutrients that have approximately the same vertical structure. 
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Figure 2-7: (Upper) Silica observed from the 36°N hydrographic section. The snapshot of 
the GCM reconstructed silica (Lower) shows a similar, but smoother field. 
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T b+ T R A36N T r A36N 




Figure 2-6: Silica variability at 36°N. The shaded area corresponds to the variability in 
each layer. The line gives the average transports through each layer in kmolSis -1 (see 
text). The left plot gives the total transport (Tr + T b ) while the right plot gives the 
thermal wind contribution (Tr). 

2.3 Measurement noise 

I now turn on to the measurement noise, n me asurement q^e measurement noise has two 
types of components: the first one comes from the temperature and salinity measurements 
while the second one is due to the approximation in the transports calculations. I first 
explore the temperature and salinity noise, which is mainly due to internal waves; then 
I discuss two main approximations that are used: geostrophy and those in the bottom 
triangles. The error of the average Ekman transport is discussed in the end. 

2.3.1 Internal waves 

Measurement noise in the data from a hydrographic section has at least three different 
sources: instrumental, navigation and internal waves. Instrumental errors are probably the 
smallest. CTDs have reached an absolute accuracy of approximately 10 -3o C for tempera¬ 
ture and better than 4 x 10~ 3 for salinity when carefully calibrated (Saunders, 1986). (The 
temperature requirements of the WOCE standards is of 0.005°C accuracy and 0.002°C pre¬ 
cision (Joyce and Corry, 1994).) The corresponding noise in the density is 3 x 10~ 4 kg/m 3 
from temperature noise (l(T 3o C) and 3 x 10 _3 kg/m 3 from salinity noise (4 x 10 -3 ). As 
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we will see below, these deviations are smaller than the noise caused by internal waves. 
Improperly calibrated CTD’s may bias whole data sets by ±0.01°C and ±0.01 psu (See 
Marotzke and Willebrand, 1996, p.64 for a review) and cause density shifts over a whole 
section, with no effect on geostrophic transport calculations. Navigation is also a minor 
source of error, as it is now very accurate. Moreover, except near high currents (that in¬ 
duce a drift of the rosette during the station) or steep topography (in the bottom triangle), 
the errors in station positions to cancel out in transport computations. 

The main noise in temperature and salinity measurements is due to internal waves. The 
transport due to the geostrophic velocity referenced to the bottom between two locations 
along the section, S w and S e , may be written 

iz L dC ii d £ (i ' x) ' 

where h is the bottom depth, p the density, g the gravity and / the Coriolis parameter. If 
if the bottom is flat (h independent of x), T only depends upon the density measurements 
at the two end points: 


rS e r 0 rS e rO 

/ dx dzv(x,z) = — dx 
J S w J—h J S w J—h 


T — 


Pof 


f dz f dC[p(C,^)-p(C,5e)]. 
J-h J-h 


( 2 . 11 ) 


To estimate the noise in T, denote the observed transport as T ohs and the true transport 
as T true . Write 


jiobs _ true ^ 

where nj- is the observational noise. From (2.11), may be written in terms of the noise 
in the density, n p : 

n T = —,f dz f d( [n p (C, S w ) - n p (C, S e )], (2.12) 

POJ J-h J-h 

so that the variance of tit is 

( n r) = jj dzdz> J J ^ d C d C ([ n P (C, S e ) - n p ((, S^)] [ n p (C, S e ) - n p (C', 5^)]) , 

where the bracket () is the expectation. 
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Internal waves are one of the best explored fields of physical oceanography and in 
principle one can compute the full covariance of n p using the three dimensional spectra 
that have been estimated from observations. For my “order of magnitude” error study 
I used the Garrett-Munk spectrum (Munk, 1981). As observations show no evidence of 
horizontal coherence over distances larger than 10km (Garrett and Munk, 1972), one can 
assume that n p (S w ) and n p (S e ) are uncorrelated (S w and S e can be taken either between 
two stations or at the two ends of a section). We then have 

{nl)= (ii) x 2 x M) - n > &)■ (2 - 13) 

To compute the integrals, the water column is divided into segments Az 2 - of 500 m or 
more so that the internal wave noise is approximately uncorrelated between each segment. 
Eq. (2.12) is evaluated as 


M k 

n T = Azk AZi X t U P ( Zl ' Se ) ~ n P ( Z i’ Sw )] 

P0J k =1 i=l 


x a r K [n p (Se) - n p (5w)] , 


where the column vector n p contains n p (zi) and the vector-matrix product a T K corre¬ 
sponds to the two sums. A z k is the length of segment k. Then, (2.13) may be evaluated 

<nb=2*(^)V K <",, n J> KT a, (2.14) 

where the covariance matrix (n p nj) is diagonal. The elements of (n p nj), namely (n 2 p ) fa), 
are chosen for each depth following Marotzke and Willebrand (1996b) with slight modifica¬ 
tions. First, the variance of the isopycnal displacements is estimated using the Munk (1981) 
formula, 


where JVq = 5.2 x 10 3 s 1 and N is the local Brunt-Vaisala frequency. The conversion 
into density variance gives 

(r»;> = 53^ x NqN (z) 3 . (2.15) 
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N varies from 0.8 x 10~ 3 s _1 in the deep layers to 5x 10 -3 s -1 in the thermocline at mid¬ 
latitudes (e.g., Gill, 1982, p.52). The corresponding noise n p goes from 0.0013 kg/m 3 rms to 
0.02 kg/m 3 rms. Because it is tiny in the deep ocean, n p was given a lower bound equivalent 
to the CTD noise (0.003 kg m -3 rms). At a latitude of 30° and for h = 4000m, the rms 
noise in the transport (Eq.2.14) was integrated using 500 m segments. The resulting noise 
in the total transport is 3Svrms. This result is found insensitive to a change in the 
segment of integration, and a depth of h = 5000 m produces the same order of magnitude. 
The noise at different latitudes is obtained by scaling the ±3 Sv rms with the appropriate 
ratio of Coriolis parameter. 

Fieux et al. (1996a) observed large internal fluctuations in the transports between re¬ 
peat stations taken a few hours apart in the Bali region. The authors suspected that 
these fluctuations were due to internal wave activity. They found variations of order 9 Sv, 
comparable to the noise of ±10 Sv predicted by Eq.(2.14) at this latitude (10°). 

2.3.2 Geostrophy 

In the ocean interior, and except for the Ekman layer, geostrophy is recognized as a good 
approximation. Even in the western boundary current, equation scaling shows that the 
along-stream velocity is in geostrophic equilibrium (see Pedlosky, 1987, Chapter 8). The 
very first test of geostrophy is probably the one of Wiist (1924) who compared current- 
meter data with the geostrophic velocities in the Florida Current and found a qualitative 
agreement (see also Wunsch, 1996, p.47). Recently, Johns et al. (1989) directly tested 
geostrophy in the Gulf-Stream at 36°N using Pegasus profiler currents. They obtained an 
agreement within 3 Sv between the absolute and geostrophic transports above 2000 dbars, 
only a 3% error in the total transport of 93 Sv. Using a satellite image of the Gulf-Stream 
position at the time of the experiment, the authors estimated the correction from the 
cyclostrophic term to be about 2 Sv, which further reduced the discrepancy to 1 Sv. This 
result is most encouraging, although the measurement errors are probably much larger. 

In addition to the cyclostrophic term, ageostrophic motions include the energetic in¬ 
ertial oscillations. The extent to which inertial motion is able to carry energy or other 
tracers when integrated over a section is not trivial in presence of the topography and the 
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planetary vorticity gradient. Ripa(1997) gives the analytical description of inertial oscil¬ 
lations in the presence of the planetary vorticity gradient. For a surface inertial velocity of 
order V = lOcms -1 , the inertial circles have a radius R = V/f = 1000 m at mid-latitudes 
and a westward “Stokes” drift appears (Ripa, 1997) 

{i) = ~\p r2 ~ - 5 x 10 ~ U x 10 6 = O (5 x 10~ 4 cms' 1 ), (2.16) 

where f3 is the planetary vorticity gradient. The magnitude is much smaller than the 
geostrophic horizontal advection and purely zonal with no implications for the meridional 
tracer transports. 

A second effect of the planetary vorticity gradient is to shift the inertial frequency 
slightly above / (e.g., Munk, 1980). The pressure field is affected ageostrophically, cor¬ 
rupting the hydrographic measurements. Such a peak above the inertial frequency has 
been observed in some of the internal wave potential energy spectra (Fu, 1981). However, 
the author noted that the peak may be spurious and a result of inertial motion of the 
moorings themselves. I do not know how much this ageostrophic effect aliases transport 
calculations, and it is assumed that the implied error does not exceed the internal wave 
error of Section 2.3.1. 

In conclusion, the effect of inertial motion is much smaller than other uncertainties in 
the budget. Still, a ±2Svrms uncertainty may be attributed to the cyclostrophic term 
across the Gulf-Stream at 36°N. 

2.3.3 Bottom triangle 

In the presence of sloping topography one has to decide how to compute geostrophic 
velocities and property transports between two stations near the bottom. The situation 
for an actual hydrographic line is sketched on Figure 2-8a. For computing geostrophic 
transports, one usually extrapolates the temperature, salinity and other properties beneath 
the deepest common level of the two stations (noted hereafter DCL). For instance, at the 
two first stations of Fig. 2-8a, the temperature was measured on the black dots in the 
triangle below the DCL. It must be extrapolated along the bottom slope or, equivalently, 
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Figure 2-8: Bottom triangle problem. The topography is shaded, the first line corresponds 
to the shallow station, the second line to the deep station. The black dots are temperature 
measurements, a) in the real ocean; b) in the model; c) in the subsampled model. DCL 
stands for the deepest common level. 

on the imaginary points in white that are below the DCL at the shallow station. 

The temperatures at the white points are denoted T s h a i while the temperatures at the 
deep station and below the DCL are noted Tdeep■ The salinities Sdeep and S s hai that are 
extrapolated with the same method as temperatures and the geostrophic velocities are 
computed using the extrapolated T s h a i and the data Tdeep and S s h a i, Sdeep- The transport 
is computed by integration of the velocity in the triangle. 

Uncertainties arise from both the extrapolation of T s h a i and the estimation of the area 
in the bottom triangle. Here, I estimate the uncertainties related to the extrapolation only 
by comparing six methods: 

• Method 1: null velocity below the DCL. The transport in the triangle is ignored. 

• Method 2: no shear. T s h a i is set equal to Tdeep at each depth so that the geostrophic 
velocity is constant below the DCL. 

• Method 3: constant slope. The slope of the isotherm just above DCL is calculated. 
T s hai is computed using the same slope (vertical variations to horizontal variations). 
The magnitude of the slope is limited to an upper bound of 1 (regardless of the 
vertical to horizontal distance ratio) to exclude excessively large shears. 

• Method 4: plane fit (suggested by Wunsch, 1996). A linear function is determined 
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by regression between the temperature data, the depth and distance on the triangle 
(black dots): 


T(x,z) = ax+ bz + c, (2.17) 

where x is the distance to the deep station. The relation gives T s ^ a ;. 

• Method 5: polynomial fit. Same as the plane fit, but a higher order polynomial is 
used in the 2 direction: 


T(x,z) = ax + bz + cz 2 + d. (2.18) 

The polynomial is useful in regions of strong T curvature (e.g., near the thermocline). 

• Method 6: horizontal extrapolation. Temperature is extrapolated linearly from 
the next deeper station. 

Other extrapolation methods such as mapping have been used in the literature (e.g., 
Roemmich and Wunsch, 1985), all bearing their own advantages and drawbacks. To de¬ 
termine the possible error caused by the extrapolation, the relative geostrophic transports 
were computed using each one of the six methods listed above and across three sections 
in the Atlantic: 24°N (Morocco to east of the Bahamas Bank); Florida Strait (here¬ 
after Fist); 36°N and 48°N. The geostrophic velocities were referenced to deep isopycnals: 
^4000 = 45.81 at 24°N and 36°N and <72000 = 36.87 at 48°N (from Macdonald, 1995) or 
to the bottom if shallower than the target isopycnal. For those “shallow stations”, the 
transport calculation was done with the velocity referenced 1) at the depth of the deepest 
station and, 2) at the DCL of the two stations. The total relative transports are displayed 
on Figure 2-9 for each method and in the two cases. 

In the first case (Fig. 2-9, upper panel) the relative transports vary from 6 to 31 Sv at 
24°N depending upon the method. While it is relatively stable in the Florida Strait (±2 Sv) 
it exhibits large variations at 36°N (+45 Sv) and at 48°N (+5 Sv). Figure2-10, lower panel, 
shows the transports cumulated from the eastern boundary at 36°N for each method. The 
disci epancies appear almost exclusively in the regions where the velocity is referenced to 
the bottom. Because the velocity there is referenced to the deepest station depth, any 
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Figure 2-9: Net relative transport across each section and for different bottom wedge 
methods. The calculation is much more sensitive to the method if the reference in shallow 
waters is taken at the deepest stations (upper panel) than at the last common depth (lower 
panel). 
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error in the bottom triangle is propagated through the whole water column. An extreme 
example occurs along the eastern slope (Fig.2-10): Between 1000 and 2000dbars, a strong 
density gradient due to the Mediterranean outflow coincides with a steep bottom slope. 
As a result, the shear between the DCL and the bottom depth varies by 1 ms -1 depending 
upon the method, introducing a 30 Sv anomaly with the “constant slope” method. While 
this extreme anomaly in the velocity field is readily identified and corrected, smaller errors 
are not obvious while their effect may be important. 

A straightforward way to minimize the error induced by the bottom wedge extrapola¬ 
tion is to reference the “shallow pairs” to the DCL instead of the deepest station. Then, 
the uncertainty in the bottom triangles has no effect on the water column above DCL. 
The transports in this case are displayed on Fig.2-9, lower panel. (Even if the reference 
isopycnal was between the DCL and the deepest station, the reference was taken at the 
DCL.) Within 1 Sv, and except for the polynomial fit at 24°N, the relative transports are 
independent of the method at 24°N, in the Florida Strait, and at 36°N. At 48°N, the 
scatter is more important (4 to 9 Sv) because of the presence the deep western boundary 
(DWBC) across several bottom triangles. 

The choice of an extrapolation method is somewhat arbitrary. In principle, one could 
use the GCM output, where the answer is known, to select the best method. Unfortunately, 
the GCM configuration is inadequate: Because of the non-slip condition at horizontal 
boundaries, the actual velocity in the triangle-which is centered on a velocity point-is 
null (Fig.2-8b). Alternatively, the GCM may be sub-sampled at every other grid points 
(Fig. 2-8c), but with the present model resolution, the resulting station spacing would be 
much larger than it is in actual hydrographic data. 

Experiments on many hydrographic sections showed that methods 2-5 (constant veloc¬ 
ity, constant slope, plane fit and polynomial fit) propagate the temperature and velocity 
structure downwards in the water column. If the depth difference between the two sta¬ 
tions is large, unrealistic temperature gradients and velocity shears can be introduced at 
depth in the bottom triangle (e.g., if the shallow station is in the thermocline and the deep 
station is at 3000dbars). Ignoring the transports in the triangles may be safer but not 
realistic, in particular when it implies a large discontinuous shear at the DCL. The “con- 
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Figure 2-10: Relative transport at 36°N integrated from east. The reference surface is 
indicated by the dash-dot line on the upper panel. The transports in the bottom triangles 
are computed with different methods. The “shallow” pairs are referenced to the deepest 
station, causing the large discrepancies. 















stant velocity” method is relatively safe, but it implies an abrupt flattening of isopycnals, 
which is not necessarily realistic either. On the other hand, the horizontal extrapolation 
has the advantage of reproducing the velocity shear from the next deeper station pair, 
which is by definition realistic. This method may be unstable in the rare case of two 
closely spaced stations followed by a distant station along a rising bottom slope. Then, 
the noise in the close stations is amplified through the extrapolation to the distant station. 
For those particular cases, the extrapolation may be done manually, using for instance the 
“constant velocity” method. For stations where no horizontal extrapolation is possible 
(i.e., the station is isolated in a trench), the polynomial fit turns out to produce the most 
realistic-looking temperatures as long as the extrapolation depth interval is small, because 
the technique permits curvature in the extrapolated temperatures. 

In conclusion, my estimation of the error due to the bottom wedge problem on a given 
transoceanic section is of order 1 Sv as long as the reference level of shallow stations is 
taken at the deepest common level of the two stations. In the presence of a boundary 
current along sloping topography, this error is probably, 2 to 4 Sv. Because this estimate 
derives from a comparison between various extrapolation methods and a null bottom 
triangle transport, it is likely to be an overestimate of the actual error. Further tests may 
be carried out to estimate the bottom triangle error (e.g., subsampling a high resolution 
section or numerical model; comparing calculations with LADCP measurements). 

2.3.4 Mean Ekman transport uncertainty 

As pointed in Section 2.2.2 it is appropriate to use the mean Ekman transport rather 
than the instantaneous value in property budgets. In Section 2.2.2, it was shown that 
the classical Ekman formula (2.5, page 28) was adequate for estimating the net Ekman 
transport in the GCM. However, due to the large variability in the Ekman layer, Eq.(2.5) 
has no definite support from observation. But this relation between the wind stress and its 
implied mass transport has never been proven wrong. The few existing measurements are 
encouraging (e.g., Chereskin, 1995; Price et al., 1988; Schludlich and Price, 1998; Weller 
and Pueddemann, 1996). A convincing test of the net meridional Ekman transport comes 
from Chereskin and Roemmich (1991) who compared, along a hydrographic section (11°N 
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in the Atlantic), the total Ekman transport (Eq.2.5) with a direct calculation of the surface 
ageostrophic component based on hydrographic and ADCP measurements. The different 
estimates agreed within their estimated error bars, 12 ± 5 Sv from the ADCP against 
13.5 ± 0.3 Sv from (2.5). Wijffels et al. (1994) arrived to a similar conclusion at 10°N in 
the Pacific Ocean. 

The accuracy of the wind stress itself is uncertain. It depends on l)the accuracy of 
the estimated wind itself and, 2) the validity of the momentum transfer parameteriza¬ 
tion at the sea surface. Mestas-Nufiez et al. (1994) compared the monthly wind stresses 
from the climatology of Trenberth et al. (1990) with the Seasat-A Satellite Scaterrometer 
(SASS) stresses. They found that the Trenberth stresses were overestimated by 46% in 
most oceanic regions, except at the Equator and at high southern latitudes where the 
errors became larger. Most of the difference was explained because Trenberth et al. (1990) 
used the 1000-mbars winds to compute the stress instead of the 10-m winds. Huang and 
Shukla(1996) arrived at a similar conclusion in comparing ECMWF to ship wind obser¬ 
vations in the tropical Atlantic. Note that the ECMWF wind stresses that were used in 
Section 2.2.2 are obtained from a Monin-Obukhov-type boundary layer parametrization 
(ECMWF forecast model physical parametrization, 1991) which is probably more accurate 
than the drag-coefficient-based analysis of Trenberth et al. (1990) . 

In the present hydrographic inverse model, the mean Ekman transports were computed 
from the NCEP re-analysis (Kalnay et al., 1996). It is difficult to provide a reliable number 
for the a priori uncertainty on the average Ekman transport because of the uncertainty 
in both the wind stress and the calculation itself (Eq.2.5). Given the discrepancies cited 
above (46%), an uncertainty of 50% was attributed to the mean Ekman transport. For 
sections of low Ekman transport, the uncertainty was given a lower bound of 1 Sv. 

Because of possible correlations between velocity and temperature in the Ekman layer, 
the error in the energy flux calculation in the Ekman layer can be much larger than that 
in the mass flux itself. Variability in the net energy flux through a trans-oceanic section 
is dominated by the vertical temperature gradients while the gyre component is a minor 
contribution (Jayne, 1999). A sensible test of the average energy transport would be to 
compare vT and vT (where T is temperature and the overbar is the horizontal average) 
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from ADCP data as in the study by Chereskin and Roemmich (1991). But this is beyond 
the present work, and I assume that the uncertainty in vT includes any possible eddy 
transports v'T' . Nutrient fluxes are much less affected by horizontal fluctuations in the 
Ekman layer because nutrients tend to be depleted near the surface 
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2.4 Reference level velocity variability 


In principle, the a priori reference velocity covariance, (bb T ) = (b(x)b(x')), could be fully 
determined from the GCM. However, the calculation is complicated because (bb r ) is 
likely to be non-isotropic and geographically dependent. Moreover, there is no guarantee 
that the GCM produces the true oceanic covariance structure because it does not entirely 
resolve oceanic eddies. Consequently, I make the assumption of no correlation between 
elements, that is, a diagonal covariance matrix. To provide at least a lower bound for the 
diagonal elements of (bb T ), the GCM output is examined. 

At the three latitudes 24°N, 36°N and 48°N, the absolute geostrophic velocity was sam¬ 
pled at the reference isopycnal every third day in the GCM. (As above, I took the reference 
isopycnals <74 = 45.81 at 24°N and 36°N and 02 = 36.87 at 48°N from Macdonald, 1995.) 

Figure 2-11, upper panel shows the average reference isopycnal position at 36° N ((74 = 
45.81, thick line) over 8 years. The envelope defining one standard deviation of the isopy¬ 
cnal displacements is indicated by the thin line. Deviations are lower than ±100 m, sug¬ 
gesting that the variations in the reference velocity are thus mainly due to variations in 
the velocity field rather than related to reference isopycnal displacement over time. 

The time-mean reference velocities (Figure 2-11, lower panel, thick line) are generally 
smaller than 0.2 cm s~ ! in the oceanic interior while the variability, indicated by the shaded 
area, is approximately ±0.5 cm s _ 1 rms. The mean reference velocities increase to 20 cm s -1 
near the western boundary and to 5cms _1 in shallow waters to the east. The situation 
is similar at 24°N (not shown) while at 48°N the amplitude of both the mean and the 
variability are twice as large. At this latitude, Macdonald (1995) chose a relatively shallow 
isopycnal of reference (a? = 36.87,1200-2000 m). Nevertheless, the reference velocities on 
a deeper isopycnal ( 0-4 = 45.81, 3000-3600 m) have the same amplitude and variability. 
Mean velocities of up to 6 cms _1 are found in the deep western boundary current at 48° N 
near the continental slope. 

It is not clear whether the mean reference velocity size, its variability, or its rms is the 
right a priori covariance for b. On one hand, the variability contains the fast responses to 
Ekman transport fluctuations. Those are filtered out in the mass balance by specifying the 
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Figure 2-11: Upper panel: reference isopycnal (a 4 = 45.81) for section 36°N. The isopycnal 
oscillates by less than 100 m around its mean position (thick line). One standard deviation 
is indicated by the thin lines. Lower panel: mean reference velocities (thick line) and 
variance (lstd. dev, shaded area). 
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time-mean Ekman transport (see the introduction, Section 2.1). On the other hand, the 
model variability in the reference velocities also contains the eddy part, both baroclinic and 
barotropic: We have seen that the relative transport Tr ( t ) is stable when integrated on 
large scales. However, Tr contains baroclinic eddies that have a signature in the reference 
velocities (unless all eddy velocities somehow vanish at the reference depth). As a result, 
the estimated mean transport Tj, may contain some eddy signals to be consistent with the 
instantaneous realization Tr (t). 

In general though, the reference velocity and its variability have similar magnitude. 
Because the GCM reproduces about half the actual oceanic variability (see Section 2.2), 
an a priori value for the reference velocity variance in the interior of 1 cm 2 s~ 2 , as used 
by Macdonald (1995) was adopted. In shallow waters and near boundary currents, this a 
priori value is increased according to current meter statistics when available. Except in 
boundary currents where an accurate reference surface could be determined from current 
meter measurements, the relative velocity variations within a range of about 500dbars 
was taken as the a priori value to reflect the uncertainty in the reference surface choice, 
with a minimum of 1 cm 2 s~ 2 . 


61 




2.5 Summary 


This evaluation of the a priori uncertainty in mass conservation equations of linear inverse 
models shows that: 

1. The variability in the absolute geostrophic velocity field is dominated by fast, depth- 
independent responses to Ekman transport fluctuations. Consequently, linear inverse 
models can consistently combine different sections as long as the mean Ekman trans¬ 
port is used for all sections. 

2. The baroclinic shear is relatively stable on large scales (at least in the GCM that was 
examined) so that the relative baroclinic transport measured at a particular time is 
representative of the time-mean. Still, the error implied by baroclinic variations is 
not negligible, with a total of ±5 Sv at mid-latitudes for the top-to-bottom transport 
and ±0.1 ±2 Sv in the selected individual layers. The finite time required by a ship 
to cross the ocean does not seem to introduce any significant bias. 

3. The measurement noise is dominated by internal wave vertical heave of isopycnals, 
implying a ±3 Sv error in the transport estimates at mid-latitudes. The geostrophic 
assumption and the bottom triangle issues induce smaller error except in very specific 
situations such as a boundary current with large curvature or a deep boundary 
current against a slope. 

Table 2.2 summarizes the a priori uncertainties on the top-to-bottom flow and on the 
flow in isopycnal layers above and below 2000 dbar. (The layers were defined in Sec¬ 
tion 2.2.3.) The uncertainty in the top-to-bottom flow is dominated by the baroclinic 
variability of ±5Sv. Assuming that the different uncertainties are independent, the total 
rms uncertainty is 6 Sv. In a conservation equation involving two sections, the noise level, 
assuming uncorrelated errors between the sections, would be 8.5 Sv rms, suggesting a 4-fold 
increase of the 2 Sv rms uncertainty adopted by Rintoul (1988) and Macdonald (1995). 

Using such a large error for the net flux through a section is somewhat counter-intuitive 
as it is known, independently from inverse models and hydrographic data, that mass is 
conserved within the uncertainty on the freshwater flux divergence, ±1 Sv or less. But 


62 



because the model is linear, with no adjustment to density, the model error due to the 
variability in the density field must be included to obtain reliable estimates and uncer¬ 
tainties for the time-mean property fluxes. Figure 2-12 illustrates the effect of imposing 
stringent mass conservation between two sections in the presence of baroclinic noise in a 
two-dimensional case. The true average is an overturning of 5Sv (upper panel). Because 
of baroclinic variability in the upper layer, a shear of only 5 Sv is observed in the first 
section (red dashed line). The second section is supposed to be perfectly measured. In 
the lower layer, the baroclinic variability is much weaker. If one enforces top-to-bottom 
mass conservation, the derived circulation (middle panel of Figure 2-12) has half the over¬ 
turning strength and spurious vertical transfers are implied by mass conservation between 
the two sections. The weak overturning is probably the true circulation at the time of 
the hydrographic section; however, it is a misrepresentation of the time-mean, leading to 
erroneous vertical transfers. But according to the present GCM simulation, the baroclinic 
shear is more stable in the deep layers. By ensuring a better conservation in the deep layer 
at the expense of the top-to-bottom conservation, one gets a circulation which is represen¬ 
tative of the time average but it contains mass residuals in the upper layer (lower panel of 
Fig. 2-12). This solution is correct as long as the mass residuals are indistinguishable from 
zero within uncertainties. A numerical illustration in a two-dimensional ocean is given in 
Appendix B. 

Along with mass, the model error for the silica flux was estimated because silica is a 
useful tracer to constrain the circulation (e.g., Rintoul and Wunsch, 1991). The uncer¬ 
tainty due to the baroclinic variability in the GCM output was estimated at ±50 kmol Si s _1 
top-to-bottom and at ±10 kmol Sis -1 in individual layers. 

The a priori reference velocity magnitude was evaluated from the GCM as lying be¬ 
tween ±0.5 and ±lcms _1 . In shallow waters and near boundary currents, the values are 
much larger. Practically, an a priori value of ±lcms _1 will be used as a lower bound 
and modified according to 1) current-meter measurements where available and 2) the local 
shear to account for the uncertainty due to the reference surface choice. In principle it 
is possible to use the spatial structure of the reference level velocity as prior information. 
For example, one can minimize the solution smoothness instead of its variance, which is 
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Figure 2-12: Enhancing conservation in the deep layers at the expense of the top-to-bottom 
mass budget: two-dimensional example. Panel 1) shows the true average circulation, an 
overturning of 5 Sv at both sections. The red dotted line corresponds to the observation of 
the section to the left. Panel 2) and 3) are two estimates of the circulation that conserve 
mass (1) top-to-bottom and (2) in the deep layer. 
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Top-to bottom 

Layer above 2000 m 

Layer below 2000 

Variability 

5 

1-2 

0.5 

Internal waves 

3 

- 

0.2 

Geostrophy (Gulf-Stream) 

2 

- 

- 

Bottom wedge 

1 

- 

- 

Total error 

6 

1-2 

0.5 


Table 2.2: Error budget across the 36°N Atlantic section. The a priori uncertainties 
are provided on the total flux and in layers above and below 2000 db (Sv rms). The 
geostrophic error comes from the Gulf-Stream curvature at 36°N. 


equivalent to requiring a solution in which the large scales are more energetic than the 
small scales by choosing accordingly the covariance (bb T ). However, the spectrum varies 
substantially with location and a white noise assumption is the most practical prior choice. 

The results are dependent upon the skill of the GCM in representing the true ocean. 
On the one hand, the GCM velocity variability is lower than the true variability by a factor 
of 2 to 4 (Stammer et al., 1996) although we do not know what are the respective fractions 
of baroclinic and barotropic variability that are represented. On the other hand, the GCM 
is known to contain a drift term. The drift was not corrected and will introduce spurious 
variations in the density field (Marotzke, personal communication, 1998). Nevertheless, 
the uncertainties are in agreement with an analysis of the few sections that were repeated in 
the Atlantic within a factor of two, which is sufficiently accurate for an a priori uncertainty 
estimate. 

Another component of the model error is the Eulerian nature of the transports calcu¬ 
lation. Eulerian tracer transports are an approximation to the true Lagrangian property 
transports. This approximation may be a source of noise in tracer transport calculations 
and budgets. The problem is too vast to be addressed in the present context and I am 
left with the assumption that the error implied does not greatly exceed the error from the 
variability. 

The uncertainty attributable to a mass conservation equation involving several sec¬ 
tions depends upon the correlation between the errors in the flux across each section. 
Measurement errors are probably uncorrelated from one section to another. In the model, 
the correlation between the net relative transport through 24°N, 36°N and 48°N in the 
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GCM was calculated over the 8-year time series. The maximum correlation was found 
between 36°N and 48°N with a value of 0.24, which is not significant as the 95% non-zero 
correlation is 0.35. (Each time series of relative transport had a decorrelation scale of 
100 days, leaving only 30 independent samples.) A coherence analysis showed no signif¬ 
icant correlation at any time scale. Consequently, one can reasonably assume that the 
uncertainties are uncorrelated from one section to the next and sum the error variances of 
sections entering a mass conservation equation. 

Conversely, variations in different layers of a section are probably correlated. Roem- 
mich and Wunsch (1985) found systematic deviations in comparing data from sections 
taken years apart. This correlation is likely to depend upon many parameters like depth, 
presence of a boundary current in the layer, latitude, etc. For simplicity, and due to the 
lack of information on those correlations, it was decided to leave the a priori uncertainties 
uncorrelated between layers. 

2.6 Application to the global inversion 

In principle, with a global GCM, one could repeat the error budget analysis in each oceanic 
region. But to determine the a priori uncertainties that are needed, it was judged sufficient 
to adapt the results in the North Atlantic to other oceans and to low latitudes. 

Geostrophy has never been proved to fail either from scaling arguments or from ob¬ 
servations (Eriksen, 1982), even within a few degrees of the Equator. To extrapolate 
a general rule based upon the present study, the uncertainty from measurement noise 
(mainly internal waves) was scaled by the inverse of the Coriolis parameter. The measure¬ 
ment noise goes from ±2 Sv at 50°N to ±20 Sv at 5°N (Table 2.3, second column). Because 
of geostrophy, this transport error is independent of the distance between stations. 

The model error due to variability does not seem to depend directly upon the Coriolis 
parameter in any of the three North Atlantic sections in the GCM and the extrapolation 
at low latitudes were obtained by assuming, as a zero-order guess, that the variability had 
about the same contribution as measurement errors on the total budget (Table 2.3, third 
column). The variability in individual layers was scaled with the Coriolis parameter. The 
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Table 2.3: A priori uncertainty (Sv rms) in the mass conservation equation for a single 
section. The results from the GCM are extrapolated at low latitudes. 

baroclinic variability is potentially higher at low latitudes because of the fast baroclinic 
waves existing near the Equator, so that the variability in the wind field, which is associated 
with large Ekman transports, has a faster impact on the density field. Wijffels et al. (1996, 
Appendix B) used XBT time series to estimate baroclinic variability in the upper 400 m 
in the North Pacific Ocean at 10°N at roughly ±7 Sv, which agrees with my a priori value 
in the upper layers at 10°N (Table2.3). However, it will be verified, a posteriori, that this 
extrapolation to low latitudes is adequate. 

Table 2.3 was used to calculate the a priori uncertainties through each section, each 
layer and in each box. Those are listed in Appendix C, Table C-l and C-2. In the case of 
outcropping layers, the uncertainty contribution from the section was made proportional 
to the actual length of the lower interface of the layer. 

At 36°N in the Atlantic, the variability in the baroclinic silica fluxes was estimated 
to about ±50 kmolSis -1 . The a priori uncertainty on the total silica conservation in a 
North Atlantic box between two sections is thus of \/50 2 + 50 2 + 50 2 = 86 kmolSis -1 , 
where the third term (±50kmol Si s -1 ) accounts for net sources and sinks (discussed 
later in Section4.3.1). The total a priori silica uncertainty in the North Atlantic was 
thus taken to ASi N ' AtlantIC = ±100 kmolSis -1 . This top-to-bottom silica uncertainty was 


extrapolated to other latitudes and oceans according to the average concentration of silica 
and the local mass variability AMass local : 


, 0,0 N (Si local ) AMass locaI 
51 X (Si 36 ° N > X AMass 36 ° N ’ 


(2.19) 


where A is the a priori uncertainty and (Si) is the concentration. (For a given uncertainty 
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in the baroclinic mass flux AMass local , it is assumed that a larger average silica concen¬ 
tration implies larger silica differences over the water column and thus a larger variability 
in the silica flux). 

The a priori uncertainty in silica conservation equations in deep layers used by Mac¬ 
donald (1998) was of the same order as the one obtained from variability, that is, 16 to 
40 kmolSis 1 . On the other hand, Robbins and Toole (1997) did not address oceanic vari¬ 
ability and used a more stringent constraint for the silica budget in the Indian Ocean, 
(±100 kmols -1 compared with ±700 kmols -1 that is obtained from Eq.2.19 for a budget 
over the Indian Ocean, including the Pacific-Indonesian Throughflow section and the 32°S 
section). 

Salinity and phosphate-oxygen “PO” concentrations do not vary greatly with location, 
and their uncertainties were derived from the uncertainty in the mass flux, as described 
in the anomaly equation formulation, Section 1.2.4 (page 22). 
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Chapter 3 


Physical circulation 


The method described in Chapter 1 along with the a priori statistics from Chapter 2 (Sec¬ 
tion 2.6) are used to estimate the ocean circulation. In this chapter, the physical circula¬ 
tion is discussed individually for the Indian, Pacific, Atlantic and Southern Ocean with 
a focus on the meridional overturning circulation (MOC). (Here the expression “merid¬ 
ional overturning circulation” is used to denote the net zonally integrated stream function 
of the depth/meridional plane. There is no implication that it is restricted to a single 
top-to-bottom cell; the vertical structure is to be determined. Its strength will however, 
sometimes be measured by the net inflow of deep waters.) A summary of the global cir¬ 
culation follows, including global averages of dianeutral velocities and diffusivities. Each 
discussion starts with an introductory literature review (first section); it is followed by 
a presentation of the regional model constraints (second section), and a circulation de¬ 
scription (third section). Global features of the freshwater divergences, heat, nutrient and 
oxygen fluxes are discussed in the following chapter. 

The vertical structure of the net flux across the hydrographic sections is presented 
graphically as an overturning streamfunction which is defined as the cumulative sum of 
the transport in individual layers starting from the bottom. [This convention contrasts 
with some of the literature where the results are usually presented as individual layer 
fluxes.] But it has the advantage of making the useful quantities readily accessible to 
the eye, the sum of the fluxes over several layers being just the difference of value of the 
streamfunction between the first and last layer. The direction of the flow is thus given by 
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the slope with depth of the overturning stream function. 

The results are compared with the existing literature. However, when the error bars 
were not provided by the author, the comparison remains qualitative. I assume that their 
error bar is zero, and qualify the number as “compatible” when inside my error bars and 
as “incompatible” when outside. 

I start with a relatively extensive discussion about the circulation in the Indian Ocean 
as this work represents the first time that the high-quality WOCE data are used to estimate 
its basin-scale. Shorter presentations of the derived circulation in the Pacific, Atlantic and 
Southern Ocean follow, with a focus on basin-scale integrated transports. 

A number of solutions were explored by changing the model configuration. The “stan¬ 
dard”, or “preferred” solution employs property anomalies and the constraints listed in 
Table 1.2 (page 17), that is, top-to-bottom silica is conserved as described in Section2.6; 
salt is conserved top-to-bottom and within individual layers; heat, and the tracer combi¬ 
nation “PO”= I 7 OXPO 4 +O 2 (Redfield et al., 1963; Anderson and Sarmiento, 1994), are 
conserved only in layers below the surface. (The surface layer is non-conservative because 
of heat and oxygen exchanges with the atmosphere.) The freshwater flux was not con¬ 
strained, but was solved for through the inversion. Again the a priori uncertainties are 
consistent with the error budget of Section 2.6, taking into account the oceanic circulation 
variability. Conservation equations for heat, salt and PO are written as anomaly equa¬ 
tions, so that the mass divergence time the average property concentration is subtracted 
from the original conservation equation (see Chapter 1). Some constraints (e.g., “PO” in 
the South Atlantic) were given a lower weight when large mass residuals were obtained in 
the first-try solution. Alternative solutions are described in the discussion for each ocean. 

Although water mass distributions guided my choice of the “first guess” circulation, the 
flow was not explicitly constrained to force the net layer fluxes to be advected away from 
their supposed sources. In several places, the estimated flow, although indistinguishable 
from zero within error bars, does not follow the large-scale property distributions in the 
sense suggested by Wiist. Away from the core of boundary currents, evidence is relatively 
weak that the large-scale property distributions of a time-dependent flow require the mean 
advective flux to follow the water mass pathways. Here, near-conservation on the large- 
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scale is enforced, which is all that appears compelling to do when using such a model. 


3.1 The Indian Ocean 

3.1.1 Introduction 

My discussion starts, somewhat arbitrarily, with the Indian Ocean. The Indian Ocean 
has at least four unique features compared to all other oceans: 1) a strong monsoonal 
variability; 2) a southward average Ekman transport on both sides of the Equator; 3) a 
low latitude northern boundary with 4) a net influx of warm water from the Indonesian 
Throughflow. This section is also the subject of a publication (Ganachaud et al., 1999) 
that is focused on the meridional overturning circulation. 

Since the studies of Warren (1981) and Fu (1986), the meridional overturning circula¬ 
tion of the Indian Ocean has been the subject of a number of recent papers (Toole and 
Warren, 1993; Robbins and Toole, 1997 ; Lee and Marotzke, 1997, 1998; Macdonald, 1998; 
Zhang and Marotzke, 1999). Although the meridional overturning rate and structure are a 
basic descriptive feature of the general circulation as well as a determinant of the heat and 
other property fluxes, there remains significant disagreement both about its qualitative 
and quantitative features in the Indian Ocean. 

From a hydrographic section at 18°S, Warren (1981) estimated a northward flow of 
19 Sv in the deep western boundary currents. Fu(1986), using the same section, but 
combining it with four other sections in the Indian Ocean and doing a formal inverse 
calculation, found weak net deep water inflow of 3.6 Sv, with the bulk of the net northward 
flow occurring at intermediate levels (1000-2000 decibars). In contrast, Macdonald (1998), 
again with the same 18°S section but combined with a global data set containing more 
recent sections than were available to Fu, (notably at 32°S), obtained 10 Sv of deep water 
northward flow at 18°S with a 5 Sv southward flow at intermediate levels. 

Using the pre-WOCE 1987 section at 32°S, and tracer properties to suggest a zero- 
velocity surface, Toole and Warren (1993) estimated the surprisingly large bottom water 
flux below 2000 decibars of 27±10Sv geostrophically. Robbins and Toole (1997) recog¬ 
nized that this large flux would cause problems with the silica budget, and added silica 
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conservation as a constraint. They obtained, through an inverse calculation, a weaker 
deep influx of 12±3 Sv. At the two standard deviation level, this value is inconsistent with 
the non-silica-conserving Toole and Warren (1993)-type circulation, 29±5Sv, according 
to the Robbins and Toole (1997) error analysis. Assuming the error estimates are ap¬ 
propriately computed, the inconsistency implies systematic (non-random) errors present 
in the models. In both circulations, the intermediate waters (1000-2000 decibars) moved 
southward. Bryden and Beal (personal communication, 1999) showed that the inferred 
overturning was sensitive to the assumed strength of the Agulhas Current, and notably to 
the presence of a northward flowing undercurrent as detected by ADCP measurements. 

Lee and Marotzke (1997, 1998), using a GCM, climatological hydrography and surface 
fluxes, inferred a mean northward overturning flow of 14 Sv in the upper 1000 m with 
a weak south-to-north deep water flow. Their solution is consistent with several other 
GCM calculations that exhibit weak deep inflow, but a vigorous shallow overturning (see 
Zhang and Marotzke, 1999, for a review). Such a shallow overturning conflicts with the 
direction of the flow as would be suggested by a classical water mass analysis (e.g., Toole 
and Warren, 1993). 

3.1.2 Indian Ocean Problems 

The Indian Ocean is complex, and a full discussion of the details of the circulation as 
previously described would be very lengthy. Here I only summarize some of the major 
uncertainties most directly affecting the inversion. 

Monsoonal Effects 

Because of the intense monsoonal shifts, estimation of the time-average circulation from 
synoptic sections is problematic; my main excuse for nonetheless doing so is that below 
the top few hundred meters, where my uncertainties are large, and away from the Somali 
Current, there is no evidence of a significant annual shift in hydrography or baroclinic flow 
field (e.g., Warren and Johnson, 1992, Fieux et al., 1986). (Fieux et al. (1986) found a 
variation of 30% in the deep western boundary current transport at 10°N, which may be 
due to internal wave noise (Section 2.3.1)). 
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The most spectacular time-dependent feature is the establishment and reversal of the 
Somali Current whose net transport at the Equator ranges from 0 to 21 Sv (including 
undercurrents) (Schott et al., 1990). However, the GCM of Lee and Marotzke (1998) sug¬ 
gested that the upper ocean variability decreases toward the south (away from the Mozam¬ 
bique Channel). They estimated a monsoonal variability in the depth-latitude stream- 
function of ±20 Sv with associated changes in the meridional energy flux of ±1.8 PW. But 
most of the streamfunction variations were explained by a barotropic response to Ekman 
transport variations, with a relatively small contribution from the vertical shear (±2 Sv or 
less). (The Somali Current region was exceptional; there the vertical shear contribution 
was large, of order 10 Sv.) Barotropic fluctuations are not aliased in the hydrography 
so that a model using the mean Ekman flux can be combined consistently with sections 
taken at different seasons. In the present analysis, internal waves and variability in the 
baroclinic field are accounted for in the uncertainties (Chapter 2). 

Non-linear interactions of the time-dependent depth-independent transports could con¬ 
ceivably generate significant rectified mean flows. Lee and Marotzke (1998) could, however, 
find no such effect in the energy. 

Agulhas Current and Mozambique Channel 

The gyral subtropical circulation is closed by the Agulhas Current, believed to transport 
between 66 and 70 Sv (Beal and Bryden, 1997). Stramma and Lutjeharms (1997) suggest 
that most of the Agulhas waters are fed by the subtropical gyre (e.g. their Figure 7), but 
with a fraction (20-25 Sv) from either the Mozambique Channel or the East Madagas¬ 
car Current. Of these Stramma and Lutjeharms (1997) and Saetre and DaSilva (1984) 
suggested the latter was the more important, in contrast with Lutjeharms (1976) who con¬ 
cluded that the Mozambique Current is the major source of Agulhas waters. Saetre and 
Da Silva (1984) found no continuous upper ocean current in the Mozambique Channel and 
little contribution to the Agulhas Current. This inference was in agreement with the study 
of Menache (1963) who believed that the Mozambique Current curved eastward and flowed 
back to the north along the west coast of Madagascar. Recent CTD data and satellite 
imagery showed that the Agulhas Current is partly fed by the Mozambique Current, and 
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partly by a westward flow south of Madagascar (Griindlingh, 1993). 

Fu’s (1986) inverse box model produced a net southward flux of 6Sv in the Mozam¬ 
bique Channel, similar to that of Stramma and Lutjeharms (1997). The global inversion 
of Macdonald (1998) did not derive a significant flux there. 

The flow through the Mozambique Channel appears to be highly variable on seasonal 
time scales. Measurements from tide gauges combined with hydrography (Donguy and 
Piton, 1991) show an annual cycle with ±10 Sv amplitude in the Channel. Altimetric data 
(not shown) examined by Marotzke (personal communication, 1999) show fluctuations 
of ±25 cm in the sea surface height difference between South Madagascar and Africa 
(25°S). Assuming a linear decrease of the velocity down to 1000 decibars, the implied 
variations in the net geostrophic transport are ±20 Sv (This calculation gives only an 
order of magnitude estimate, but the 1000 decibar vertical scale and the linear decrease 
are consistent with the solution described below. Baroclinic Kelvin waves may account for 
these rapid fluctuations. Seasonal variability is probably also present in the density field, 
as observed in the GCM of Lee and Marotzke (1998) in the Somali region.) 

Pacific-Indonesian Throughflow (PIT) 

The Pacific-Indonesian Throughflow (hereafter PIT) is highly variable and poorly sampled, 
but a potentially important element of the Indian and Pacific oceans circulation. If the 
PIT average transport were 10 Sv at 20°C and assuming that it ultimately leaves the 
Indian Ocean at 10°C, it would produce a convergence of 0.4 PW in the Indian Ocean heat 
transport budget. This flux convergence is of the same size as the total climatological heat 
input to the Indian Ocean by air-sea exchanges from Oberhuber (1988). 

The published estimates of the net PIT mass transport range from 18±7 Sv westward 
to 2.6 Sv eastward (Fieux et ah, 1996), with a host of estimates between these two ex¬ 
tremes, depending upon the measurement time and method (Godfrey, 1996). Macdonald 
(1998) constrained her inverse model to a net PIT of 10±10Sv and obtained, after inver¬ 
sion, a flux of 9±7 Sv. If not specifically constrained, her model produced a net PIT of 
ll±14Sv. Interestingly, her sensitivity experiments (PIT= OSv; PIT= 20 Sv) also showed 
that variations in the PIT strength had no effect on her estimated Agulhas Current and 
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Drake Passage fluxes. Nevertheless, such changes affect dramatically property budgets and 
derived air-sea exchanges in the Indian and Pacific oceans. Zhang and Marotzke (1999) 
estimated the PIT from climatological Indian Ocean hydrography and a GCM, and found 
that a small mean inflow (2.7 Sv) was sufficient to close the Indian Ocean heat and salt 
balance. 

The strength of the PIT is plausibly dependent upon the phase of ENSO (e.g., Waj- 
sowicz, 1994, 1997, who used a numerical model). Prom XBT time series and historical 
data, Meyers et al. (1995) found a total mean transport in the upper 400 m of 5Sv to 
the west with a ±2.5 Sv variation in the annual cycle. They estimated a top-to-bottom 
mean transport of 7Sv. Moorings in the Makassar Strait, which is believed to be the 
main path of the Throughflow, showed an annual average transport of 9 Sv in 1997 with 
variations from 5 Sv during El Nino to 12 Sv during La Nina (Gordon and Susanto, 1998). 
For present purposes, I use a Throughflow of 7±7Sv as the initial best estimate. (But to 
account for variability and measurement noise from my error budget, the actual constraint 
on the PIT between Indonesia and Australia is 7±15Sv, as explained below.) 

3.1.3 Model constraints 

Figure 3-1 shows the Indian Ocean sector of the model which is made of the WOCE sections 
12 (Dec. 1995-Jan. 1996, e.g., Johnson et al.,1998); 13 (Apr.-June 1995); 14 (June 1995); 
110 (Nov.-Dee. 1995, Bray et al. (1997); the pre-WOCE section 15 (Nov. 1987; Toole and 
Warren, 1993; Robbins and Toole, 1997) and the Java Australia Dynamics Experiment 
(JADE) section (Aug. 1989; Fieux et al., 1994; Coatanoan et al., 1999). (The 1989 JADE 
data were used instead of the 1992 data because they included nutrient measurements. 
However, the 1989 data were combined with a station of the 1992 data to complete the 
section, as discussed below (Section 3.1.6.)) 12, 13 and 14 were obtained from the WHPO 
data center with permission of the respective Pis B. Warren, G. Johnson, W. Nowlin and 
J. Toole. 110 was provided by G. Packard with permission of the Pis J. Toole, N. Bray and 
J. Sprintall. The data from the Franco-Indonesian JADE cooperation were provided by 
M. Fieux and A. G. Ilahude (Pis); B. Coste and C. Coatanoan (nutrients). 

Some unrealistically large geostrophic velocities were found in the thermocline waters 
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of 12 due to narrow station spacing and small Coriolis parameter. An attempt was made 
to apply a horizontal filter to the temperature and salinity fields in the thermocline waters. 
The filter wavenumber cut-off was set so that large velocities (>2 m/s) were removed. The 
filtering did not change the integrated mass transports so I decided to leave the original 
temperature and salinity fields to avoid the arbitrary treatment of edges in the filtering 
process. Some velocities are unrealistic when station spacing is narrow, but the implied 
transport is accounted for in the error budget. 

Constrained boxes are defined for the Indian Ocean as the “subtropical box” (15, 14, 
13), the “Mozambique Channel box” (14, I2W), the “Tropical box” (13, 12, 110) and the 
“North Indian box” (I2W, 12, 110, J89). Note that the North Indian box includes the area 
between 110 and J89 because the 110 section could not be completed up to the Indonesian 
coast. The layers, defined by neutral surfaces (Table3.1), were chosen to facilitate the 
comparison with the existing literature, to match the Southern Ocean layers of the global 
model, and to have a relatively homogeneous thickness over the basin. Mass conservation 
was required consistent with the error budget of Chapter 2 (see TableC-2, page 194). The 
resulting a priori uncertainties that characterize the net transports range from ±7Sv at 
32°S to ±15 Sv on the JADE 89 PIT section. The Ekman transports from the NCEP 
reanalysis (Kalnay et ah, 1996) with their a priori uncertainties are listed in Table 3.2. 

Because of the large baroclinic variability in the East African Coastal Current (EACC) 
on section I2W, a large adjustment (±20 Sv) was allowed to the flow in the first layer. 

The Indian Ocean topography is complex (Figure 3-1) and the deep layer fluxes were 
broken down and constrained according to the bathymetric features. These “bathymetric 
constraints” are listed in Table 3.3. Unlike the RT97 inversion, the net bottom water flow 
into the Crozet and Perth Basins is not constrained to be northward. This change is 
discussed further below. The flux into the Mozambique basin was weakly constrained to 
be northward (l±2Sv). The flux below 1300 decibars in the PIT was constrained to be 
small as this level corresponds to the limit of the waters of Indonesian origins (Fieux et 
al., 1994). 
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Figure 3-1: Indian Ocean sections and topography. The 4000 m isobath is contoured and 
areas shallower than 3000 m are shaded. Box names are indicated in the shaded rectangles. 
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Layer 

Upper 

Interface 

RT97 

Layer 

15 Interface Property 
from RT97 

Water Mass 

Name 

1 

Surface 

1 


Surface Water 

2 

7 n =25 

3 

7 n =26.2 

2 


Subantarctic Mode Water 

4 

7 n =26.9 

3 

o e = 26.80; 9 = 9.85 

Upper Antarctic Mode Water 

5 

7 n =27.36 

4 

o 6 = 27.20; 9 = 4.65 

Lower Antarctic Mode Water 

6 

7 n =27.7 

5 

<72 = 36.65; 9 = 3.26 

Upper Deep Water 

7 

7 " =27.96 

6 

<r 2 = 36.92; 9 = 2.20 

Lower Deep Water 

8 

7 n =28.07 

9 

7 n =28.11 

7 

<74 = 45.89; 9 = 1.28 

Upper Bottom Water 

10 

7 n =28.15 

11 

7 n =28.23 

8 

<74 = 46.01; 9 = 0.43 

Lower Bottom Water 


Table 3.1: Indian Ocean layers ( 7 " upper interface in kg m -3 ). For comparison, the 
Robbins and Toole (1997) layers are indicated along with the potential density and tem¬ 
peratures at each interface (from their table 1). The italicized layer numbers separated 
by the line indicate a slightly different interface ( 7 " = 26.20 instead of 26.50 in RT97). 


Section 

Initial 

Final 

% correction 

15 

0.3±1 

0.5±1 

- 

14 

-0.5±1 

-0.5±1 

- 

13 

-11.5±5 

—10±4 

- 8 % 

I2W 

-4.5±20 

9±10 

-300% 

12 

-11.6±5 

-11±5 

- 6 % 

no 

-2.9±1 

-2.9±1 

0 % 

J89 

-0.65±1 

- 0 . 6±1 

- 


Table 3.2: Initial Ekman transport and that derived from the model (Final). The adjust¬ 
ments apply to the total transport in the first layers, including in principle any deviation 
from the time-mean geostrophic transport (e.g., monsoonal variability). The relative am¬ 
plitude of the adjustment is indicated for transports larger than 1 Sv. 
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Name 

Section 

Layer Depth 


Longitude 

Natal Valley 

15 

9-10 (3400) 

0±0.5 

30°E-34°E 

Mozamb. B. 

15 

9-11 (3500) 

1±2 

35°E-46°E 

Madagascar B. 

15 

9-11 (3600) 

0±5 

46°E-56°E 

Central B. 

15 

9-11 (4000) 

Oil 

78°E-89°E 

Mozamb. Sill depth 

14 

8-11 (2800) 

0±0.5 

Full section 

PIT sill, mass 

J89 

6-10 (1300) 

0±6 

Full section 



each layer 

0±3 


silica 


total 

0±140 

kmolSis -1 



each layer 

0±70 

kmolSis -1 

Name 

Sections 

Layer Depth 


Longitude 

Mozamb. B. 

15/14 

8-11 (2800) 


35° E-46° E/Full 14 

Perth Basin 


9-11 (3500) 


96° E-end/87° E-end 

Mascarene B. 

PITT 

8-11 (3000) 


49 0 E-58°E/49 o E-59°E 

West Aust. B. 

Mrsmtm 


0±4 

96° E-end / 89° E-end 


Table 3.3: “Bathymetric” constraints on the flux in layers. The constraints are imposed 
on either the total flux through a single section (upper part, positive northward) or on the 
divergence between two sections (lower part). The layer and depth below which the flux 
is constrained are indicated along with the uncertainty on the net flux (single section) or 
divergence (between sections) (in Sv). “B.” denotes “basin”. 

Reference Surfaces and Dianeutral terms 

The initial zero-velocity-surfaces, or “reference surfaces”, are indicated in Table3.4. The 
15 reference surface was taken from the RT97/TW93 initial guess with a modification 
close to the coast in the Agulhas Current region to account for the presence of an Agul¬ 
has Undercurrent (Beal and Bryden, 1997, Fig. 3 ). In general, a deep reference surface 
was preferred to a shallow one to decrease the potentially large velocity errors at depth, 
which can result in an eddy-induced shear region. The northward moving bottom waters 
in the abyssal Indian Ocean lie in the potential temperature range -0.8 to 1.4° C in the 
western basin and 0 to 1°C in the eastern basin (Schmitz, 1996b, p.122 ). Therefore, 
the interface between the Circumpolar Deep Water and the “Indian Ocean Deep Water” 
defined by r y n = 28.11 (9 = 1.2° at 20°S) was chosen as a “first guess” for the 12, 13 and 
110 reference surface, (a similar surface was chosen by other authors, e.g. Johnson et al., 
(1998, Table 1.)) The model is not crucially dependent upon this choice: an experiment 
which started with a bottom reference surface at 12, 13, 110 and 15 (away from the Ag- 
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Section 

Reference Surface 

Remark 

15 

station-specific 

TW93/Beal and Bryden(1997) 

14 

2500db 

Mozambique sill depth 

I2W 

bottom 

See text 

12,13,110 

7 n = 28.11 

RT97/TW93 

J89 

bottom 

See text 


Table 3.4: Initial zero-velocity-surface, or reference surfaces. 


ulhas Undercurrent) produced a similar overturning after inversion. Strong shear down 
to 3000 decibars in the composite JADE section led to the choice of a bottom reference 
level, with the bathymetric constraint ensuring no net flow in individual layers below the 
sill depth ( 7 71 = 27.7, 1300 decibars) by adjusting the reference velocities accordingly. 

In common with previous inverse solutions, the small scales of the reference velocity (&;) 
are not resolved, and the posterior uncertainties on the individual 6 * are almost unchanged 
from the prior values. Following the error budget (Section2.5), the a priori uncertainties 
were given a minimum value of ±1 cms -1 and increased according to the local shear and 
when current-meter data suggested it (Dickson, 1989, page E-9, E-ll ; Schott et al., 
1988). 

A large range of adjustment was allowed to dianeutral transfers w* (±10~ 3 cms _1 ) and 
mixing coefficients k* (±100c.m 2 s _1 ). These a priori values are near the upper limit of, 
or larger than, published estimates (e.g. TW93, RT97; Polzin et ah, 1997; Munk and 
Wunsch, 1998), leaving the model free to determine the coefficients from the constraints. 
However, both upwelling and mixing turn out to be well-determined by heat and salt 
anomaly conservation requirements, and the final values have smaller amplitudes than the 
a priori ranges. No positivity requirement was set for k*. 

3.1.4 Indian Ocean Circulation: standard solution 

The resolution matrix of the Gauss-Markov estimate shows that heat and salt anomaly 
conservation in the layers are well resolved, i.e., provide independent information. The 
anomaly equations for the conservative tracer “NO” (9.1 xN0 3 +0 2 ; Broecker (1974)) are, 
in general, indistinguishable form the “PO” equations and I choose to conserve “PO” 
because the nitrogen cycle is the more complex one. Mass conservation equations are only 
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a) Subtropical 32S-20S 


b) Tropical 20S-8S 


c) North (8S-JADE 89-Coast) 


o 

1000 

2000 

3000 

4000 

5000 

6000 


-5 



-1.7±5.9x 10 kg/s 


-0.033±5.8x 10 s kg/s 


3.8± 8x 10 kg/s 


Figure 3-2: Standard solution mass residuals in the Indian Ocean in a) the subtropical box; 
b) the tropical box; and c) the North Indian box, which also includes the area between 110 
and JADE 89. The line gives the residual in invidual layers (in 10 9 kgs _1 ) and the shaded 
area give the one standard deviation uncertainty. The top-to-bottom residual is indicated 
below the graph. 


weakly resolved. The solution elements (reference velocities; Ekman transport; freshwater 
flux; dianeutral mass transfers and mixing) were all within the a priori range. Except for 
silica, all constraints were met within uncertainties, and the solution appears to be both 
dynamically and statistically acceptable. Dissolved silica in the Indian Basin is conserved 
only at two standard deviations, with a net convergence of —620±360 kmol Si s , possibly 
associated with the high PIT flux (there is an equivalent divergence in the central Pacific, 
see Section4.3.1). Mass residuals in individual layers are smaller than 1 Sv and all are 
indistinguishable from zero within two standard deviations (Fig. 3-2). The net residual in 
the Mozambique box (not shown) is also small at 0 ± 7 Sv. The 13 and 110 intersection 
near Australia (Fig. 3-1) generates a small closed box. Experience shows such boxes should 
not be constrained, owing to the potentially large hydrographic aliases at open-ocean 
intersections. Nevertheless, the residuals in this small area were confirmed, a posteriori, 
to be small. 

Small adjustments were made to the initial NCEP Ekman transports (0% to 8%, Ta¬ 
ble 3.2), except in I2W (300%) in the Mozambique Channel where the variability in the 
baroclinic transports is large. (The adjustment is not attributable solely to the Ekman 
transport, as it includes everything accommodated in the first layer(s) to meet the con- 
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servation constraints.) 

The reference velocities, absolute velocities and transports cumulated from West (or 
from North along JADE 89, I2W and 110) from the standard solution are given in Ap¬ 
pendix D (pages 196-209) but my discussion will focus on the net integrated fluxes and 
divergences. 

Net Fluxes 

The integrated flow for the standard solution, or net “overturning streamfunction” and the 
net flux across each section are shown on Figure 3-3. (Again, overturning streamfunctions 
are used, the slope of the curve indicating the direction of the flow, positive to the north 
for zonal sections, and to the east for JADE). The uncertainties, calculated from the error 
covariance of the solution vector, are given by the shaded area. 

The net meridional mass fluxes, which include the Ekman transports, are interpreted 
as equal to the PIT flux and are given below each graph of Figure 3-3. These net mass flux 
estimates are the same, across each latitude within error bars. At 32°S, the PIT flux is 
estimated as 16±5Sv (Fig. 3-3a). Through J89, the net flux is 15±3Sv westward (Fig. 3- 
3f), in the upper range of the a priori constraint on the PIT. (The final uncertainty of the 
PIT is surprisingly small, given the prior uncertainty. Its size results from the bathymetric 
constraints below the sill depth (Table 3.3). An uncertainty of ±5 Sv is believed to be more 
realistic, as discussed later). In the Mozambique Channel, the flux is southward and of 
the same magnitude as the PIT, 14 ± 6 Sv (14 and I2W, Fig. 3-3d, e), implying a weak and 
uncertain net transport between Madagascar and Australia of —4 ± 8 Sv (13, not shown). 
The net transport through 110 is westward but very uncertain, 2± 14 Sv, so that the mass 
flux between Java and the 12-110 intersection is not known. 

Deep Inflow 

The estimated Indian Ocean overturning, or deep inflow below 2000dbars at 32°S is 
10.6 ±4Sv (15, Fig.3-3a). Similar values are obtained at 20°S (i.e., 14+13) and 8°S 
(I2W+I2+I10), although the vertical structure changes due to dianeutral exchanges (de¬ 
scribed below). The net overturning (i.e., deep inflow) is about the same magnitude as 
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a) 15 (32S) 


b) 13+14 (20S) 


c) I2W+I2+I10 (8S) 



-16.3± 5.1 x 10 9 kg/s 
d) 14 (Mozambique S) 




-18.3+ 6.4x 10 9 kg/s 
e) I2W (Mozambique N) 


—14.1 ± 6.3x 10 9 kg/s 




f) JADE 89 (116E) 



-14.8± 3x 10 kg/s 


Figure 3-3: Overturning streamfunction for the standard solution across latitudes a) 32°S; 
b) 20°S; c)8°S; d) across the southern and e) northern Mozambique Channel; f) across the 
PIT. Each step of the streamfunction gives the direction and magnitude of the flow in the 
layer above (positive to the north/east). The layer interfaces are indicated to the right of 
each graph. 
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the one obtained by RT97, the flow has a different structure. In my solution, the Agulhas 
Undercurrent carries 13 ±5 Sv to the north in layers 6 to 7 west of 31°E with an immediate 
southward recirculation to its east, producing an insignificant transport of 3 ± 6 Sv to the 
north between the coast of Africa and 35°E in those layers. This feature was not present 
in the previous inversions. Due to the lack of resolution of the model short scales, the hor¬ 
izontal partition of the flow is not accurately known. However, in the Madagascar-Crozet 
basin and below 2000 dbars, there is no indication of an anticyclonic circulation as in RT97 
(their Fig. 3), but a quasi-uniform northward flow east of 52°E (Fig.D-2, page 197). In 
the Perth basin, there is a weak and uncertain anticyclonic circulation and the net flux 
below 2000 dbars is indistinguishable from zero (l±8Sv to the north). 

The vertical flow structure differs with the RT97 as well (Fig. 3-3). Below j n = 28.11, 
my northward flux is of 8 db 3.5 Sv, whereas the deep inflow in the RT97 solution below 
this surface was 12 i 3 Sv. But in my standard solution there is 2.6 Sv of northward flow 
just above = 28.11, so that the total deep inflow below 2000 dbars is similar to the 
that of RT97. RT97 constrained the deep flow to move away from its sources in individual 
basins, resulting in an equivalent constraint on the Circumpolar Deep Water (CDW) flux 
of 18±5Sv below q 71 = 28.11 (their Table 3, with the errors computed by assuming the 
constraints in the different basins are independent). The magnitudes and weights of these 
constraints were imposed subjectively (P. Robbins, personal communication, 1998) and 
account for their larger resultant inflow at depths. 

At 32°S, but in the bottom layers alone, (layers 9-11, 3200 dbars to bottom, not shown) 
the flux indicates northward-directed boundary currents along the ridges. The net flux 
across the Indian basin and below 3200 dbars is 8 ± 4 Sv to the north. Half of this occurs 
in the Madagascar-Crozet basins while the other half occurs in the Perth Basin. But the 
character of the horizontal structure of the deep flow is uncertain dependent upon depth. 

At 20° S between Madagascar and Australia, the flux in layers below 2000 dbars, 11 ± 
4Sv, occurs mainly in the Mascarene Basin (Fig. D-6, 7±6Sv northward-of which 4±5Sv 
occurs west of Mauritius (57°E)). In both the Central Indian and West Australian basins, 
the flows are noisy and uncertain (respectively 0.5 ± 7Sv and 3 ± 7Sv northward below 
2000 dbars), reflecting again the lack of horizontal resolution. Farther north, the fluxes in 
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individual basins become more and more uncertain (Figs D-8-D-12). 

The total deep water inflow into the Indian Ocean, 10.6 ± 4 Sv is thus consistent with 
the RT97 result (12±3 Sv) and larger than the GCM results. This overturning is the sum 
of larger local horizontal fluxes that partially cancel, e.g., those in the Agulhas Under¬ 
current; a 15-Sv eddy near 39°E, and the Madagascar-Crozet Basin circulation (Fig. 197). 
The present model would produce an uncertainty of ±8 Sv with total mass and silica con¬ 
servation only (cf., RT97). Addition of the other sections and constraints (heat, layers, 
etc.) accounts for the reduced uncertainty. 

Intermediate and Surface Flows 

The cumulative flow from the west across 32°S and above 2000 dbars (layers 1-6) is shown in 
Fig. D-2(upper). In the surface layers, at 32°S (15), the Agulhas Current carries 74±7Sv 
to the south in layers 1 to 5 west of 32°E, which is thus consistent with the Beal and 
Bryden (1997) estimate of 69 ± 2 Sv. No obvious horizontal recirculation is observed, 
consistent with the RT97 results; the gyre northward wind-driven flow occurs over the 
whole interior. 

At 25°S in the Mozambique Channel (14, Fig. D-4d), a strong anticyclonic circulation 
results in a powerful southward current of 23 ± 3 Sv at the western boundary in layers 1 
to 5 between 36°E and 37°E. 

Dianeutral Transfers and Mixing 

Dianeutral mass transfers (Fig. 3-4 upper panels) indicate downwelling at the base of the 
surface layers of the subtropical and North Indian boxes. Most of the deep upwelling ap¬ 
pears to take place in the North Indian box, with about 10 ± 5 Sv (3 ± 1.5 x 10~ 5 cms _1 ). 
The deep dianeutral transfers in the subtropical and tropical boxes are not significantly dif¬ 
ferent from zero. The upwelling in the North is associated net dianeutral mixing equivalent 
to diffusivity values of 4 to 10 cm 2 s” 1 with larger values in the deepest layer approaching 
30cm 2 s _1 . (Owing to possible large errors in the horizontal area of neutral surfaces at 
depths, this latter value is very uncertain). In this box, downgradient mixing is found 
consistently over all the deep interfaces, although there was no requirement that the mix- 
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Figure 3-4: Dianeutral transfers (upper) and diffusivities (lower) in the Indian Ocean. The 
line indicates the transfer or diffusivity at the bottom of the layer. 


ing coefficient be positive. The average diffusivity at interfaces below layer 5 and above 
the last interface is of 8.6 ± 4cm 2 s~* in the North Indian box. This value is much larger 
than the 1 cm 2 s _1 global mixing required by a global 1-dimensional balance (Munk and 
Wunsch, 1998). The North Indian Ocean topography is relatively rough, which can lead to 
enhanced mixing, as suggested by microstructure measurements (e.g., Polzin et al., 1997). 
(In the Mozambique Channel, dianeutral exchanges are not significant owing to the large 
seasonal baroclinic variability there.) 
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Standard circulation summary 

The standard circulation is thus associated with an Agulhas Current of 74 ± 7 Sv and with 
a relatively strong PIT of 15 ± 5 Sv. The overturning rate that I obtain is similar to the 
RT97 estimate. It is significantly larger than the Fu (1986) estimate (3.6 Sv) and smaller 
than the estimate of Macdonald (1998, 17 ± 5Sv). On the other hand, a much weaker 
inflow is found consistently in GCM simulations (Zhang and Marotzke, 1999). But GCMs 
also have their limitations (weak vertical mixing, dependence upon climatological forcing, 
and low resolution of the deep boundary currents), and so their systematic error could 
easily explain my more vigorous result. Nevertheless, to test the hypothesis of small deep 
inflow, I tried to constrain the flow below 2000 dbars to 0 ± 2 Sv. The model produced a 
2±2 Sv northward flow below 2000 dbars by returning part of its 4±3 Sv northward bottom 
inflow between 2500 and 3500 dbars, but the shallow meridional overturning of Lee and 
Marotzke (1997,1998) was not found. There was no major violation in the conservation 
equations in the Indian Ocean. However the dianeutral transfers were found significantly 
different from zero and oriented downwards in the subtropical and tropical regions, asso¬ 
ciated with negative diffusivities. Consequently, this solution was not retained. 

In the standard solution, dianeutral advective and diffusive transfers are consistently 
determined and have magnitudes in the higher end of the previously published range. 

3.1.5 The Mozambique channel 

Obtaining a time-average circulation in the Mozambique channel is a challenge because 
the 14 and I2W sections were taken at different seasons and the I2W section is close to 
the Equator (4°S), enhancing the noise in the geostrophic balance. Initial inversions 
produced a large and significant convergence in the intermediate and surface layers of the 
Mozambique Channel box. This convergence was found to be insensitive to the reference 
surface: the model adjusted the reference velocities to produce the same integrated flow 
for an initial sill-depth reference (2500 dbars) and an initial bottom reference. Because 
the convergence was significant, the original model was rejected. The error budget that 
I adapted from the Atlantic is too optimistic for this seasonally variable region. One 
could increase the a priori noise to account for the large variability in this region. The 
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alternative that I adopted was to allow a large ageostrophic adjustment in the shallowest 
layer of I2W to accommodate the variability. This adjustment formally appears as the 
“Ekman” transport adjustment, Table 3.2. 

After adjustment, mass is balanced within lSv in individual layers (not shown). Be¬ 
cause of the high variability, the vertical structure of the flow in the first layers is probably 
unrealistic. Nevertheless, the 14 ± 6 Sv net flux through the Channel (Fig.3-3d, e) is rel¬ 
atively robust and consistent within error bars with the different experiments that have 
been done in this region. DiMarco et al. (1998) obtained a bit larger southward trans¬ 
port (18 to 21 Sv) assuming a zero velocity at the sill depth and below. But the relative 
transport below 2500dbars was not included in their calculation (S. F. DiMarco, personal 
communication, 1998). To produce no flow below sill depth on I2W, important reference 
velocity adjustments were required, resulting in the smaller transport that I report. 

3.1.6 The Pacific Indian Throughflow (PIT) 

The mass flux through the PIT (15 ± 3 Sv) is on the high side of the commonly accepted 
values (e.g., Meyers et al., 1995; Godfrey, 1996), but consistent with the preferred Fieux 
et al. (1994) value from the same JADE data. The silica flux through the JADE section 
is larger by a factor of 2 (400 ± 200kmolSis _1 ) than what was estimated by RT97 from 
climatological silica concentration in the PIT region and a throughflow of 7 Sv. The section 
was constructed by adding a station from the 1992 JADE section to the 1989 JADE data 
to close the box at the northern end. Such an addition is artificial and may introduce a 
bias because of the large seasonal variability in the Java Current. Sensitivity experiments 
revealed however that the presence of the added station was essential to permit enough 
range of adjustment to the model (a large a priori uncertainty and reference velocity 
adjustment were allowed at the added station pair). In particular, removing this station 
produced large silica fluxes below the PIT sill depth and increased the net silica flux 
through the section. 

The PIT (Fig. 3-3f) shows a weak, uncertain net mass flux of intermediate waters 
flowing eastward (3 Sv, layers 3-4, 300 to 800 dbars) in those layers. There is an apparent 
contradiction with the tongue of silica emanating from the PIT in the Wyrtki (1970) atlas. 
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Fieux et al. (1994) discuss this eastward intermediate flux from those same 1989 data and 
find no inconsistency with the local water properties, the net flux being the sum of large 
and opposite fluxes of different waters types. (But as already described above, property 
tongues do not necessarily imply net mass flux. Any zero-mean time-dependent mass flux 
between the Pacific and the Indian Oceans which leaves behind fluid particles on the return 
flow would generate such a tongue (see Wunsch, 1996, p.79)). In any case, the silica flux in 
these layers is indistinguishable from zero. This intermediate water flux is enhanced when 
one removes the 1992 JADE station near Bali, and forcing the net intermediate silica flux 
to be westward causes large eastward silica and mass fluxes below sill depth on the JADE 
section, an unrealistic circulation. It was thus decided to leave the intermediate silica flux 
unconstrained. 

The local formal uncertainty of ±3 Sv on the PIT is probably optimistic because of 
the large baroclinic variability there (e.g., Meyers et al., 1995). A more realistic PIT 
uncertainty is given by the net mass flux uncertainty at 32°S, that is, ±5Sv. 

3.1.7 Full Property Equations and PO constraints 

A number of experiments were performed with the full property equations (that is, salt, 
heat and PO not expressed in terms of property anomalies). Use of full property equations 
has several drawbacks: The system is of lower rank, and it is difficult to weight the 
heat conservation equation properly. An uncertainty of ±0.1 PW in the heat equation 
is plausible in the thermocline; it is an overestimate in the deep layers characterized by 
smaller temperature gradients. In the anomaly equation formulation, the uncertainty can 
be attributed more consistently (Section 1.2.4), but the solution is, as always, directly 
dependent upon the reliability of the estimated error structure of the equations. 

The solution with the full property equations (not shown) indicates a similar deep 
water inflow at 32°S but with increased transport 15-18 Sv as one goes north. For the 
practical reason mentioned above, heat conservation was not required. To balance mass, 
the model creates downward dianeutral fluxes in the subtropical and tropical boxes, with a 
larger upwelling occurring in the North Indian box. Several experiments showed that with 
the full property equations, the model could not produce a consistent circulation without 
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those downward dianeutral fluxes. Although the circulations were not statistically different 
from the standard solution, those solutions were discarded. 

An experiment with anomaly equations (including heat) but no PO constraints indi¬ 
cated a similar circulation with downward fluxes in the subtropical and tropical boxes. 

3.1.8 Indian Ocean summary 

Figure 3-3 shows my best-estimate of the net meridional (a function of depth and latitude) 
overturning streamfunction of the Indian Ocean as determined from the WOCE hydro- 
graphic sections both locally and globally. This estimate is consistent with both dynamics 
and simple biogeochemical principles (silica and “PO” conservation). Deep upwelling is 
returned primarily in the 800-1500 dbars depth range, rather than near-surface. 

1. A meridional overturning strength (measured by the rate of deep inflow) of 11 ±4 Sv 
at all three latitudes (32°S, 20°S and 8°S) is consistently obtained, compatible with 
the RT97 estimate at 32°S. Nevertheless, the global circulation requires modifications 
in the zonal and vertical partition of the deep flow found by RT97. Results from 
general circulation models have consistently found a weak deep inflow (Zhang and 
Marotzke, 1999) which is barely compatible, within an uncertainty of two standard 
deviations, with my results. 

2. My estimate of the average Pacific-Indonesian Throughflow during the WOCE period 
(early 90s) is 15 ±5 Sv westward. The result is consistent with the global circulation, 
but I believe largely dictated by the use of the 1989 JADE cruise. The southward 
flux in the Mozambique Channel is of the same magnitude, implying a weak or null 
net flow between Madagascar and Australia. 

3. Upwelling and down-gradient dianeutral mixing are found in the northern part of the 
Indian Ocean below the surface layers, with insignificant dianeutral transfers south of 
about 8°S. The deduced “equivalent-horizontal average” diffusivity ranges between 
2 andl0cm 2 s _1 while the upwelling velocity ranges between 1 and 3 x 10~ 5 cms -1 . 
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3.2 The Pacific Ocean 


3.2.1 Introduction 

While the Indian Ocean discussion was relatively extensive, I choose to keep the Pacific 
Ocean circulation description brief, the main focus being on the meridional overturning 
circulation (MOC). A large literature describes the Pacific circulation (e.g., Talley and 
Roemmich, 1991), and other studies are based on the same sections as the ones used here 
in the North Pacific (Roemmich and McCallister, 1989, hereafter RMC89; Macdonald, 
1998) and in the South Pacific (Tsimplis et al., 1998; Wijffels et al., 1999, Wunsch et ah, 
1983). 

Deep water inflow 

In the South Pacific, regional estimates of the deep water inflow are relatively consistent. 
Wunsch et ah (1983) estimated net northward deep water flow of about 12 Sv using 
the Scorpio sections (Stommel et ah, 1973; Warren, 1973) at 28°S and 43°S. At 17°S, a 
similar northward bottom flow of 11 to 14 Sv below 9 = 1.1°C was found from current 
meter and hydrographic measurements (Roemmich et ah,1996; Banks et ah, 1995). Recent 
measurements in the western Southwest Pacific basin at 32°S evidenced a net northward 
flow of 15 Sv across a current meter array between the Tonga-Kermadec Ridge and 169°W 
and below 3250dbars, with a ±10Sv variability (Whitworth et ah, 1999). Tsimplis et 
ah (1998), using the same WOCE sections that I use in the South Pacific, obtained a similar 
inflow of about 12 Sv at both latitudes. However, the global model of Macdonald (1998) 
yielded a moderate bottom water inflow of 8±5 Sv at 43°S from the older Scorpio data. 

At 10°N, Wijffels et ah (1996) estimated a deep inflow of 8Sv. In the North Pacific 
north of 24°N, RMC89 estimated the circulation using an inverse model combing zonal 
and meridional sections and requiring top-to-bottom mass conservation. They obtained 
a net transport of bottom water of 10 Sv to the north across 24°N and of only 0.7 Sv to 
the north across 47° N. In this same region, the Macdonald (1998) model yielded a weaker 
bottom water flux of 5±3 Sv at 24°N, although it could be consistently constrained to a 
10-Sv flow. She found no significant flow at 47°N. As will be seen, a significant northward 
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bottom water flow at 24°N is incompatible with conservation of heat and salt in the North 
Pacific. 

3.2.2 Model constraints 
Data and treatment 

Figure 3-5 shows the Pacific Ocean sector of the model which is made of the pre-WOCE 
sections PI (Aug. 1985) and P3 (Apr.-May 1985) and the WOCE sections P21 (Apr.- 
Junl994) and P6 (May-Jul. 1992). I also attempted to use P4 (10°N, Feb.-May 1989) but 
it produced large downward mass transfers at intermediate levels, probably because of 
the baroclinic variability in the Mindanao Current (S. Wijffels, personal communication, 
1998) which made the section incompatible with the rest of the circulation. It was thus 
decided not to use P4 in the present model. 

The hydrographic data were provided by the Pis L. Talley (PI), D.Roemmich and 
J. Swift (P3), M. McCartney (P6, P21, P4), H.Bryden (P6, P21, P4), Brady (P4) and 
J. Toole (P4, P6) either directly or through the WHPO center. I chose not to use P2 
(Bando, Fukasawa and Okuda; Nov.-Feb 1994) because of the large time offset from the 
surrounding sections P3 and PI. 

A gap over a few stations in the phosphate data set in P6 was filled by using a local 
nitrate-phosphate Redfield relation. As described above, relative geostrophic velocity was 
computed using horizontal extrapolation in the bottom triangles. In trenches where no 
horizontal extrapolation was possible, a polynomial was fit to the temperature. Any 
unrealistic-appearing shear in the triangles was set to zero shear. 

Layers and constraints 

The Pacific boxes defined by these sections are respectively the South subtropical box (P6- 
P21); the central Pacific box (P21-P3-JADE89); the North Pacific box (P3-P1) and the 
Aleutian box (north of PI). Conservation equations were written for deep-water layers in 
the Aleutian box because water properties suggest that there is no deep water formation 
in the North Pacific (e.g., Warren, 1983), the deepest winter convection affecting layers 
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Figure 3-5: Pacific Ocean sections and topography. The 4000 m isobath is contoured and 
areas shallower than 3000 m are shaded. Box names are indicated in the shaded rectangles. 
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Layer 

Upper 

interface 

24°N Properties 

Water mass 

name 

1 

surface 

a 9 = 24.0; 9 = 22.7 

Surface Water 

2 

7 n - 24.3 

a e = 24.3; 9 = 21.9 

3 

7 n = 26 

ag = 26.0; 9 = 12.6 

North Pac. Int. Water 

4 

7 n = 26.5 

ag = 26.4; 9 = 9.1 

5 

7 n = 26.8 

ag = 26.7; 9 = 7.1 

6 

7 " = 27 

ag = 26.9; 9 = 5.9 

7 

7 71 = 27.15 

ag = 27.0; 9 = 5.2 

8 

7 n = 27.35 

ag = 27.2; 9 = 4.4 

Antarctic Int. Water 

9 

7 ™ = 27.6 

<r 2 = 36.56; 9 = 3.5 

10 

7 " = 27.8 

er 2 = 36.77; 9 = 2.5 

11 

7 71 = 27.9 

<7 2 = 36.88; 9 = 2.0 

12 

7 n = 27.95 

<j 2 = 36.92; 9 = 1.75 

North Pac. Deep Water 

13 

7 n = 28 

<r 2 = 36.97; 9 = 1.5 

14 

7 " = 28.05 

<7 2 = 37.01; 9 = 1.3 

15 

7 n = 28.09 

9 = 1.15 

16 

7 n = 28.1 

0 = 1.10 

17 

7 ” =28.11 

<74 = 45.885; 9 = 1.06 

Antarctic Bot. Water 

18 

7 n = 28.125 

<74 = 45.900; 9 — 1.02 

19 

7 n = 28.15 

9 = 0.85 

20 

7 n = 28.18 

9 = 0.85 


Table 3.5: Pacific Ocean neutral surface layers T (in kgm~ 3 ). Interfaces 15 and 17 have 
been used in the North Pacific only. The Macdonald (1998) water mass definition are 
indicated (from her table 17). 

above lOOOdbars at 47°N (e.g. Reid, 1965, pp. 50-51; Reid, 1973; Talley et al., 1998, Fig. 
10). Heat and salt are thus conserved in layers north of 47°N below lOOOdbars, allowing 
for dianeutral advective and diffusive fluxes. (With no constraints north of 47°N, the 
model had an unphysical albeit non-significant southward net flow of —2 ± 2 Sv below 
3000dbars. Inclusion of the Aleutian box reduced this flow.) 

The layers were suggested by John Toole (personal communication, 1998), optimized to 
match the Macdonald (1998) layers and to keep similar thicknesses over the whole Pacific 
(Table3.5). The initial Ekman transports from the NCEP re-analysis are given in the 
second column of Table 3.6. 

As in other oceans, the circulation is required to conserve mass, salt, silica (top-to- 
bottom); and “PO” and heat away from the surface. The net northward flux across P 6 and 
P21 was constrained to the Bering strait (0.8 Sv, Coachman and Aagaard, 1988) +PIT 
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Section 

Initial 

Final 


P6 

QE3 

1.3±1 




-20±6 

15% 

P3 

7.5±3 

9±3 

13% 

PI 

-4±2 

-3.9T1.6 

-3% 


Table 3.6: Pacific initial Ekman transport and correction 

flux (7Sv, see the discussion in the Indian Ocean section), resulting in a net constraint of 
7.8±10 Sv [10=\/7 2 + 7 2 , where the first ‘7’ is the a priori uncertainty in the throughflow 
and the second ‘7’ is the variability in the baroclinic field according to the error budget]. 
The total flux across PI and P3 was constrained to 0.8±3Sv and 0.8±6Sv respectively. 

The deep water flow across P6 between 179°W and 168°W was constrained with the 
current meter measurements of Whitworth et al. (1999) with 16±5Sv northward below 
7 n = 27.95; 6 = 1.75°. (The uncertainty of ±5Sv corresponds to their estimated vari¬ 
ability in the baroclinic transport from repeat sections, which is half the total variability 
as measured from the current meters. Importantly, the authors found the P6 baroclinic 
transport to be representative of the average of 3 repeat sections.) Additional northward 
constraints were applied for each depth range, from the Whitworth et al. (1999) measure¬ 
ments and following Wijffels et al. (1999): 2000-3250 dbars: 1 ± 3 Sv; 3250-4250: 9 ± 3 Sv; 
4250-bottom: 6 ± 2.5 Sv. (Because of my layer definition, the depth range separation is 
taken at 4250 dbars instead of the original 4750 dbars for which Whitworth et al. (1999) 
provided the transports (their Fig. 17). The uncertainties are assumed to take into account 
possible shifts in the transport.) Uncertainties correspond to half the total variability men¬ 
tioned by Whitworth et al. (1999). The deep flow through P21 was similarly constrained 
with 12±4Sv below = 27.95, according to the analysis of Roemmich et al. (1996) in 
the Samoa Passage region. (Note that while Whitworth et al. (1999) estimated only the 
transport on the Western part of the deep water, Roemmich et al. (1996) estimated the 
total net northward transport of water colder than 9 = 1.1°C.) 

The deep flow was also constrained according to observed bathymetric features (Ta¬ 
ble 3.7). A relatively large tolerance of ±5Sv was given to the South West Basin, Yupan- 
qui Basin (20°S-105°W), and Peru-Chile Basin to allow dianeutral exchanges with layers 
above. 
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Name 

Section 

Layer depth 

Net flux 

Longitude 

N. Tasmanian 



0 ± 2 Sv 

153°E-156°E 

S. NewCal Trough 

P 6 

7 n = 27.9 (2000) 

0 ± 2 Sv 

162°E-168°E 

Yupanqui Basin 

P21 

7 n = 28.05 (3700) 

0± 5 Sv 

94°W-78°W 

Bauer Basin 

P21 

7 n = 28.05 (3900) 

0 ± 2 Sv 

107°W-100°W 

N. Peru Basin 

P21 

7 n = 28.05 (3900) 

0 ± 2 Sv 

94°W-75°W 

Okinawa Trough 

P3 

7 n = 27.6 (1200) 

0 ± 2 Sv 

125°E-127°E 

N. Philippines Bas. 

P3 

7 n = 28.05 (3100) 

0± 2 Sv 

127°E-141°E 

Kuril Trench 

PI 


0 ± 2 Sv 

145°E-147°E 

Name 

Box 

Layer depth 

Net Div. 

Longitudes 

Tasmanian 

P6/P21W 

7 n = 28.1 (3900) 

0± 2 Sv 

153°E-157 0 E/153 0 E-156°E 

NewCal Trough 

P6/P21W 

7 n = 27.9 (2000) 

0± 2 Sv 

162°E-168 0 E/162 0 E-164°E 

S. Fiji Basin 

P6/P21 

7 n = 28.0 (2700) 

0 ± 2 Sv 

173 0 E-179°E/165 0 E-168°E 

South West Basin 

P6/P21 

7 n = 28.05 (3600) 

0±5Sv 

179°W-122 0 W/173 0 W-114°W 

Peru-Chile Basin 

P6/P21 

7 71 = 28.05 (3900) 

0± 5 Sv 

112 o W-71 0 W/77 o W-75°W 


Table 3.7: Constraints imposed on the deep layer flux by bathymetric features. The 
constraints are imposed on either the net flux through a single section for a longitude 
range (upper part of the table, positive northward) or on the net divergence (Div.) in a 
box between two sections (lower part of the table). The layer and its approximate depth 
below which the flux is constrained are indicated, along with the uncertainty on the net 
flux (section) or divergence (box) (in Sv). 


Reference surfaces and dianeutral terms 

The reference surfaces, indicated in Table 3.8, were chosen roughly at the top of the 
northward-flowing Circumpolar Deep Water. At 32° S, near the deep western boundary 
current against the slope of the Kermadec ridge, the current meter results of Whitworth 
et al. (1999) showed that the limit of the northward-flowing waters did not coincide with 
any “water-mass boundary” as classically used as reference surfaces. The initial velocity 
field was thus adjusted to match their 2-year mean absolute velocity. At 17°S against the 
Tonga Ridge, the reference surface was moved just above the high shear region of the deep 
western boundary current. Away from the boundary currents, I chose 9 = 1 . 1 ° or 1.3° as 
a reference (i.e., Banks et al., 1995; Tsimplis et al., 1998). However, an experiment with 
reference surface on 7 " = 27.95 (9 ~ 2°) showed that integrated mass transports after 
inversion are relatively insensitive to reference surface choice in this area. Along sections 
P3 and PI, the reference surfaces were chosen based on RMC89. Dianeutral terms and 
reference velocity adjustment were chosen as for the Indian Ocean, page 80. 
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Section 

Reference surface 

Remark 

P6 

7 n = 28.10; 0-1.1° 

match to Whitworth et al. (1999) in DWBC 

P21 

7 n = 28.10; 0- 1.3° 

Banks et al. (1995) 

P3 

7 n = 27.95; 0 - 1.6° 

2000dbars Roemmich and McCallister (1989) 

PI 

7 n =28.09; 0-1.1° 

4000 dbars-initial velocity similar to RMC89 


Table 3.8: Reference surfaces in the Pacific Ocean. 


a) South Subtropical 32S-17S 
o 


b) Central Pacific 17S-24N c) North Pacific (24N-47N) 


1000 

2000 

3000 

4000 

5000 

6000 


-5 


-5 



-0.15±4.5x 10 9 kg/s 


-2.1±4.9x 10 9 kg/s 


1.3±3.2x 10 9 kg/s 


Figure 3-6: Standard solution mass residuals in a) the South Pacific box, 30°S-17°S; b) the 
Central Pacific box, 17°S-24°N and c) the North Pacific box, 24°N-47°N. The line gives 
the residual in invidual layers (in 10 9 kgs -1 ) and the shaded area give the one standard 
deviation uncertainty. The top-to-bottom residual is indicated below the graph. 


3.2.3 Pacific Ocean circulation 

Almost all residuals of the equations of the “standard” solution were within the one stan¬ 
dard deviation a priori range. The Whitworth et al. (1999) constraint of net flow of 16±5 Sv 
across their current meter array was met with 11 ± 3 Sv; however, constraints on the ver¬ 
tical partition of the flow were not met within one standard deviation, but within two 
standard deviations of the a priori: the model produced northward transport higher than 
the ones prescribed below 4250dbars, 10 ± 2Sv, and lower just above, 4 ± 1 Sv.) Mass 
residuals in individual layers are smaller than 0.5 Sv and all indistinguishable from zero 
within two standard deviations (Fig. 3-6). Small adjustments were made to the initial 
Ekman transport (3 to 15% , Table 3.6). 

Again, the discussion that follows is focused on the large scale integrated fluxes. How¬ 
ever, reference velocities, absolute velocities and water-mass transports accumulated from 
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a) P6 (32S) 



16.1±5.1x 10 9 kg/s 
d) P3 (24N) 



—1.21 ± 3.3x 10 9 kg/s 


b) P21W+P21 (17S) 



15.6± 5.1 xIO 9 kg/s 


e) PI (47N) 



0.221 ± 1.5x 10 9 kg/s 


c) P21W only 



-8.52± 11 x 10 9 kg/s 


Figure 3-7: Pacific Ocean standard circulation. The mass transport across each section 
is shown, integrated from the bottom (the overturning streamfunction ). The shaded area 
corresponds to the uncertainty (one standard deviation). The net transport is indicated 
below each graph. While panel b) gives the total transport across 17°S, panel c) gives the 
transport between New Caledonia and Australia only (P21W). 


west to east are given in Appendix D (pages 210-219) for reference. The integrated flow 
for the standard solution for the Pacific Ocean circulation is shown on Figure 3-7. (Again, 
overturning streamfunctions are used, the slope of the curve indicating the direction of 
the flow). 


Net fluxes 

The net fluxes across 32°S and 17°S corresponding to the Pacific-Indonesian Throughflow 
(PIT) plus Bering Strait transport are estimated to be 16±5Sv going north. The total 
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flow between New Caledonia and Australia is uncertain and to the south, -8±11 Sv (Fig. 3- 
7c), implying a 24±12 Sv northward flux between New Caledonia and Chile. Tsimplis et 
al. (1998) obtained similar transports. At 32°S in the Tasman Sea, however, the total 
transport between Australia and 170°E is also southward, with about —10 Sv, in apparent 
contradiction with Wijffels et al. (1999) who found a closed circulation in the Tasman Sea 
implying no net transport. 

Deep inflow 

In the South Pacific, net northward flow of deep water is found at approximately 8Sv 
(Fig. 3-7a, b). At 32°S, most of the flow occurs in the bottom layer, below 7™ = 28.18 
(1 9 ~ 0.8°) at 8 ± 2Sv, while at 17°S the bulk of the northward flow occurs higher in the 
water column, with 7± 1.3 Sv below 7” = 28.1 (9 ~ 1.2°). Dianeutral transfers (described 
below) account for this shift. The deep western boundary current against the Tonga- 
Kermadec ridge at 32°S has a strength of 17±3Sv below 7” = 27.95 (9 ~ 1.9°) (page 211). 
(The DWBC was actually flowing south at the time of section P6, but probably as a result 
of a strong barotropic fluctuation, thus not visible in the hydrographic data.) The total 
flow between the Tonga-Kermadec Trench and 168°W is weaker but still consistent within 
error bars with the Whitworth et al. (1998) current meter results with 11 ± 3Sv between 
179°W and 168°W below 7 n = 27.95 (9 ~ 1.9°). 

In the North Pacific at 24°N (Fig. 3-7d) the deep inflow is indistinguishable from zero 
in the bottom layers below 7” = 28.09 (9 ~ 1.15°), thus much weaker than the RMC89 
flow of 10 Sv. Still, there is a total eddy 1 transport of northward-flowing water of 7±6Sv 
below 9 = 1.15°, but this same amount returns south in the same isopycnal depth range. 

The deep inflow into the Pacific Ocean estimated with this model is thus lower than 
most of the published estimates mentioned above that found 10 to 13 Sv from regional 
analysis (e.g., Wijffels et al., 1999; Whitworth et al. 1999; Tsimplis et al., 1998; Wunsch 
et al., 1983) and 5 to 8 Sv from the global Macdonald model. The RMC89 circulation 
conserved mass only for the top-to-bottom integrated transport within boxes enclosed 

'Northward eddy transport is defined here as the sum, over all station pairs, of the flow that is oriented 
northward within the considered depth range. 
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by zonal and meridional sections: no conservation statements were applied in individual 
layers. (They used the same 24°N and 47°N sections as I do but added a zonal section 
at 36°N and several meridional sections.) Their estimated 10 Sv inflow implies a large 
dianeutral velocity of 3xl0“ 4 cms _1 (RMC89, page 201) and an associated diapycnal dif- 
fusivity of around 60 cm 2 s _ 1 . Such a diffusivity is in the upper range of local measurements 
in the presence of rough topography (e.g., Polzin et al., 1997) whereas the North Pacific 
topography is relatively smooth. 

The net deep inflow in the standard solution is limited mainly by the conservation 
equations in the North Pacific deep layers, and in particular by the salinity anomaly 
conservation equations. An experiment with heat conservation only, produced a 5 Sv deep 
inflow with strong upwelling (w* ~ 4 x 10~ 5 cms _1 ) and vertical mixing (60cm 2 s _1 ) at 
4000 dbars between 24°N and 47°N. An experiment with a deep reference surface (-y n = 
28.11, which is close to the upper limit of the deep inflow found by RMC89) produced a 
qualitatively similar circulation. It was tried to constrain the deep infow to 10 ± 5 Sv to 
the north, but the model could not meet the constraint and the circulation and property 
fluxes remained almost unchanged. The RMC89 circulation in the North Pacific could be 
reproduced only by suppressing conservation in individual layers. The Macdonald (1998) 
circulation, which produced a more moderate inflow of 5±3Sv at 24°N, had no heat 
conservation constraints, and the fact that anomaly equations were not used gave the salt 
conservation equation a much lower weight in her result. Salt and heat conservation thus 
greatly limit the net inflow of bottom water into the North Pacific basin north of 24°N, 
and in turn also limit the inflow in the South Pacific. As with the case of the Indian 
Ocean meridional overturning, it is not known whether this low inflow is an actual feature 
or whether it is model dependent, for instance due to poor representation of vertical 
exchanges. 

Because no net deep flow enter the North Pacific, the deep inflow in the South Pacific 
Ocean is all upwelled between 17°S and 24°N. The estimate of the South Pacific deep inflow 
is related with the weight of the PO constraint. An experiment with no PO constraints 
yielded a 14.8±3 Sv deep inflow associated with larger mass residuals at intermediate and 
deep levels. This solution is physically acceptable as well. However, the two solutions 
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are in apparent contradiction (the uncertainties on the deep inflow of the two solutions 
hardly overlap within two standard deviations at 10 Sv). This contradiction suggests that 
the a priori weights on some of the constraints may be too high. In particular, there is a 
trade-off between the constraints on the deep inflow derived from the current meter data 
of Whitworth et al. (1999), the PO constraint, and mass residuals in the Central Pacific. 
The current meter data may be hampered by the shortness of the time series (22 months), 
while the weight on the PO constraints was chosen analogously to the salinity, according 
to the “ad-hoc” method described on Section 1.2.4 (page 22). I chose the standard solution 
as the one of low-overturning because of the smaller mass residuals, but this choice is, to 
some extent, arbitrary. 

Deep, intermediate and surface flow 

The horizontal flow structure in the density class defining the different water masses de¬ 
fined in Table 3.5 are given in Appendix D (Figures D-16-D-24). The deep water flow 
in the South Pacific is directed southward as found by Wunsch et al. (1983), Macdon¬ 
ald (1998) and Tsimplis et al. (1998). Most of this flow occurs between 170°W and 140°W 
at 32°S (Fig.D-16). Such a deep flow contrasts with the speculative interior circulation 
that Whitworth et al. (1999) inferred based on a Stommel and Arons (1960) schematic 
and current meter measurements near the Kermadec Trench. But the present circulation 
is consistent with tracer conservation, including silica conservation over the Pacific Basin 
(Section 4.3.1). 

In the South Pacific at 30°S (P 6 , Fig.D-16), the East Australian Current closes the 
subtropical gyre with a strength of -36±10 Sv southward in the region between Australia 
and 154.5°E and from the surface to 7 " = 27.95, similar to previous estimates (e.g., Mac¬ 
donald, 1998, page 348; Chiswell et al., 1997). (The formal uncertainty of the boundary 
current transports are underestimated because the error budget did not take into account 
baroclinic variability on small scales.) At 25°S, the estimated East Australia Current 
strength is somewhat smaller at -25±6 Sv (from Australia to 155°E). Both estimates are 
also consistent with those of Tsimplis et al. (1998). At both latitudes, an immediate re¬ 
circulation returns about 80% of the flow to the north within 500 km of the current. The 
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Peru-Chile Current is apparent on the eastern side of P 6 (Fig. D-16), carrying 8.5±4 Sv to 
the north at 30°S (surface to j n = 27.95; 75°W to Peru). At 16°S (P21E), the flow is to 
the south against the coast (-8±2 Sv), with a 22±2 Sv northward jet offshore (Fig. D-20) 
and eddy-like motion farther west. 

In the North Pacific (Figs. D-22-D-24), the intermediate and surface flow field is very 
similar to that of Macdonald (1998) at both sections, the differences arising mainly in 
the deep layers. The subtropical gyre at 24°N is closed by the Kuroshio (68±7Sv to the 
north between the surface and y n = 27.95 and from the western boundary to 130°E). A 
similar flow of 60 Sv was reported from earlier data (Worthington and Kawai, 1972, page 
373). About 2/3 of this flow is returned south between 132°E and 134°E and 1/3 over 
the rest of the basin. The Kuroshio transport in the East China Sea alone, 28 ± 6 Sv, is 
consistent with estimates from Macdonald (1998) and RMC89. However, recent current 
meter observations from the WOCE PCM -1 array 400 km upstream suggest a somewhat 
smaller flux of 22±2.5 Sv across the strait between Taiwan and the Ryukyu Islands (Lee et 
al., 1999; Johns et al., 1999). It is not known whether the difference is due to interannual 
variability (e.g., Ichikawa and Beardsley, 1993), or to additional mass influx through the 
Ryukyu Islands. At 47°N, a northward flow is found close to the coast (31±7Sv from 
the surface and 7 ” = 27.95 and west of 146.5°E), adjacent to the Oyashio that carries 
-37±2Sv (146.5°E to 148.4°E) south. On average, southward flow is expected against the 
coast just north of this region, but the PI section reaches the coast close to the Kuroshio 
Front, which generates warm eddies (Kawai, 1972, page 248). The northward flow that is 
found against the coast is probably associated with such an eddy (i.e., Talley et al., 1991). 

Dianeutral transfers and diffusivity 

Dianeutral mass transfers indicate upwelling at the surface in the Central Pacific (Fig. 3 - 
8 b, 45±6 Sv) and downwelling in the subtropical gyres (-16±5 Sv in the South Pacific; 
-17±3 Sv in the North Pacific) (Fig.3-8a,c). At depth, and in the South and Central Pacific, 
there is upwelling below 2000 dbars associated with systematic and significant diffusivities 
ranging from 1 to 8 ± 2 cmV 1 between 2000 and 4500 dbars (Fig. 3 - 8 a,b). Diffusivities 
across the deepest layer interface are larger at about 40cm 2 s _1 . The large diffusivities 
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a) South Subtropical 32S-17S 



b) Central Pacific 17S-24N c) North Pacific (24N-47N) 




x10 9 kg/s x10 9 kg/s 





Figure 3-8: Dianeutral transfers (upper) and diffusivities (lower) in the Pacific Ocean. 
The line indicates the transfer or diffusivity at the bottom of the layer. 
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are related to small horizontal interface areas. But the area calculation is uncertain for 
the deep layers as it was derived using Levitus (1994). A back-of-the-envelope heat bal¬ 
ance by RMC89 and Roemmich et al. (1996) suggested possibly dianeutral diffusivities 
of 60 to 1000cm 2 s _1 in the deep layers downstream from the Samoan Passage. In the 
deep North Pacific, dianeutral diffusivities are positive, but not significant. North of 47°N 
(not shown), dianeutral transfers are not distinguishable from zero. Average values for 
dianeutral diffusivities and velocities are given at the end of this chapter. 

The magnitudes of the vertical transfers and diffusivities are within the range of pre¬ 
viously published values (e.g., Tsimplis et al., 1998; Macdonald, 1998; Wijffels, 1993). 
Tsimplis et al. (1998) obtained a smaller upwelling of a few tenths of a Sverdrup below 
3500 dbars and downwelling above 3500 dbars. Macdonald (1998) obtained no significant 
vertical exchanges there. On the other hand, she estimated significant upwelling in the 
North Pacific which I do not because of the weaker deep inflow discussed above. 

3.2.4 Pacific circulation summary 

The standard estimate of the Pacific circulation is associated with an inflow of CDW 
that is weaker than most previous estimates. The small size of the inflow is related to 
conservation of heat and salt in the deep North Pacific layers, which precludes significant 
net inflow across 24°N. As a result, the inflow from the South is limited also, and most of 
the circumpolar deep water (CDW) moving northward across 32°N upwells in the Central 
Pacific. By removing the PO conservation constraint, a solution with higher CDW inflow in 
the South Pacific was obtained, but the associated mass residuals increased in the Central 
Pacific. I chose, arbitrarily to some extent, to retain the PO-conserving circulation as my 
standard solution. The bulk of the bottom water inflow returns at deep levels so that the 
deep overturning cell is isolated, in the integrated sense, from the intermediate and surface 
water circulation. The South Pacific thermocline flow is directed northward, suggesting 
that the density class found in the Pacific-Indonesian Throughflow ultimately originates 
in the South Pacific Ocean south of 32°S. (However, the pathway is complicated, e.g., Fine 
et al. (1994)). Such a circulation was also inferred by previous authors (e.g., Macdonald, 
1998; Wijffels et al., 1996; RMC89; Wunsch et al., 1983). 
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Dianeutral diffusivities in the deep layers indicate significant and consistent values in 
the range 2-8cm 2 s _1 , with higher values in the bottom layers. The model transports of 
surface boundary currents are consistent with previous estimates. 
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3.3 The Atlantic Ocean 


3.3.1 Introduction 

The North Atlantic Ocean is the best documented in comparison with the other oceans. 
The meridional overturning circulation consists of 13 to 30 Sv northward flow of thermo- 
cline water coming from the South Atlantic that eventually sinks at high latitudes to form 
the southward flowing North Atlantic Deep Water (NADW). The NADW flows to the 
south, mainly in the Deep Western Boundary Current (DWBC) whose transport increases 
owing to the contribution of Atlantic Antarctic Intermediate Water (AAIW) from above 
and Antarctic Bottom Water (AABW) from below. The net NADW export to the South¬ 
ern Ocean is approximately 20 Sv according to most authors. AABW is formed in the 
Weddell Sea and around the continental shelf of Antarctica and flows north in the At¬ 
lantic at a rate of about 8Sv at 30°S, according to current meter measurements (Holfort 
et al., 1999b, 8Sv; Hogg et al., 1999, 7 Sv). Part of this flow is entrained into south¬ 
ward flowing NADW within the Brazil Basin; about 2 Sv enters the Guiana Basin (Hall 
et al., 1997) and 1.2 Sv enters the Eastern Atlantic through the Romanche Fracture Zone 
(Mercier and Speer, 1998). 

Elements of this large-scale circulation has been estimated quantitatively in several 
studies, including the inverse box model approaches of Wunsch (1978, western North At¬ 
lantic); Wunsch and Grant (1982, whole North Atlantic); Hall and Bryden(1982, 24°N); 
Roemmich and Wunsch (1985) and Rintoul and Wunsch (1991) (subtropical North At¬ 
lantic). In the South Atlantic, the circulation was estimated by Fu(1986); Lux (1997, 
equatorial region); Holfort et al. (1999); and Saunders and King (1995, 45°S). But the 
most comprehensive circulation picture probably comes from the global model of Mac¬ 
donald (1998). Besides box section models, a variety of other inverse models has been 
used, e.g., Martel and Wunsch (1993, linear, finite-difference, steady); Mercier et al. (1993, 
non-linear) and Schlitzer (1988, 1993, coarse resolution, intermediate complexity model), 
each bearing its own skills and drawbacks. 

In addition to those estimates, recent qualitative pictures of the North Atlantic circu¬ 
lation and its connection with other basins have emerged from water mass analysis and 
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current meter measurements by Reid (1994, 1996) and a literature review by Schmitz and 
McCartney (1993) and Schmitz (1996a). (The two later reviews include quantitative esti¬ 
mates from many authors, but the resulting estimate of the circulation is presented with 
no error bars.) 

The different water mass transports presented in these papers are discussed below in 
a quantitative comparison with my results. 

3.3.2 Model constraints 

Data and treatment 

Figure 3-9 shows the Atlantic Ocean sector of the model which is made of the WOCE sec¬ 
tions A2 (Gauss 226, Jul. 1993, Koltermann et ah, 1999); A5 (Jul. 1992, Parillaet al., 1994); 
A6 (Feb. 1993) and A7 (Jan 1993) (Arhan et al., 1998; Oudot et al., 1998); A8 (Mar. 1994), 
A9 (Feb. 1991) and A10 (Dec. 1993) (Siedler et al., 1996); and All (Dec. 1993)(Saunders 
and King, 1995). I was not able to incorporate the available non-WOCE-standard data 
at 36°N (see page 14). At 26°N in the Florida Strait (a continuation of section A5), the 
1992 data did not contain nutrients, so that CTD and nutrient data from the 1981 AT109 
section of Roemmich and Wunsch (1985) were used. The low latitude and small station 
spacing enhanced the noise in geostrophic velocities across the A6 and A7 sections, in 
particular against the coast of Brazil. To decrease the noise, a horizontal Gaussian filter 
with cutoff at 70 km was applied to the temperature and salinities in the upper 1000 dbars. 
Because the A6 section does not include the part of the North Brazil Current that flows 
on the continental shelf, a 5-Sv transport was added to the section (see below). 

The hydrographic data were provided by the Pis P. Koltermann (A2); G.Parillaand 
R. Millard (A5); M. Arhan, A.Morliere, C. Oudot, C. Colin and Y. Gouriou (A6, A7); 
G. Siedler, T. Mueller and R. Onken (A8, A9, A10); and P. Saunders (All), either directly 
or through the WHPO Center when public. 
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Figure 3-9: Atlantic Ocean sections and topography. The 4000 m isobath is contoured and 
areas shallower than 3000 m are shaded. The light rectangles indicate the box labels. 
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Layers and constraints 

The Atlantic boxes (Figure 3-9) are respectively the North Atlantic box (A2-Florida Strait- 
A5); the North Tropical box (A5-A6); the Equatorial box (A6-A7); the South Tropical 
boxes “Ascension” (A7-A8) and “Helena” (A8-A9); the South subtropical “Rio” box (A9- 
A10) and the South subtropical “Tristan da Cunha” box (A10-A11). In addition, a super 
“Central Atlantic” box was defined between A5 and A9 to impose tighter conservation 
constraints between those two sections. Without this super box, coupling between the 
flow north and south of the Equator is loose because of large noise in the equatorial region 
transport calculations. 

In the North Atlantic, the neutral surface layers were chosen to match the Macdon¬ 
ald (1998) layers and water mass definition (Table 3.9) for comparison with her results and 
other authors (e.g., Rintoul and Wunsch, 1991). The Macdonald nominal layers at 24°N 
(<?j) are given in the second column for reference. (Note that Macdonald changed the 
reference pressure when the potential density surface departed far from it. Her resulting 
layers consistently matched neutral surfaces, in contrast with the original a layers which 
occasionally led to large departure from neutral surfaces, e.g., a 500dbar deviation was 
found near the Mediterranean waters.) 

In the South Atlantic, the layer and water mass definitions were based on both Mac¬ 
donald and Speer et al. (1996). The nominal interface between NADW and AABW, 
7 n = 28.11, remains close to 0 = 2° between 30°S and 7°N, but varies from 1.4° to 
1.8° at All (45°S-30°S). 

The initial Ekman transports derived from the NCEP re-analysis are given on the 
second column of Table 3.11. 

Conservation equations are applied as described in Chapter 1, and weighted according 
to TableC-2 (page 194). The “PO” conservation constraints were given lower weights in 
boxes “Rio”, “Tristan da Cunha” and superbox A5-A9 because with the initial weight from 
Section 1.2.4 (page 22), the derived circulation contained mass residuals in deep layers and 
systematic negative diffusivities in several boxes. The net flux across all latitudes was con¬ 
strained to the Bering Strait volume flux and salt fluxes. Coachman and Aagaard (1988) 
estimated the Bering Strait flow to -0.8±0.6Sv for the mass and 26.7xl0 6 kgs -1 for the 
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salt. However, the uncertainties applied on the net flux constraints across individual sec¬ 
tions are much larger, according to the error budget (TableC-l, page 193). The Florida 
Strait transport was constrained to 31±lSv, following the measurements of Schott et 
al. (1988a) 2 . Section A6 does not include the part of the North Brazil Current (NBC) that 
is on the continental shelf. Current meters deployed near 7°N produced an estimate of 
this transport of 4.8 Sv (Johns et al., 1998). Similarly to Lux (1997), the 5-Sv continental 
NBC was included in my tracer budgets by using the average tracer concentration in the 
first lOOdbars of the westernmost station. Because the shelf water is warm (26°C), the 
influence of the NBC transport is important (0.5 PW). 

In the south, the AABW flux was constrained to 4±0.4Sv northward in the Vema 
Channel (Hogg et al., 1982) and to 2.9±lSv in the Hunter Channel (Zenk et al., 1999). 
The total flux was constrained to 7±2 Sv. Further North, the AABW flux was constrained 
to 5.5±2Sv in the western part of the basin at 11°S (A8), according to McCartney and 
Curry (1993); to 3±6Sv at 4.5°S (A7) and 2±4 Sv at 7.5°N (A6), in the western basin 
according to the current meter measurements made in the same basin but at the Equa¬ 
tor by Hall et al. (1997). (There too, the weights were defined by the error budget.) 
The AABW flux through the Romanche and Chain Fracture zones was constrained to 
1.4±3.6Sv (Mercier and Speer, 1998). 

The Deep Western Boundary Current was weakly constrained at 7.5°N to 20±10Sv, 
according the SOFAR float measurements of Richardson and Schmitz (1993) (between 
900 and 2800dbars and within 700 km of the coast). At 19°S, the DWBC transport was 
constrained according to the current meter measurements of Weatherly et al. (1999), that 
is, -15 ± 7Sv (1200-2000 dbars); -11 ± 4Sv (2000-2800 dbars); and -13 ± 6Sv (2800- 
3600dbars). The “bathymetric constraints”, essentially due to the Mid-Atlantic Ridge, 
are listed on Table 3.12. 

2 An uncertainty of ±lSv is lower than the one suggested in the preceeding error budget. But I kept 
this original uncertainty from independent measurements because, from the many estimates derived from 
hydrographic data in the Florida Strait, all seem to be consistent with the direct measurement to this 
accuracy. 
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Layer 

Upper 

interface 

24°N (A5) Properties 
ax from Macdonald (1998) 

Water mass 
label 

1 


9 = 17.9 

Surface Water 

2 

BBiM 

og = 26.40; 9 = 15.8 

Thermocline 

Water 

3 


a e = 26.80; 9 = 13.6 

4 

miw?smwm 

<jg = 27.10; 9 = 11.3 

5 


a e = 27.30; 9 = 8.80 

Intermediate 

Water 

6 

7 n = 27.62 

ag = 27.50; 9 = 6.90 

7 


a e = 27.70; 9 = 5.20 

NADW 

8 


<r 2 = 36.87; 9 = 3.60 

9 


<r 2 = 36.94; 9 = 3.20 

10 


02 = 36.98; 9 = 3.00 

11 

■MggaftEEM 

<r 2 = 37.02; 9 = 2.76 

12 


a 4 = 45.81; 9 = 2.50 

13 


<r 4 = 45.85; 9 = 2.24 

14 

wiMmmjim 

<j 4 = 45.87; 9 = 2.11 

15 


<r 4 = 45.90; 9 = 1.93 

Bottom 

Water 

16 


<r 4 = 45.91; 9 = 1.85 

17 

EH 

<74 = 45.92; 9 — 1.81 


Table 3.9: North Atlantic Ocean neutral surface layers T (in kgm -3 ). The Macdon¬ 
ald (1998) water mass definition are indicated (from her Table 11). NADW stands for 
North Atlantic Deep Water. 



Upper 

interface 

11°S (A8) Temperature 

Water mass 
label 

1 

surface 

CM 

II 

Surface 

Water 

2 

7 n = 26.2 

0= 16.6 

3 

7 n = 26.8 

9 = 11.8 

4 

7 n = 27.3 

9 = 6.20 

AAIW 

5 

7 n = 27.5 

9 = 4.30 

Upper 

NADW 

6 

7 n = 27.72 

9 = 3.90 

7 

7 n = 27.903 

9 = 3.70 

8 


9 = 3.30 

Middle 

NADW 

9 


6 = 2.75 

10 

iWfrffFif 

9 = 2.40 

Lower NADW 

11 

7 n = 28.11 

9 = 2.00 

AABW 

12 

7 n = 28.20 

t-H 

II 

05 

13 

7 n = 28.23 

9 = 0.60 

14 

7 n = 28.25 

6 = 0.40 

15 

7 n = 28.27 

9 = 0.22 

16 

7 n = 28.296 

0 = 0.11 


Table 3.10: South Atlantic Ocean neutral surface layers T (in kgm 3 ). Speer et al. (1996) 
layer definition are used. Water mass definition are indicated (from their Table 4). 
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Section 

Initial 

Final 

% correction 

A2 

-2.7±1 

-2.7T1 


Florida Strait 

0 

0 

- 

A5 

+4.6T2 

3±1.7 

-32% 

A 6 

+8.2T4 

8.3T3.6 

1 % 

A7 

-17±8 

-19T6.8 

13% 

A 8 

-10±5 

-8.6±4 

-14% 

A9 

-5.54:3 

-6.3T2.6 

13% 

A10 

-0.5±1 

- 0 . 6±1 

- 

All 

+6±3 

5.8T2.8 

-3% 


Table 3.11: Atlantic initial Ekman transport and correction 


Name 

Section 

Layer depth 

Net flux 

Longitude 

Walvis R. North 

A10 


0 ±lSv 

7.3°W-1.7°E 

Walvis R. South 

A10 

7 " = 28.11 (4200) 

0 ±lSv 

2.2°E-13.4°E 

Name 

Box 

Layer depth 

Net Div. 

Longitudes 

East N. Atlantic 

A2/A5 

7 n = 28.0986 (3500) 

0±3Sv 

27 0 W-11°W/45 0 W-16°W 

East S. Atlantic 1 

A8/A9 

7 n = 28.11 (4000) 

0± 2 Sv 

13.2 o W-13.4°E/12.4 o W-10.8°E 

East S. Atlantic 2 

A9/A10 

7 n = 28.11 (4000) 

0 ± 2 Sv 

12.4°W-10.8 o E/13.6 o W-14.7°E 

East S. Atlantic 3 

A10/A11 

7 n = 28.20 (4500) 

0 ± 2 Sv 

13.6 0 W-14.7°E/15 o W-15.6°E 


Table 3.12: Constraints imposed on the deep layer flux by bathymetric features. The 
constraints are imposed on either the net flux through a single section for a longitude 
range (upper part, positive northward) or on the net divergence (Div.) in a box between 
two sections (lower part). The layer and depth below which the flux is constrained are 
indicated along with the uncertainty on the net flux (section) or divergence (box) (in Sv). 
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Section 

Reference surface 

Remark 

A2 

bottom (LCD) 

match to Clarke et al. (1998) in the DWBC 

FI. Strait 

bottom (LCD) 

25±20cms _1 initial (Schott et al., 1988a, Figure 2 ) 

A5 

7 n = 28.1295; 0 ~ 1.90° 

NADW/AABW 

A 6 

Pair-Specific 

Lux (1997) 

A7 

Pair-Specific 

A 8 

Pair-Specific 

Speer et al. (1996) 

A9 

Pair-Specific 

Weatherly et al. (1999) 

A10 

7 n = 28.11; 6 ~ 1.90° 

NADW/AABW 

All 

7 n = 28.11; 6 ~ 1.40° - 1.8° 

NADW/AABW - Saunders and King (1995) 


Table 3.13: Reference surfaces in the Atlantic Ocean. LCD stands for Last Common 
Depth of station pairs. MAR=Mid-Atlantic Ridge. 


Reference surfaces and dianeutral terms 

The reference surfaces, indicated in Table 3.13, were chosen approximately at the interface 
between AABW and NADW, or according to current meter measurements or previous 
studies. At 19°S, the reference surface was taken at the bottom in the first 300 km of the 
western end, based on the current meter measurements of Weatherly et al. (1999). Farther 
east, the reference was taken at 3000 dbars to avoid a high shear region below that depth 
west of the Mid-Atlantic Ridge. Still at 19°S but in the eastern Basin, the reference was 
taken at the bottom. Integrated fluxes after inversion were insensitive to this choice versus 
the 7 ” = 28.11 surface. 

3.3.3 Atlantic Ocean Circulation 

Except for the DWBC transport across 19°S between 2800 and 3000 dbars and the Ro- 
manche Fracture Zone transports, the standard model solution met all Atlantic constraints 
within one standard deviation of the a priori noise range. The model’s lower DWBC 
transport at 19°S was found to be weaker than the Weatherly et al. (1999) constraint 
(3Sv southward), but still within 2 standard deviations. Mass residuals are indicated on 
Fig. 3-10. These are notable only in the upper 1000 dbars, but indistinguishable from 
zero everywhere. There is a little more noise in the deep layers of the Rio and Tristan da 
Cunha boxes (Fig. 3-10f, g), probably due to enhanced time variability there. The max¬ 
imum adjustments in the net fluxes in the surface layer from the initial guess (including 
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the Ekman transport) is in the 20%-30% range (Table3.11). 

Absolute velocities and transports integrated from the west are given in Appendix D 
(pages 220-235). The NADW flow does not appear to take place exclusively in a well- 
defined deep western boundary current across every section (e.g. at 24°N, page 223), as 
suggested by ocean circulation theory and current meter measurements. This was already 
pointed out by other authors who used different sections (Macdonald, 1998; Rintoul and 
Wunsch, 1991). It is not clear whether such an interior flow is an accident of a particular 
section time or reference surface choice. Nevertheless, in my model, net integrated prop¬ 
erty fluxes appear to be insensitive to the details of the horizontal circulation. That is, 
forcing the velocity field with results from a current meter array can change the horizontal 
structure as was done at A2, but does not necessarily affect the integrated fluxes. 

The zonally integrated mass flux streamfunctions of the standard circulation are shown 
on Figure 3-11 for each section. The net fluxes (numbers below the plots) are indistin¬ 
guishable from the Bering Strait outflow but are much more uncertain because of the 
baroclinic variability in the upper layers. Integrated transports over chosen sets of layers 
are indicated on Table 3.14. (The layers of integration do not correspond exactly to the 
approximate water masses defined above. Rather, the flow is summed over consecutive 
layers for which the integrated transport is flowing in the same direction.) 

AABW 

A net northward flux of AABW (defined as the total flow in the northward-flowing bottom 
layers) is indicated in the South Atlantic, gradually decreasing from 6±1 Sv at 30°S to 
l±lSv at 24°N. At 45°S, the AABW flux is more uncertain with 5 ± 2Sv. Within my 
error bars, the estimated transport in the North Atlantic is consistent with the Schmitz 
and McCartney (1993) circulation schematic (their Figure 13). In the South Atlantic, the 
AABW flow is slightly smaller than but fully consistent with that suggested by Hogg et 
al. (1999). 
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Figure 3-10: Standard solution mass residuals in the Atlantic Ocean. The line gives 
the residual in invidual layers (in 10 9 kgs _1 ) and the shaded area give the one standard 
deviation uncertainty. The top-to-bottom residual is indicated below the graph. 
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Figure 3-11: Atlantic Ocean standard circulation. The mass transport across each section 
is shown, integrated from the bottom (the overturning streamfunction). The shaded area 
corresponds to the uncertainty (one standard deviation). The net transport is indicated 
below each graph. While panel c) gives the total transport in the Florida strait alone, 
panel d) gives the total transport across 24°N. Note the different scales on c). 
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Latitude 

“AABW” 

“NADW” 

“Intermediate+thermocline” 

48°N (A2) 

0 

-16±2 

19±2 

24° N (A5) 

1±1 


14±2 

7.5°N (A 6 ) 

3.6T1.2 


1±4 

4.5°S (A7) 

3.4T1.2 

-29±4 

42±7 

11°S (A 8 ) 

3.1T1.5 

-23±4 

28±4 

19°S (A9) 

3.4±1.1 

S£H 

26±3 

30°S (A10) 

6±1.3 

■eidsb 

18±3 

45°S (All) 

5±2 

-18±4 

8±3 


Table 3.14: Atlantic “water mass” transports (Sv). Transports are summed over consec¬ 
utive layers of same integrated flow direction (see Fig.3-11), which does not correspond 
exactly to the approximate water mass definitions of Tables 3.9 and 3.10. Ekman transport 
is not included. 


NADW 

The NADW flow increases to the south, from -16 Sv at 48°N to -23 Sv at 30°S (Table3.14). 
The increase in NADW transport can be explained by entrainment of AABW and AAIW 
(i.e., Schmitz and McCartney, 1993). The NADW flow at 45°S is marginally weaker and 
more uncertain, at 18±4 Sv. 

The present NADW transport estimates across the different sections are generally 
consistent with previous studies, e.g., Wunsch and Grant (1982, 14 Sv at 48°N; 18 Sv 
at 24°N); Roemmich and Wunsch (1985, 24°N, —20±5Sv); Fu(1986, South Atlantic, 
-20Sv); Speer et al. (1996, —14Sv between 7 " = 27.5and 7 ” = 28.11 at 11°S, versus 
—20±4Sv here); Holfort and Siedler (1999, 30°S, — 23 Sv); and Saunders and King (1995, 
45°S, — 22 Sv). Transports are also, within my error bars, relatively close to the Schmitz 
and McCartney (1993) and Schmitz( 1996a) schematics. 

Koltermannet al. (1999) estimated variable NADW transport at 48°N between -15 and 
-20 Sv and at 24°N from -17 to -24 Sv. But the model that they used relied upon Sverdrup 
balance in the interior to compute the reference velocity adjustment. Luyten et al. (1985) 
suggested that Sverdrup balance was not a good approximation in the northern Atlantic. 
It is additionally difficult to distinguish between the parts of their circulation that is due to 
model configuration (e.g., reference surface choice) and the parts due to actual variability 
in the density field. 
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The NADW transports agree quantitatively with the global Macdonald (1998) esti¬ 
mates everywhere except at 48°N and 11°S. At 48°N, Macdonald’s large NADW transport 
of —26 Sv could be due to temporal variability (e.g., Koltermann, 1999) or to a bottom tri¬ 
angle problem. (Transports are sensitive to the bottom triangles across this section when 
the reference surface is taken at the deepest common level of station pairs, see Figure 2-9, 
page 53). At 11°S, I do not know what is responsible for Macdonald’s large —28 ± 1 Sv 
transport-possibly time variability. Nevertheless, her high transport figure is within the 
present uncertainties. 

At 4.5°S, the NADW flux is the largest (29±4Sv). From this same section, Lux(1997) 
obtained a similar flow associated with large dianeutral transfers between 7.5°N and 4.5°S 
that accounted for the difference in the overturning structure between the two latitudes. 
Geostrophic transports computed at 4.5°S are most subject to noise because of the small 
Coriolis parameter there, i.e., Chapter 2). 

I did not attempt to compare in detail the deep flow with Reid (1994) because he gives 
transport values integrated from the upper limit of NADW all the way to the bottom, thus 
including AABW. But his qualitative picture of the circulation, with a large production of 
NADW (29 Sv flowing South at 60°N), and a decrease of NADW flow to the south clearly 
conflicts with the standard circulation. 

Intermediate and surface flow 

Integrated intermediate and surface flows are indicated in Table 3.14, third column (Ekman 
transport is not included). Discussion of those flows is, by conservation of mass, equivalent 
to discussing the combined NADW and AABW transport. The total northward flow in the 
upper layers, including the Ekman flux, is in the range of 14 to 20 Sv outside the equatorial 
box, thus compatible with most authors’ estimates, including Macdonald (1998, Tables 12 
and 14); Schmitz and McCartney (1993) and Schmitz (1996a). The vertical structure of 
the flow (Fig. 3-11) generally indicates an equipartition between AAIW and surface flow, 
contrary to the results of delasHeras and Schlitzer (1999), who suggest that intermediate 
water flow dominates the return flow into the North Atlantic. 
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Boundary Currents 


The horizontal flow structure in the density classes defining the different water masses 
defined in Tables3.9-3.10 is given in AppendixD (FiguresD-26-D-40). Stations with 
common flow direction were selected to estimate boundary current transports. Again, the 
formal uncertainty of the boundary current transports is underestimated because the error 
budget did not take into account baroclinic variability on small scales. 

In the North Atlantic at 48°N, there is a northward flow between 50°W and 43°W of 
47±2 Sv of waters above 7 n = 27.82 (9 — 5°). This flow, composed of the North Atlantic 
Drift plus recirculating currents, is larger than the Macdonald (1998) estimate (about 28 Sv 
from her Fig. 5) and the Schmitz and McCartney (1993) sketch that suggests 37 Sv. 

At 24°N above y n = 27.82 (9 = 5.4°), the upper part of the Antilles Current flows 
against the Bahamas Bank, with 27±3 Sv between 72°W and the Bank. This flow is more 
than twice that estimated by Macdonald (1998) and Schmitz and McCartney (1993), but 
is present in the baroclinic field before inversion. Current meter observations to the north 
at 26.5°N indicate a flux of about 5Sv (Lee et al., 1996), but over a different longitude 
range (to my knowledge, there are no current meter data along the western part of the 
24°N WOCE section). 

The A 6 section ( 6 °N in the western end) is extremely noisy against the coast of Brazil. 
A northward flow of 5±1.5Sv is found through the first station pair, that totals 10 Sv 
when the North Brazil Current on the shelf is included. But east of the first station pair, 
the eddy field dominates the transport signal, with transport variations of ±8 Sv from one 
station to the next (above 7 ” = 27.82, 9 — 4.5°). At 11°N, Macdonald (1998) obtained a 
northward flow (27 Sv) against the coast. But at that latitude, the North Brazil Current 
lies against the coast, in contrast with the highly variable retroflexion region at 6 °N (Johns 
et al. 1990). 

Farther south along the coast, the North Brazil Current can be identified with 44±7 Sv 
flowing north at 4.5°S (above 7 ” = 27.82, 9 = 4.2°, west of 34°W) and 18±4Sv at 11°S 
(above 7 " = 27.5, 9 = 3.9°, west of 34.7°W). The transport estimate at 11°S is similar 
to that obtained by Speer et al. (1996) from the same section, but only half the 40 Sv 
obtained by Macdonald (1998) from earlier data. 
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At 19°S and 30°S, one finds the southward-flowing Brazil Current, with a strength of 
-15 ± 3 Sv at 19°S (above 7 n = 27.5, Q = 3.6°, west of 36.9°W) and -14 ± 4Sv at 30°S 
(above 7 ” = 27.5, 6 = 3.5°, west of 40° W), or —9±3 Sv west of 46.5° W, depending upon the 
boundary current definition. Macdonald (1998) obtained a transport of similar magnitude 
at 23°S (—13±1 Sv). My estimate at 30°S is lower than the current meter results of Muller 
et al. (1998, —16 Sv at 30°S and west of 45° W) and Zemba (1991, —22 Sv against the coast, 
based on earlier hydrographic data). It is not clear whether the difference is significant. 
There is little correction to the initial (relative) flow field in this region (Fig.D-38), so 
the difference with Zemba likely reflects real oceanic variability. At 45°S, one identifies 
the Falkland Current, with 30 ± 10 Sv flowing north (top-to-bottom, west of 58.4°W), 
consistent with Saunders and King (1995, 40 Sv). 

Along the eastern boundary in the South Atlantic, the Benguela Current flows north¬ 
ward at 30°S (A10 and All). Across A10, the strength of the flow is of 14 ± 2Sv east 
of 11.8°E. Across All, it is of 18±1.6Sv east of 12.3°E (both above 7 " = 27.5, 0 = 3.5°, 
that is, including AAIW). This magnitude is in agreement with independent estimates 
(e.g., Shannon and Nelson, 1996, 15 Sv at 30°S; Macdonald, 23 ± 10 Sv at 27°S, Garzoli et 
al., 1996, 16 Sv at 30°S). 

In the South Atlantic, the present boundary current estimates are thus in agreement 
with estimates from the literature that are based on the same section or on modern current 
meter arrays. However in the North Atlantic, substantial differences with previous results 
are found. The differences may be simply due to changes in the circulation structure due to 
natural oceanic variability. But again the actual transport uncertainties (particularly those 
due to time variability) are much larger than the ones quoted here because of unaccounted 
baroclinic variability on small scales and the intense eddy field in western boundary regions 
that makes transport estimates extremely sensitive to the choice of integration boundaries. 

Dianeutral transfers and diffusivity 

Dianeutral transfers are shown on Fig. 3-12, upper panels. Upwelling in the bottom layers 
is found almost everywhere. At deep levels, upwelling is found mainly in the Equatorial 
(4.5°S-7.5°N), Rio (30°S-19°S) and Tristan da Cunha (45°S-30°S) boxes, totaling 5-7Sv, 
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Figure 3-12: Dianeutral transfers (upper) and diffusivities (lower) in the Atlantic Ocean. 
The line indicates the transfer or diffusivity at the bottom of the layer. 
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depending upon depth. At intermediate levels (1000-1500 dbars), upwelling is found in the 
Tristan da Cunha (45°S-30°S, Fig.3-12f) box while downwelling is found in the Equatorial 
box, with —10 ±2 Sv at 1300 dbars. This pattern is consistent with the NADW southward 
increase from entrainment of AABW and to some extent AAIW. In the other boxes, 
transfers are more uncertain at intermediate and deep levels. 

In the surface layers, a large downwelling is found over the North Atlantic box, with 
—13 ±3 Sv, then upwelling in most other boxes, except in the Tristan da Cunha box, — 5± 
2 Sv downwelling. In the North Tropical and Rio boxes, the Ekman transport convergence 
would suggest a downwelling (in Cartesian coordinates). The geostrophic flux in the 
upper layer and possibly horizontal advection across outcropping isopycnals are thus strong 
enough to cancel and reverse this tendency. 

Similarly, Macdonald (1998) found a consistent upwelling over the deep Equatorial 
layers and downwelling of AAIW. In the South Atlantic, her results are very noisy while in 
the North Atlantic she obtained upwelling, where my estimate does not produce significant 
dianeutral transfers. 

Dianeutral diffusivities in the deep layers are positive and systematic in most boxes 
(Fig. 3-12, lower panels). Their values are in the range 0-50 cm 2 s -1 , depending upon depth 
and location. The highest diffusivities are found in the bottom layers (up to 50±50 cm 2 s _1 
in the North Atlantic), where the area of the layer interfaces are small-and subject to 
uncertainties in their calculation from Levitus (1994). In the North Tropical, Equatorial 
and Rio boxes, diffusivities are significant and in the range 5-20cm 2 s _1 . The average 
diffusivity in the Rio box is of 4.6 ±3 cm 2 s -1 , which is similar to the 1-5 cm 2 s _1 calculated 
by Hogg et al. (1982; 1999a) from a heat budget in the deep Brazil Basin. In the Equatorial 
and Rio boxes, highest diffusivities are found approximately at the depth of the crest of 
the Mid-Atlantic Ridge, over which the sharp topography is known to induce enhanced 
mixing (Polzin et al., 1997). 

Fig. 3-13 shows the average diffusivities and dianeutral transfers in the central Atlantic 
superbox between 24°N and 19°S (see above). The average diffusivity in this superbox 
below 1000 dbars is 7 ± 2cm 2 s _1 ). There is no evidence for enhanced dianeutral mixing 
above 1000 dbars in the Atlantic Ocean; estimates are within the range allowed by the 


123 




24N-19S,k z 



W fiux 



Figure 3-13: Diffusivities (left) and dianeutral transfers (right) in the central Atlantic box 
(24°N-19°S). 


uncertainties-less than l-2cm 2 s _1 . 

Dianeutral velocities and diffusivities are thus well resolved in the Atlantic, with av¬ 
erage values for the diffusivities generally within the range of those previously reported. 
Basin-average diffusivities and dianeutral velocities are given at the end of this Chap¬ 
ter (page 138). 


3.3.4 Atlantic circulation summary 

The Atlantic circulation estimate is relatively well determined, compared with the larger 
uncertainties in other oceans. In particular, uncertainties of derived water mass are rela¬ 
tively small. North Atlantic Deep Water flows southward at —16±2 Sv at 48°N; it entrains 
Antarctic Bottom Waters and Antarctic Intermediate Waters on its way to the south, re¬ 
sulting is a flux of -23 ± 3 Sv at 30°S. A northward flux of 6 ± 1.3 Sv of Antarctic Bottom 
Waters is found to enter the Atlantic at 30°S west of the Mid Atlantic Ridge, progressively 
decreasing to 1 ± 1 Sv at 24°N and, OSv at 48°N. The flux at 45°S is slightly different, 
although indistinguishable within uncertainties from the 30°S flow. The circulation in the 
intermediate and upper layers provides the 14 to 17 Sv flow of warm water to balance the 
deep water export. This circulation is generally consistent with the existing literature. 
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However, pointwise differences were found, both in integrated transports and horizontal 
structure. A large part of those differences is probably attributable to natural oceanic 
variability; but because of the different methods used, it is difficult to determine how 
much. A possible extension of this work is to produce, with the same model, an estimate 
of the circulation based on the earlier data. 

The strength of the Benguela Current that I diagnosed has about the same magnitude 
as independent estimates but, as is the case for most boundary currents, its transport 
is very uncertain due to an enhanced eddy field. As a result, no attempt was made to 
distinguish between the provenance of warm and cold waters, i.e., “warm water path” 
versus “cold water path” (i.e., Gordon, 1986; Rintoul, 1991). 

The horizontal structure of the flow was analyzed succinctly and shows that locally, the 
flow direction is highly variable and generally does not correspond to the one suggested by 
water property distribution. However, the direction of the integrated flow does correspond 
roughly to that suggested by traditional water mass analysis. 

Consistent upwelling, with significantly non-zero diffusivities, were found in the deep 
layers, with k * values in the range 2 to 5cm 2 s _1 over most of the Atlantic. Strongest 
mixing and upwelling occur over different depth ranges, depending upon the region. 
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3.4 The Southern Ocean 


3.4.1 Introduction 

The Southern Ocean is the highway of the global circulation. In addition to the one¬ 
way eastward flow in the Antarctic Circumpolar Current of about 140 Sv (Whitworth et 
al., 1982), there are major interbasin exchanges between the Indian and Atlantic Oceans 
(e.g., Georgi and Toole, 1982; Gordon, 1986). Recent estimates of the Southern Ocean 
circulation from inverse calculations based on hydrographic sections include the work 
of Macdonald (1998, global), Sloyan and Rintoul (1999a,b, south of 12°S); Gille (1999, 
southeastern Pacific); and Rintoul (1991, Southern Atlantic). 

Because most layers are in contact with the atmosphere, the Southern Ocean is a place 
of important dianeutral transfers, that is, of water mass transformation (e.g., Macdonald, 
1998; Sloyan and Rintoul, 1999a; Schmitz and McCartney, 1993). In the present data set, 
all layers that are found above 4000dbars north of 30° S outcrop in the Southern Ocean. 
In that view, a large part of the diapycnal transfers in the global ocean could take place 
in the layer outcropping regions in the Southern Ocean. The Antarctic Bottom Waters 
(AABW), the densest waters in the global ocean, are formed along the shelf of Antarctica. 
The Weddell Sea has been in the past regarded as the main formation site. However, recent 
observations indicate formation in various places around Antarctica including the Amery 
Basin, Adelie Land, Ross Sea and western Weddell Sea (Fig. 3-14), e.g., Rintoul (1998); 
Sloyan and Rintoul (1999); Orsi et al. (1999) and references therein. However, the total 
amount of AABW being produced is poorly known. Existing literature cites a range of 7 to 
20 Sv. But estimates based on either heat budgets (Gordon, 1975) or inverse models forced 
by existing surface flux estimates (Sloyan and Rintoul, 1999a,b) indicate a much higher 
production rate of O(40-50Sv). The bottom water production rate will be discussed in 
light of my global model results. 
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3.4.2 Model constraints 


Data and treatment 

Figure 3-14 shows the Southern Ocean meridional lines, namely, A21 (Drake, Jan. 1990, 
Roether et al., 1993); A12 (May 1992, Lemke, 1994); 16 (Feb. 1996, Park et al. 1999); 
I9S (Dec. 1995, Hufford et al., 1997); P12 (Jan 1995, Rosenberg et ah, 1997) and P14S 
(Jan. 1996, McTaggart and Johnson, 1997). The hydrographic data were provided by 
the Pis W. Roether (A21); P. Lemke (A12); Y.-H. Park, N.Metzl and C. Jeandel (16); M. 
McCartney and A. Gordon (I9S); S.Rintoul (P12); J.Bullister and G. Johnson (P14S). 

Layers and constraints 

The Southern boxes (Fig. 3-14) are respectively the Georgia-Weddell box (A21-A11-A12); 
the South Agulhas box (A12-I6); the Southern Indian box (I6-I5-I9S); the Southern Aus¬ 
tralia box (I9S-P12); and the South Pacific box (P12-P6-A21). A box called the “Mac¬ 
quarie” box defined by P12, P14S and the western part of P 6 was explored, but the 
conservation constraints were weaker there because P14S did not reach the Antarctic con¬ 
tinent. Moreover, this box was not closed between P 6 and New Zealand, leaving a possible 
leak down to 2000 dbars to the north. An additional super box was defined in the Southern 
Indian between A12-I5-I9S. 

The neutral surfaces defining water masses were again based on the Macdonald (1998) 
isopycnal surfaces (Table3.16). But in the Southern Ocean, the formal separation be¬ 
tween AABW and CDW required a somewhat lighter neutral density, 7 " = 28.27, that 
corresponds to the sill depth of Drake Passage (i.e., Orsi et al., 1999). The initial Ek- 
man transports are relatively small (Table 3.16) and their uncertainty was given the lower 
bound of ±1 Sv. 

Because of the possible freshwater input to most layers as they outcrop, the a priori 
noise in the salt anomaly equations was increased to ±(35g/kg) x(lSv). No heat and 
PO conservation constraints were imposed. Thus only mass and total silica had an effec¬ 
tive influence. Consequently, the uncertainties on the diagnosed dianeutral transfers and 
mixing are large relative to other oceans. 
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Figure 3-14: Southern Ocean sections and topography. The 4000 m isobath is contoured 
and areas shallower than 3000 m are shaded. Box names are indicated in the shaded 
rectangles. 
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II 

0 
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II 
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O 

0 

14 

7 n = 28.27 
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15 
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II 

1 
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Table 3.15: Southern Ocean neutral surface layers 7 " (in kgm -3 ). The apparent tem¬ 
perature reversal between layers 7 and 8 is due to outcrop and salinity effects. Water 
masses: AAIW=Antarctic Intermediate Waters; UCDW/LCDW= Upper/Lower Circum¬ 
polar Deep Waters. AABW=Antarctic Bottom Water. 


The net flux through the Drake Passage (P21) was constrained to 140 ± 6 Sv east¬ 
ward, according to the current meter and hydrographic transport estimate of Whitworth 
et al. (1982) and the error budget of the preceding chapter. (140±5Sv corresponds to a 
weighted average of their estimates with respective uncertainties, i.e., Macdonald (1995); 
and an extra Sverdrup, rms, comes from the error budget.) The bathymetry in the South¬ 
ern Ocean does not suggest imposing additional transport constraints. 


Section 

Initial 

Final 

% correction 


0.06T1 

MEM 

- 


-1.3±1 


1 % 

16 

0 . 6±1 



I9S 

-0.9±1 

-0.9±1 

- 

P12 

Oil 

0i0.5 

- 

P14S 

-0.5il 

-0.46±1 

- 


Table 3.16: Southern Ocean initial Ekman transport and correction (in Sv, positive east¬ 
ward 
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Reference surfaces and dianeutral terms 


All reference surfaces were taken at the bottom (last common depths of station pairs), that 
is, where the velocities are likely to be the weakest. As for the other oceans, a large range 
is allowed in the dianeutral transfers (±10 -3 cms -1 and ±100cm 2 s -1 ) because little is 
known about their actual size for outcropping layers. Moreover, the diffusivity formalism 
may not appropriately represent water mass modification by air-sea exchange. 

3.4.3 Southern Ocean Circulation 

All residuals of the equations of the “standard” solution were within the a priori range. 
Mass residuals in individual layers are tiny and all are indistinguishable from zero (Fig. 3- 
15). But the uncertainties of the residuals are close to the a priori values, reflecting the 
low weight of the salt equations and the absence of heat and PO constraints. Residuals 
in the Macquarie box indicate a systematic convergence, with —5 ± 7 Sv in total, possibly 
due to an eastward flux between P14S and Antarctica (or just noise). Because the wind 
is nearly perpendicular to the sections, the Ekman transports (and their corrections) play 
a minor role in the Southern Ocean (Table 3.16). 

Absolute velocities and water-mass transports cumulated from the north are given in 
Appendix D (pages 236-247) for reference. The integrated flow for the standard solution 
of the Southern Ocean circulation is shown on Figure 3-16. The net fluxes through the 
“choke point” sections are close to the Drake Passage flow, 140±6 Sv eastward across A21. 
South of Australia, the net flux is increased by about the size of the Pacific-Indonesian 
throughflow, adding up to 154±9 Sv, which is similar to the regional calculation of Rintoul 
and Bullister (1999). Above 3000dbars, the transports increase in a quasi-linear way, 
giving a parabolic shape to the streamfunction. Macdonald (1998) obtained a similar 
structure based on earlier data. 

Dianeutral transfers and diffusivity 

Dianeutral mass transfers are indicated on Fig.3-17. (In the Weddell region, the first 
six layers outcrop north of the ACC.) Because of the low property conservation weights, 
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Figure 3-15: Mass residuals in the Southern Ocean boxes. 
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Figure 3-16: Southern Ocean standard circulation. The mass transport across each section 
is shown (xlO 9 kgs -1 ), integrated from the bottom. The shaded area corresponds to the 
uncertainty (one standard deviation). The net transport is indicated below each graph. 
The meridional slope of the layers are steep (all layers outcrop), and the depth on the 
graph indicates the cumulated average thickness of the layers. 
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dianeutral diffusivities are uncertain. In the Weddell sea (Fig.3-17a), there is a weak 
indication of a water mass transfer from dense water to light water of about 6 ± 7 Sv near 
the surface (Hereafter I label such a transfer “upwelling” for convenience, but bearing in 
mind that the water may as well be transferred horizontally by warming at the outcrop.) 
Similar upwellings are observed in the “Agulhas” box (Fig.3-17b) and in the Southern 
Pacific (Fig.3-17e, 13 ± 5Sv). However, there is dense water production, or downwelling, 
in the Southern Indian and Southern Australia boxes in this same density range (Fig.3- 
17c, d) so that the resulting net diapycnal flux through -y n = 27.72 in the Southern Ocean 
is not significantly different from zero at 6 ± 9 Sv. 

The surface 7 " = 27.72 corresponds to the upper part of NADW in the South At¬ 
lantic. Conversion of NADW into AAIW has been suggested in the past (see Schmitz and 
McCartney, 1993, for a review) and more recently by numerical simulations (Shriver and 
Hulburt, 1997). However, neither the inverse box models of Macdonald (1998), Sloyan and 
Rintoul (1999a,b) produced net upwelling across that interface. Here, there an upwelling 
is found, although not significant. 

In deeper layers, deep and bottom water formation is inferred in the Weddell Sea with 
—14 ± 11 Sv crossing 7 n = 28.15. In comparison, Rintoul (1991) estimated AABW pro¬ 
duction at 9 Sv based on earlier data. Deep transfers are very uncertain in all other boxes. 
However, because the uncertainties are correlated, the total southern deep water formation 
(including entrainment) across 7 ” = 28.15 is significant, at —20 ± 6 Sv. While Macdon¬ 
ald (1998) did not quantify the bottom water formation rate, Sloyan and Rintoul (1999a) 
obtained a comparable value across 7 " = 28.20. However, their water mass transfer just 
above-across 7 " = 28.0, was much larger (50 Sv). Such a large deep water production 
would imply a much larger deep inflow into the Pacific and Indian basins than the one 
that is permitted by this model, as discussed in Sections 3.1 and 3.2. 

Existing estimates of bottom water formation in the Southern Ocean bear large uncer¬ 
tainties. The high Sloyan and Rintoul estimate may be related to the COADS heat flux 
forcing they adopted, and which is very uncertain in the Southern Ocean. Other water 
mass formation estimates are available from tracers and CFC measurements [e.g., Jacobs 
et al. (1985), Orsi et al. (1999)], but those depend upon uncertain assumptions about the 


133 



entrainment rates. 

Orsi et al. (1999) estimated the total Southern Ocean production of water below 7 " = 
28.27 at 9Sv (providing no error bars). But because 7 ” = 28.27 intersects the bottom 
south of the northern boundaries of my boxes, the net integrated dianeutral transfer across 
7 " = 28.27 is null in my model. 

In conclusion, my estimate of the Southern Ocean bottom water formation is in agree¬ 
ment with most previous estimates-given the large existing uncertainties. However, I do 
not recover the large dianeutral fluxes obtained by Sloyan and Rintoul (1999a,b) which 
would be incompatible with my circulation at lower latitudes. Dianeutral diffusivities in 
the Southern Ocean are not resolved because of the weak weights on conservation equa¬ 
tions, other than mass. 
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Figure 3-17: (continued) 
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Figure 3-17: Dianeutral transfers (upper) and diffusivities (lower) in the Southern Ocean. 
The line indicates the transfer or diffusivity at the bottom of the layer. Again, layer depth 
is just an average. 



























3.4.4 Southern Ocean circulation summary 


The zonal circulation in the Southern Ocean is dominated by the Antarctic Circumpolar 
Current, with a (constrained) flux of 140 ± 6Sv at Drake Passage. The analysis of di- 
aneutral fluxes gives estimates for the deep and bottom water transformations that are 
consistent with my global circulation estimate. While the conversion of NADW into AAIW 
is rather uncertain with 6 ± 9 Sv, 20 ± 6 Sv of lower NADW is transformed into bottom 
water. According to Levitus (1988), northward Ekman transports in the range 15-25 Sv 
are found around latitude circles in the Southern Ocean, within the range my 6 ± 9 Sv 
water mass conversion estimate. But it is not clear how the Ekman transport, the non¬ 
zonality of isopycnal outcrops and other oceanic processes interact to produce the water 
mass conversion. 
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3.5 Global circulation summary 


In a very idealized summary of the integrated deep water fluxes, my best estimate of the 
average overturning circulation for the period 1985-1996 can be described as follows. A 
volume of 16±2 Sv of North Atlantic Deep Water was being produced in the northern North 
Atlantic and moved southward. During its transit through the Atlantic Ocean, Antarctic 
Bottom Water and Antarctic Intermediate Waters were entrained, increasing the NADW 
transport to 23±3 Sv as it exited the South Atlantic at 30°S. In the Southern Ocean, a total 
of 22 ±8 Sv of AABW was formed from lower Circumpolar Deep Water-which corresponds 
approximately to the lower NADW density range. Bottom water inflows (NADW+AABW 
mixture) to the Atlantic, Indian and Pacific Oceans are respectively 6 ± 1.3 Sv, 11 ±4Sv 
and 8 ± 2 Sv. In the Indian and Pacific Oceans, most of this water returns southward at 
deep and intermediate levels. 

Due to enhanced variability, the circulation in the upper layers is more uncertain. A 
northward flow of 18 ± 3 Sv of thermocline water from the South Atlantic Ocean bal¬ 
ances the NADW southward flow; and the strength of Pacific-Indonesian-Throughflow is 
estimated at 15 ± 5 Sv, which is in the higher range of previous estimates. 

Generally speaking, the circulation in the Atlantic Ocean is remarkably close to esti¬ 
mates obtained from previous data, suggesting little change in integrated mass transports 
over the past 30 years. (However, some of the changes or non-changes may be model- 
dependent, and a more stringent test would be to produce an estimate of the circulation 
from earlier data using the same treatment and model as used here.) The deep inflow into 
the Indian Ocean is consistent with previous studies as well. My standard model estimate 
of the inflow in the South Pacific Ocean is in the lower range of previously published val¬ 
ues and depends upon the weight of the PO constraints. Non-conservation of PO yields 
a higher inflow but produces larger mass residuals in the deep Pacific layers. Closely re¬ 
lated, the deep inflow to the North Pacific Ocean is weaker than previously believed, as 
a consequence of heat and salt conservation in the northern parts of those basins, which 
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W* at top (Sv) 

w* (xlO 4 cms *) 

k* (cm 2 s *) 

Atlantic bottom 

5±2 

0.6±0.5 

7±6 

Indian bottom 

4±2 

0.5±0.5 

12±16 

Pacific bottom 

9±2 

0.4±0.2 

14±4 

Southern bottom 

- 20±6 

-0.4±0.3 

- 



0.2±0.05 

3±2 

Indian deep 

9±4 

0.26±0.15 

3.5±2 

Pacific deep 

7±3 

0.08±0.06 

4.4±1.3 

Southern deep 

6±9 

0 ± 0.2 

- 


Table 3.17: Basin-averaged dianeutral velocities and diffusivities. The average is calculated 
on neutral surfaces from 7 71 = 28.1 to the bottom (generally 3800 dbars to the bottom) 
for the “bottom layers” and from y n — 27.96 to 7 ” = 28.07 for the “deep layers”. The 
second column indicates the total transport in the basin through the top interface, while 
the third and forth columns indicate the average dianeutral velocities and upwelling. 


were not generally considered in previous models. 

Associated with the circulation, significant dianeutral exchanges were found. Table 3.17 
summarizes the dianeutral exchanges basin-wise in the deep and bottom layers. Dianeutral 
velocities are in the range 0.1 — 0.6 x 10 -4 cms ' 1 while the diffusivities are in the range 
3 — 4cm 2 s ~ 1 in the deep layers, and larger near the bottom. (Diffusivities were not 
resolved in the Southern Ocean, where isopycnals outcrop.) Bottom water formation in 
the Southern Ocean could be diagnosed consistently with the global circulation, implying a 
dianeutral velocity of —0.4±0.3 x 10 -4 cms -1 . The global average dianeutral velocity over 
the Atlantic, Pacific and Indian Oceans North of 30°S is of w* = 0.56 ± 0.36 x 10 _ 4 cms _1 
in the bottom layers and of w* = 0.16 ±0.04 x 10 _ 4 cms _1 in the deep layers. The average 
diffusivity is of k* = 9 ± 5 cm 2 s -1 in the bottom layers and of k* = 3.4 ±1.3 cm 2 s -1 in the 
deep layers (these averages were not weighted by the area, however). 

This simplified circulation scheme is based on integrated mass transports calculations, 
and, as we will see in the next Chapter, does not always directly determine the actual 
property transports. 
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Chapter 4 


Global Property Fluxes and 
Divergences 


In this chapter, fluxes of freshwater, heat, nutrients and oxygen as diagnosed from the 
circulation scheme are analyzed. The analysis starts with the freshwater flux diver¬ 
gences (Section 4.1); energy flux estimates are then discussed (Section 4.2) followed by 
silica (4.3.1), nitrate and phosphate (Section 4.3.2) and oxygen (Section4.4) transports. 
Transports and residuals of nutrients and oxygen are discussed in relation to biological 
activity. In particular, the residuals allow an estimate of the export production. Net 
property transports across sections are listed in Appendix E. 

4.1 Freshwater flux divergences 

Freshwater divergences are calculated as part of the solution (Eq. 1.8 and 1.9, page 20) 
through the anomaly formulation, subject to the assumptions about the errors described 
in Chapter 2. Inferred precipitation plus runoff minus evaporation (P-E) in the standard 
solution are given for the different oceans and compared with other authors in Tables 4.1 
to 4.3. 

In the Indian Ocean (Table 4.1), freshwater divergences are indistinguishable from 
zero in the North Indian Ocean sector, north of 8°S (12-110, 0.1±0.3 Sv net precipitation). 
South of 8°S, the solution shows net evaporation in the tropical (20°S-8°S) and subtropical 
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Indian 

Between 


20° S 

8 °S 

and 


8 °S 

Coast 

Present study 

-0.35 ± 0.25 

-0.3 ±0.2 

±0.1 ±0.3 

Jourdan et al. (1997) 

Q2U 



Wijffels et al. (1992) 




Macdonald (1996) 

-0.4 



Oberhuber (1988) 

<-0.4 —> 

±0.25 


Table 4.1: Precipitation (+) minus evaporation (-) estimates for the different oceanic 
regions bounded by the model sections in the Indian Ocean. Units are in Sv. 



Pacific 

Between 


17°S 

24° N 

47°N 

and 

17°S 

24° N 

47° N 

Coast 

Present study 

-0.3 ±0.26 

±0.1 ±0.3 

±0.14 ±0.2 

±0.1 ±0.7 

Jourdan et al. (1997) 

-0.5 

-1.2 

0 

±0.3 

Wijffels et al. (1992) 

-0.2 

±0.3 

0 


Macdonald (1996) 



0 


Tsimplis et al. (1998) 

-0.4 




RMC89 



0 



Table 4.2: Same as Table 4.1 but for the Pacific Ocean. RMC89=Roemmich and McCal- 
listair (1989) 

boxes (32°S-20°S) (a minus sign denotes net evaporation). These two latter estimates are 
consistent with the estimates from satellite radiometric data of Jourdan et al. (1997), (Ta¬ 
ble 4.1, error bars are not given for these latitude bands, but they probably do not exceed 
100%). Wijffels et al. (1992) , Macdonald (1995) and Oberhuber (1988) report estimates 
in the same range. (Note that the Wijffels et al. (1992) is essentially an integration of 
the Baumgartner and Reichel (1975) estimate.) Oberhuber’s (1988) values have a greater 
contrast north and south of 10° S, but are within our error bars. Overall, the Indian 
Ocean is evaporative, with —0.6 ± 0.4 Sv net air-sea exchange. As for all property fluxes, 
the model’s large Pacific-Indonesian Throughflow may affect substantially the estimated 
freshwater divergences over the Indian Ocean. 

In the South Pacific (Table 4.2), net evaporation characterizes the subtropical region, 
in agreement with the estimates of Jourdan et al. (1997) and the inversion of Tsimplis et 
al. (1998). It encompasses the Wijffels et al. (1992) estimate within uncertainties. North 
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Atlantic 


48° N 

24° N 

7.5°N 

4.5°S 

ITS 

19°S 

30° S 


24° N 

7.5°N 

4.5°S 

11°S 

19°S 

30° S 

45° S 

Present 

0 ± 0.15 

-0.1 ±0.3 

0.1 ±0.4 

0.2 ±0.5 

-0.3 ±0.3 


0.45 ± 0.2 

JPG1997 

-0.2 

-0.4 

±0.2 

-0.15 

-0.25 

-0.25 


SBD1989 

-0.3 

-0.2 






WSB1992 

-0.3 

-0.2 

±0.2 

-0.1 

-0.2 

-0.2 


M1995 

-0.1 

0 

<— ±0.15 —■* 

-0.4 




Table 4.3: Same as Table 4.1but for the Atlantic Ocean (JPG1997= Jourdan 
et al. (1997); SBD1989=Schmitt et al. (1989); WSB1992= Wijffels et al. (1992); 
M1995=Macdonald (1995)) 


of 17°S, freshwater divergences are more uncertain. Our result in the Central Pacific 
is incompatible with the large evaporation found by Jourdan et al. (1997)-although the 
uncertainty on the later is probably at least 100%. In contrast, the model estimate for 
the Pacific Ocean north of 17°S is for net precipitation: a total of 0.3 ± 0.3 Sv. The net 
freshwater budget over the Pacific Ocean north of 32°S is balanced within the uncertainties, 
0 ± 0.4 Sv. 

The North Atlantic Ocean freshwater flux divergences are weak and uncertain (Ta¬ 
ble 4.3). Our estimates agree with Macdonald’s (1995) within one standard deviation. 
The evaporation—precipitation magnitudes reported by Jourdan et al. (1997), Schmitt et 
al. (1989) and Wijffels et al. (1992) are larger, although within two standard deviations of 
our estimate. In the South Atlantic, freshwater flux divergence estimates are larger-but 
still encompass zero within an uncertainty of two standard deviations. Both size and mag¬ 
nitude are consistent with previous estimate. Over the central Atlantic between 24°N and 
19°S, the model estimates a net evaporation of -0.1±0.2 Sv, in agreement with the Wijffels 
estimate (-0.2 Sv), but smaller than the radiometric measurements (-0.55 Sv). Overall, 
the model’s Atlantic freshwater flux divergence is generally uncertain, but predicts net 
evaporation of —0.5 ± 0.3 Sv between 30°S and 47°N. 

In the Southern Oceans, freshwater flux divergences (not shown) are even more uncer¬ 
tain due to the low weight of the salinity equations there. The net air-sea flux estimate 
south of 30°S is of precipitation at 0.8±0.9 Sv, which is approximately in balance with the 
net global evaporation between 30°S and 47°N, —1.2 ± 0.5 Sv. 
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Except in few places, freshwater divergences are relatively uncertain. Given the large 
error bars, there is overall agreement with previous studies. Figure4-1 summarizes the 
freshwater divergences (numbers in the boxes). For completeness, those divergences were 
integrated from arbitrary references to compute the so-called “freshwater transports” 
across each section. But any non-divergent field could be added to those transports, 
and only the divergence is of physical interest. Nevertheless, for comparison the integra¬ 
tion of the freshwater divergences was done as follows: 1) The “freshwater flux” through 
the Bering Strait was taken equal to 0.79 ± OSv at 47°N in the Pacific Ocean and to 
—0.95 ± OSv at 47°N in the Atlantic Ocean, that is, following the Wijffels et al. (1992) 
calculation. Through the Pacific-Indonesian Throughflow and the Drake Passage, the 
freshwater flux was set to 0 ± 0 Sv (the uncertainties are taken to zero, because of the 
arbitrariness of adding a non-divergent field). 2) Freshwater divergences were integrated, 
starting, in the Atlantic Ocean, at 30°S; in the Pacific Ocean, at 47°N; and in the Indian 
Ocean, at the northern boundary. Those starting points were chosen to minimize the 
uncertainties on the freshwater transports. In the Southern Ocean, freshwater transports 
were calculated from the southern Atlantic box. In principle, the calculation in the Pacific 
Ocean is a continuation of the integration from the Southern Ocean. But because of the 
large uncertainty in the Southern Ocean, it was preferred to integrate the Pacific Ocean 
freshwater divergences from the North Pacific, resulting in a small shift in the freshwater 
flux at 30°S in the Pacific Ocean. 
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Figure 4-1: Freshwater divergences and transports. While the divergences (positive for 
precipitation/negative for evaporation) are diagnosed from the model, net transports are 
integrated arbitrarily and rely upon stringent assumptions (see text). 


143 

































4.2 Heat transports 


The oceanic energy flux substantially regulates the Earth’s climate. Simulations suggest 
that changes in the thermohaline circulation strength have a large impact on atmospheric 
temperature gradients which control climate (e.g., Manabe and Stouffer, 1988; Marotzke, 
1996 ). Existing estimates from surface observations are highly uncertain (Gleckler and 
Weare, 1997; Josey et ah, 1999). For instance, Josey et al. (1999) found a positive bias 
of 10 to 15 Wm~ 2 over large oceanic regions. Integrated over the Indian Ocean, the bias 
would imply an overestimate of about 0.4Petawatts (1PW= 10 15 W), that is, 100% of the 
Indian Ocean heat gain as estimated, for instance, by Hastenrath and Greishar (1993). 
But even after correction, Josey et al. (1999) found a global heat imbalance of 30 W m~ 2 , 
suggesting large residual uncertainties. 

Analyses from numerical weather prediction centers provide oceanic surface fluxes that 
are often used as a boundary condition to'force ocean models. But meteorological models 
do not provide uncertainties. Using a heuristic approach, Trenberth and Solomon (1994) 
proposed an uncertainty of 30 W m~ 2 over the oceans, based on the European Centre 
for Medium Range Weather Forecasts (ECMWF) outputs. The consequent uncertainties 
in global oceanic heat transports were calculated as 0.3 to 0.5 PW. From a comparison 
between ECMWF, National Meteorological Center (NMC) and net radiative fluxes from 
satellite data, Keith (1995) did a similar calculation and obtained larger uncertainties 
around 0.6 PW for the meridional oceanic energy transport at most latitudes. Recently, 
Gamier et al. (1999) studied the new ECMWF 15-year re-analysis and found a remaining 
global heat imbalance of 8 W m -2 , implying a systematic error over the ocean. 

Arguably, the most accurate energy flux estimations come from circulation estimates 
based on hydrographic sections. Reported flux uncertainties are generally in the range 
±0.3 PW (e.g., Macdonald, 1998; Holfort and Siedler, 1999). The present estimate of the 
global heat transports indicates uncertainties that depend on latitude, ranging between 
±0.1 PW and ±0.5 PW. 
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4.2.1 Global energy transports and heat transfers 

Energy transports are given across each section. However, those are expressed in units 
of Petawatts referred to the Celsius scale for unit consistency with most of the literature. 
Warren (1999) argues, convincingly, that the terminology “heat flux” is inappropriate in 
this context (see also, Bohren and Albrecht, 1998, pp.24-28) and one should refer to the 
“energy flux”, although it does not include the kinetic energy-which is a small contribution 
in this case. 

Because the model privileges mass conservation in the deep layers at the expense 
of the top-to-bottom mass conservation (i.e., Fig. 2-12, page 64), the net flux through 
individual sections is not strictly mass-conservative in semi-closed basins. This noise 
increases towards low latitudes. But in contrast with net fluxes, estimated heat anomaly 
residuals are mass-conserving by definition, and provide a more reliable diagnostic of 
atmospheric heat exchanges. The energy flux divergence in a layer that is in contact with 
the sea surface may be written 

D g = <j) pvC p 9ds + k f OpO r> 0 ds + 

J J lower interface 

where the first integral is made over the lower and side boundaries of the layer and the 
second integral is made over the lower interface, p is in situ density; v is the velocity across 
the boundary; C p is heat capacity; 9 is potential temperature; n is dianeutral diffusion and 
ng is the noise in the heat divergence. The contribution from spurious mass divergences 
is removed from ng , defining the heat anomaly residuals (see also Section 1.2.4). 

A g = D g — 9 <j) pv ds, 

Where 9 is the average potential temperature in the layer. There is no assumption that 
the noise cancels out (i.e., ng ~ 9 § puds) in this diagnostic, and the uncertainty of A# is 
computed explicitly from the solution. In contrast with A g, Dg does not depend upon the 
temperature scale because the mass divergence has been removed. Diagnosing A g instead 
of Dg is an arbitrary choice because the uncertainty in Dg already takes into account the 
spurious mass divergence. However, several simple examples (e.g., AppendixB) suggest 
that A g is closer to the true energy divergence in the presence of spurious mass noise. 
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The uncertainty on A g and on the divergences of all tracers is computed using the full 
covariance matrix of the solution. Because of the existing correlations between errors 
the uncertainties of divergences are occasionally much lower than suggested by a visual 
differentiation of the fluxes. 

The formal uncertainties of energy transports and residuals include the part due to 
oceanic variability in the mass flux. However, the uncertainty due to temporal variations 
in the temperature field is not taken into account. Holfort and Siedler (1999) gave a rough 
estimate of ±0.2 PW for this contribution in the South Atlantic. The authors do not 
provide much detail on the method used to obtain this number, and a more accurate 
estimate would be useful, from XBT data or GCM simulations for instance. 

Figure 4-2 shows the total energy fluxes across latitudes or longitude lines (arrows, 
the white boxes giving the uncertainty). The energy flux is northward everywhere in the 
Atlantic Ocean. 

Heat residuals in each region were computed as the sum of anomaly residuals of each 
layer of the model (outcropping and interior layers). [Because of the large differences in 
the average temperature in surface and deep layers, heat anomalies must be defined using 
the average temperature of individual layers rather than the top-to-bottom value (e.g., see 
the discussion on silica anomalies Section 1.2.4, page 24).] 

There is a net and significant cooling in the North Atlantic and North Pacific be¬ 
tween 24°N and 48°N (respectively — 0.5±0.1PW and — 0.6±0.2PW, Figure4-2). For 
completeness, this figure was also produced with the divergences given in units of W m~ 2 
(Figure 4-3). At low latitudes, warming is observed in both Atlantic and Pacific. In the 
Atlantic, heat residuals in individual basins are uncertain but coherent, so that the to¬ 
tal warming between 30°S and 24°N in the Atlantic is significant, at 0.8±0.2 PW. In the 
Indian Ocean, residuals are noisy everywhere. The Indian Ocean heat budget is tightly 
related to the strength of the Pacific-Indian throughflow, as discussed below. In the central 
Pacific Ocean, there is a net heat gain of 1.75±0.4PW. 

In the Southern Ocean there is a marginal indication of a net cooling in the Indian 
sector and south of Australia. But in the southern Pacific, there is a significant heating 
(0.45±0.2PW) in Layer 3, which outcrops south of 50°S (south of Tasmania). 
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Figure 4-2: Global temperature transports and heat residuals. Each arrow denotes the 
temperature transports between continents. The white boxes correspond to the uncer¬ 
tainties (one standard deviation). Between sections, heat residual anomalies are indicated 
by the full boxes. 
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The map of heat transfers derived from the ECMWF re-analysis (Gamier et al., 1999) 
bears qualitative similarities to our estimates, with most of the oceanic heating occurring 
in eastern and equatorial upwelling regions, while most of the oceanic cooling occurs in 
the western boundary current regions of the northern hemisphere and at high latitudes. In 
the southern hemisphere, cooling was also indicated in the western boundary currents but 
the geographic distribution of the currents is such that there is an equivalent heating in 
the same latitude bands, so that the total heat transfers in each of our boxes are relatively 
small. (Generally, oceanic heating or cooling seems to be correlated with ocean currents 
more than it is with latitude.) In particular, the ECMWF heat transfers indicate oceanic 
heating between 40°S and 65°S, thus consistent with the southern Pacific heating that is 
calculated here. The Comprehensive Ocean-Atmosphere Data Set (COADS, daSilva et 
al., 1996) also indicates global heating around 50°S, but a strong cooling between 43°S 
and 15°S as well. (Robbins et al. (personal communication, 1999) suggest that COADS 
estimates is erroneous in this region because of the lack of data.) 

In summary, given the large range of existing energy flux estimates and uncertainties, 
we find no qualitative incompatibilities with our estimate. We now turn to the quantitative 
comparison. 

Atlantic energy fluxes 

Figure 4-4 compares our estimate of the Atlantic energy fluxes (black stars) with other 
published estimates. The Trenberth and Solomon (1994) uncertainties were derived heuris- 
tically, assuming that spurious heat residuals on land were representative of the energy 
flux uncertainties everywhere. In this view, they provide a lower bound on the actual 
uncertainty of the ECMWF model. The Semtner and Chervin estimate corresponds to 
the “POCM-4B:” integration (Semtner and Chervin, 1992; Stammer et al., 1996; Jayne, 
1999). 

There are no incompatibilities, within one standard deviation, with the previous hydro- 
graphic estimates. The present net energy flux across 7.5°N is probably underestimated 
because of the high variability in the western part of this section (e.g., Johns et al., 
1998). Heat residuals, because they are corrected for the spurious mass imbalances that 
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Figure 4-4: Atlantic energy transports (black stars with thick error bars) and residuals 
between selected latitudes (numbers above the northward-pointing arrows, positive for 
heating). (The energy flux at 7.5°N is probably noisy due to variability in the North Brazil 
Current region.) Errorbars from Holfort and Siedler (1999) are of ±0.3 PW; Rintoul and 
Wunsch (1991), ±0.2 PW; Hall and Bryden (1982), ±3PW; and Saunders and King (1995), 
±0.1 PW. (Adapted and completed from Macdonald and Wunsch (1996) and Heras and 
Schlitzer (1999).) 
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are present in the surface layers, indicate an uncertain ocean warming between 4.5°S and 
7.5°N (0.35 ± 0.5 PW). 

Heat residuals integrated over large latitude bands are indicated on Figure 4-4 by the 
numbers on the blue and red arrows. While residuals in individual boxes are uncertain 
and noisy at low latitudes, the total warming over the central Atlantic (30°S-24°N) is 
significant, with 0.7±0.2 PW warming. Poleward of this latitude band, net cooling is 
observed, with — 0.3±0.15 PW in the South Atlantic (All to 30°S); —0.5 ±0.1 PW in the 
North Atlantic (24°N-48°N) and —0.6 ±0.1 PW in the northern North Atlantic (north of 
48°N). 

Except for the Keith (1995) estimate in the North Atlantic, most other estimates from 
GCM, COADS and National Weather Prediction centers are lower than ours in the At¬ 
lantic Ocean. In the North Atlantic, there is a fair agreement about the heat loss between 
24°N and 48°N (except for Hu, 1997 and Semtner and Chervin). The warming inferred 
in the central Atlantic is compatible with most estimates, except Hu (1997) and Hasten- 
rath (1982), which have a lower heat gain. Between All (plotted at 45°S on Fig. 4-4) and 
30°S, no climatology indicates the same cooling that I found. But because section All is 
not entirely zonal (Fig.3-9), the discrepancy may be due to the differences in area of heat 
flux integration. Saunders and King (1995) obtained a slightly lower energy transport 
across All: 0.5 PW. 

Indian Ocean heat budget and Pacific-Indonesian Throughflow 

As a consequence of the large inflow of warm water from the PIT (1.36 ± 0.15 PW energy 
transport across JADE), the net estimated heat gain over the Indian Ocean is indis¬ 
tinguishable from zero (0.1 ± 0.2 PW). The heat residual anomalies in individual boxes 
(Fig. 4-2) are also not significantly different from zero. Only in the subtropical box (de¬ 
tail not shown) is there a marginal indication of oceanic heat loss in the surface layer 
(0.5±0.4PW) and heat gain in the second layer (0.4±0.6PW) which outcrops almost 
completely by 32° S. 

The ECMWF energy flux estimates of Siefridt et al. (1999) and Gamier et al. (1999) 
effectively show no significant heating over the northern part of the Indian basin; neither 
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does the NCEP/NCAR re-analysis of Kalnay et al. (1996, shown in Gamier et al. 1999). 
However, both indicate a cooling of 0.4 PW between 20°S and 32°S. On the other hand, 
the COADS data of da Silva et al. (1996) and the climatology of Oberhuber (1988) and Hsi- 
ung(1985) indicate a warming of about 1.2 PW in the Indian Ocean north of 20°S. The 
GCM assimilation computations of Lee and Marotzke (1997, 1998; hereafter LM97,98) 
which are forced with the Oberhuber (1988) climatology, indicate a consistent heating be¬ 
tween 20°N and 14°S (0.8 PW) and a cooling in the subtropical region (14°S-32°S, 0.3 PW). 
Several other GCM studies based on different surface forcing show heating magnitudes 
north of 14°S that are similar to those of LM97,98 (i.e., their Fig. 19). All GCM-based 
inversions have a relatively weak or non-existent PIT, so that local atmospheric heating 
accounts for most of the temperature flux leaving the basin at 32° S. 

The estimate of heat divergence over the Indian Ocean is thus largely determined by the 
size of the PIT (e.g., Godfrey, 1996 ). Based on a PIT varying randomly from 1 to 11 Sv, 
Robbins and Toole (1997) estimated the energy flux divergence for the Indian Ocean north 
of 32°S at 0.4±0.2 PW, similar to the result of Hastenrath and Greischar (1993) which was 
based on a climatological combination of ship observation, radiation measurements and 
sub-surface temperature. 

My lower heat gain over the Indian basin is the consequence of the relatively large PIT 
temperature flux combined with a relatively weak overturning circulation. Similarly, the 
energy transport in the South Pacific Ocean is oriented northward because of the large 
PIT (Fig. 4-2). Most other estimates obtained a weak (and 100% uncertain) southward 
transport (see Tsimplis et al. (1998), Fig.17, for a review). Tsimplis et al. (1998) actually 
obtained northward energy transport at 32°S assuming a 10-Sv PIT. 

Indo-Pacific energy flux 

To circumvent the uncertainty related to the PIT energy transport, the net meridional 
temperature transport is traditionally summed over the Pacific and Indian Oceans. Fig¬ 
ure 4-5 shows the Indo-Pacific energy fluxes. I find a net heat gain of 1.6 ±0.4 PW between 
30°S and 24°N, with a —0.6 ± 0.2 PW heat loss (largely in the North Pacific) between 
24°N and 48°N. Between 30°S and 20°S, the model has a marginal indication of cooling, 
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Figure 4-5: Indo-Pacific energy transports and residuals. See Fig. 4-4 for the legend. 


—0.5 ± 0.5 PW. 

Those estimates of the net energy transports are in agreement with Macdonald and 
Wunsch (1996) and Bryden et al. (1991) in the North Pacific within one standard deviation. 
In the south Indo-Pacific, our southward energy flux is smaller than the Macdonald and 
Wunsch (1996) estimate, but still compatible within two standard deviations. 

The heat divergences obtained by Garnier et al. (1999) and Semtner and Chervin are 
similar to ours; but all other estimates suggest a larger heat gain of up to 3.5 PW in the 
tropical band (i.e., Hastenrath, 1982), which are incompatible with our estimates. 

Global energy fluxes 

Figure 4-6 shows the global energy transports and residuals. Globally, the ocean heated 
(2.3±0.4 PW) between 30°S and 24°N and cooled poleward of those latitudes, with —1.7± 
0.2 PW in the northern hemisphere and —0.7 ± 0.3 PW in the southern oceans. The 
Macdonald and Wunsch (1996) estimate indicated a similar warming, as did the Hu (1997) 
isopycnal model. Most other estimates from numerical models and in situ data indicate a 
larger warming of around 3.5 PW in the tropical band. But the large uncertainties in the 
ECMWF and NMC model products provided by Keith (1995, green shaded area) indicate 


153 





60S 40S 20S O 20N 40N 60N 

latitude 


Figure 4-6: Same as Fig.4-4, but for the global energy transports and residuals. The green 
band indicates the Keith (1995) uncertainties. 


154 





that there are no incompatibilities with our estimates. Gamier et al. (1999), from the 
same ECMWF model but using the 15-year re-analysis, obtained a much lower oceanic 
warming at low latitudes than Keith (about 2.5PW). The Semtner and Chervin estimate 
is also smaller, with about 1.75 PW warming between 30°S and 24°N. 

Conclusions for energy fluxes 

Apart from hydrographic estimates, the existing estimates of oceanic energy transports 
and ocean-atmosphere energy fluxes span a large range of conflicting values. The degree 
of agreement between previous estimates and my model varies with location. Generally 
speaking, the ECMWF 15-year re-analysis (Gamier et al., 1999) is in best agreement with 
my energy flux estimates over most oceanic regions. Many past estimates were provided 
with no uncertainties, but are incompatible, within my model’s uncertainties, with the 
present estimates. However, according to the uncertainty estimate of Keith (1995), all 
meteorologically-based energy flux estimates are consistent with ours within one stan¬ 
dard deviation over the northern hemisphere and within two standard deviations over the 
southern hemisphere. Furthermore, according to Keith (1995), this uncertainty is a lower 
bound on the actual uncertainty from NWP centers. 

On the other hand, our estimate is compatible within error bars with the published 
estimates from hydrographic sections. Nevertheless, the cooling that we obtain in the 
Southern Ocean is weaker than most other estimates, resulting in a smaller poleward 
energy flux at 30°S. 
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4.3 Export production and nutrient cycle 

Quantification of the global biogenic production in the ocean is one of the major issues 
in climate prediction. Photosynthesis associated with primary production in the surface 
layers uses carbon dioxide to form unicellular algae. Carbon is thus fixed, forming or¬ 
ganic matter in the surface layers. A large part of the organic matter produced in the 
surface waters, (also called the total primary production) is recycled in the euphotic zone 
(regenerated production), with the balance exported towards the deep ocean (export pro¬ 
duction). This export represents a flux of carbon towards the deep ocean which is believed 
to control and regulate air-sea CO 2 exchanges with a direct implication on climate. Few 
methods exist to estimate this “export production”. Satellite measurements of ocean color 
provide an estimate of the biomass in upper layers (e.g., Antoine et al., 1996), from which 
the primary production can be derived. Export production is related to the primary pro¬ 
duction. (The difference between the export and primary production is the “regenerated 
production”.) But not only is the satellite-derived primary production uncertain, but also 
its relationship to export production is poorly known, and probably changes with time 
and location. Export production can apparently reach up to 50-80% of the primary pro¬ 
duction during episodic blooms or in high productivity areas; otherwise, it varies in the 
range 5-10% (Buesseler, 1998). 

The existing measurements of export production from sediment traps vary by one or 
two orders of magnitude, depending upon region and regime. Lampitt and Antia (1997) 
reviewed the existing measurements normalized at 2000 m. Using their normalization 
function, the annual averages exported production at 150 m is in the range 0.3 to 3molC 
m 2 yr _1 in all oceanic regions except polar oceans (there is at least 100% uncertainty 
from this conversion, K. Buesseler, personal communication, 1999). In polar oceans, the 
variations are more extreme, with values reported between 0.08 and 5molC m _2 yr _1 . 
But different biases can affect sediment trap results (Siegel and Deuser, 1997), leading to 
underestimation of the true particle flux by a factor of 3 to 10 (Buesseler, 1991). Export 
production estimates based upon the 234 Th method (Michaels et al., 1994; Buesseler, 1998) 
are about ten times higher, ranging from 1 to 10 mol C m _2 yr“ 1 for all oceans excluding 
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the polar regions, and from 0.3 to 30molCm~ 2 yr _1 in the polar oceans. 

On the other hand, export production is directly related to dissolved nutrient fluxes. 
A flux of particles from the surface to the deep layers creates a sink of dissolved nutrients 
in the upper layers by consumption and a source in the lower layers through dissolution 
(also called remineralization). In the hypothetical steady-state, those divergences are 
balanced by horizontal and vertical transports. Similar divergences may occur as a result 
of horizontal advection of dissolved organic matter. 

In this section, flux divergences of dissolved nutrients are estimated and discussed in 
terms of particle fluxes and biogeochemical cycles. Redfield ratios of dissolved nutrients are 
used to convert the results between the different dissolved properties to compare my results 
with the literature (e.g., Redfield et al., 1963). I use the Anderson and Sarmiento (1994) 
stoichiometric relation P/N/C/—O 2 = 1/16/117/170. The discussion is started with the 
silica cycle, where I justify the assumption of quasi-conservation of dissolved silica from 
surface to bottom over the water column. Following that, nitrate and phosphate flux diver¬ 
gences are analyzed. Oxygen flux divergences, which are tightly related to the biological 
cycle, are discussed at the end. 

4.3.1 Silica 

Silica is a major nutrient that feeds biological production in the oceans. Siliceous plank¬ 
tonic organisms, mainly diatoms and radiolaria, consume dissolved silica to build their 
skeletons (e.g., Broecker and Peng, 1982 ). In the opaline form, silica goes through the 
food chain and ends up in sinking fecal pellets and skeletons. Most of the opal is dissolved 
in deep waters, although a fraction of the particles are buried on the sea floor. Opal burial 
occurs mostly in highly productive areas, forming remarkable opaline bands on the sea 
floor under upwelling areas (Lisitzin, 1972; Aston et al., 1983; Archer, 1993). 

The dissolution process within the water column is poorly known. Edmond (1974) 
suggested that most dissolution takes place on the sea floor. Using sediment traps near 
Guinea, Takahashi and Honjo(1981) also found that radiolarian shells did not dissolve 
substantially between about 1000 m and the sea floor. Sayles et al. (1996) reached a similar 
conclusion in the Sargasso Sea, where they found that below the first 150 m, where 66% 
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of the dissolution had taken place, no apparent dissolution occurred before particles reach 
the sea floor. Particles particles reaching the bottom remain on the sea floor for a few 
hundred years; most eventually remineralize with only a small fraction fully buried and 
lost to the ocean system. River runoffs are believed to provide an equivalent gain of silica 
to the ocean (Treguer et al., 1995) ). 

These results are uncertain. Not only are they localized to a few places, but also sedi¬ 
ment trap measurements can be biased by a factor of 10 (e.g., Michaels et al., 1994). Such 
studies are complicated by both the diversity of siliceous species and the vertical nature 
of the data (i.e. along the wire of the mooring) in a situation where horizontal advection 
and inhomogeneous surface production may bias the measurement (e.g., Buesseler, 1991; 
Siegel and Deuser, 1997). 

Treguer et al. (1995) and Nelson et al. (1995) established a global budget for the silica 
in the oceans. Their estimate for the total silica input to the oceans, mainly through river 
runoff and erosion rate, was close to 200 kmol Si s _1 . They estimated the upper bound for 
the global primary production of siliceous organisms at 9000 kmol Si s _1 , using existing 
14 C production estimates combined with silica/carbon ratio in siliceous organisms. Their 
lower bound estimate, 6400 kmol Si s _1 , was obtained from a model. 

Should dissolved silica be conserved over the water column? 

As mentioned in Chapter 1, quasi conservation of dissolved silica over the water column 
is required as a constraint on the physical circulation. The rational for this constraint is 
given here. 

In the upper 100 m (the euphotic zone), dissolved silica is almost totally depleted 
by diatoms and radiolarians. According to Nelson et al. (1995), 50% of the global opal 
production is recycled in the euphotic zone, allowing a maximum global particle flux of 
4500 kmol Si s _1 to enter the deep layers (those numbers are highly uncertain). As noted 
above, some of these particles are buried on the sea floor. The total burial/production 
ratio is, on average, 3% and varies strongly with location. The burial is believed to be 
almost negligible in the large oligotrophic regions (e.g., Treguer et al. 1995) . Similarly, in 
high-productivity regions such as the Southern Oceans, recent isotopic measurements from 
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the JGOFS Antarctic program indicate that little burial has occurred for at least the last 
5000 years (a maximum of 5% of the total production, F. Sayles, personal communication, 
1999). This result contradicts previous beliefs that up to 25% of the total production was 
buried there (DeMaster, 1991). 

Globally, the net sources and sinks of silica thus appear to be much smaller than 
the biogenic fluxes. The global sink of silica estimated from sediment data is about 
220 kmol Sis -1 , which is close to the total input from river runoff, at 160 kmol Si s -1 (Treguer 
et al., 1995). Other sources of silica, such as hydrothermal activity, basaltic erosion and 
eolian inputs, are probably even smaller (40 kmol Si s -1 in total, Treguer et al., 1995) 
although little is known about them. Regionally, those sources and sinks may be more 
important, but still smaller than biological fluxes. 

Two additional processes may lead to non-conservation of dissolved silica between 
hydrographic sections. Namely, the erosion of the sea floor during a transient period, or 
the horizontal advective transport of siliceous particles across an oceanic section. 

Is sea floor erosion able to provide a large flux ? In the steady state, it is believed 
that silica is buried in the sediments, i.e., there is a net opal flux to sediments. The 
burial can lead to a sediment accumulation up to 17 cm/Kyr in the Southern Oceans 
(DeMaster, 1981). In my view, a net erosion would have to correspond to the reversal 
of the accumulation during a hypothetical episodic (on climate time scales) event. A 
100 kmol Si s -1 source over an area such as the North Atlantic between 24°N and 36°N 
corresponds to an area average erosion rate of 0.4 mol Si yr -1 m -2 . In this case, taking 
a sediment density of 0.4 gem -3 (DeMaster, 1981), with a 1.5% opal content (Sayles et 
al., 1996, p.20), the sea floor erosion rate would be 


0.4mol Si yr 2 m 2 x 60 — —. x 

mol bi 


1 g Tot 
0.015 g Si 


1 cm 3 
0.4 g Tot 


= 4m /Kyr. 


An erosion rate of 4 m /Kyr (meter per Kilo-year) seems unreasonably large. Besides, 
direct measurements of the benthic silica flux with benthic chambers indicates maximum 
fluxes of 0.036 mol Si yr -1 m -2 near Bermuda (Sayles et al., 1996). Their calculation based 
on the Si gradient in the pore water samples near the sediment-water interface led to the 
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same result. In the Equatorial Pacific, similar measurements indicated larger fluxes (up 
to 0.25 mol SiC >2 yr _1 m~ 2 , McManus et al., 1995). The largest dissolution rates are found 
near Antarctica, with up to 1.5 mol SiC >2 yr _1 m~ 2 (F. Sayles, personal communication, 
1999). Still, measurements indicate that those fluxes from the sea floor derive from the 
opal rain, thus implying net deposition and no net erosion to the sea floor. As a result, 
the net source of silica near the sea floor should be approximately in balance with the opal 
rain. 

Can siliceous particles be advected across sections ? Horizontal advection of siliceous 
particles may create an apparent source or sink in the dissolved silica budget. Values of 
vertically integrated rates of biogenic silica production vary between 0.2 mol Si yr -1 m~ 2 
(in the Sargasso Sea) and 18 mol Si yr -1 m -2 (in coastal upwelling regions; for a summary, 
see Nelson et al., 1995). 

If the particles rained with a downward velocity of 100 m/day (Banse, 1990 ), the 
average particulate silica concentration over the water column would range from 5 x 
10~ 3 pmol/kg to 500 x 10~~ 3 /imol/kg. Multiplied, for instance, by the Florida Current 
flow (30 Sv), the silica flux would range 0.15 to 15 kmol Si s -1 . Taking into account that 
50% of the particles are permanently recycled in the upper layer would increase the par¬ 
ticle concentration in the fast-moving layers-and thus the flux by a few fold. Honjo and 
Manganini (1993) argued for a slower downward particle velocity of 50 m/day, which, for 
the same production, would double the particle concentration and the resulting horizontal 
flux. 

The particle flux across a whole section is difficult to estimate, but it is likely to be 
the same order of magnitude as the above estimated Florida Current flux because the net 
flux is likely dominated by regions of simultaneous high productivity and strong currents. 
Particle advection thus appears to be, within our knowledge, smaller than our reference 
uncertainty of 100kmol Si s -1 (i.e., Section 2.2.4), although the magnitude is not so weak 
that it can be completely ignored. Note that within the steady state assumption, strong 
particle advection would be associated with a large silica depletion on one side of the 
section and an enrichment on the other side. 

Seasonal effect may also be an issue. Williams and Follows (1998) found that important 
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seasonal correlation between nitrate and Ekman transports at the edge of the subtropi¬ 
cal gyre may produce “eddy” fluxes that are not captured by analyses based on annual 
averages. But for nitrate this effect is weak in comparison to the magnitude of our uncer¬ 
tainties (see the discussion below, page 168), and it is anticipated to be small for the silica 
too. 

To first order therefore, sea floor erosion and particle advection do not appear to create 
large sources or sinks of dissolved silica. Other sources or sinks are a priori relatively 
small. Within present knowledge, dissolved silica can thus be taken as conserved within 
ilOOkmol Si s -1 over large oceanic areas such as the North Atlantic between 24°N and 
36° N. (Practically, uncertainties are larger due to oceanic variability and measurement 
noise, Section 2.6.) (For practical reasons, no constraint was imposed to the global silica 
budget.) 

Silica flux divergences 

Silica, and other nutrient flux divergences were examined in each box. Most of the resid¬ 
uals in individual layers have a 100% uncertainty. However, when a sign change was 
noticed between the surface and the deep layers, nutrient residuals were summed sepa¬ 
rately surface/deep. In boxes where residuals within individual layers were more than 
150% uncertain, the summation was done over all layers (top-to-bottom). Figure 4-7 
shows those silica fluxes and divergences. (Note that for technical reasons, and in contrast 
with all other tracers, silica divergences are not corrected for mass divergences, i.e., the 
anomaly formulation was not used (Section 1.2.4, page 24). As a result, some of the silica 
divergences may be associated with mass imbalances). 

Silica fluxes (arrows) are directed southward within the Atlantic Ocean, and poleward 
elsewhere. In the Atlantic, the flux is of about -150 ± 75kmol Si s -1 . In other regions, 
meridional fluxes are more uncertain. Silica residuals are indistinguishable from zero 
within one standard deviation in most places, i.e., the conservation constraint was met. 
But there are exceptions in the central Pacific and subtropical and tropical Indian Oceans 
where the residuals are between one and two standard deviations. In the central Pacific 
Ocean, there is an apparent source of dissolved silica (+0.35±0.15 mol Si02 yr _1 m~ 2 ) while 
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Figure 4-7: Global dissolved silica transports and residuals. The length of each arrow 
corresponds to the silica transports between continents. The white boxes behind each 
arrow indicate the uncertainty (one standard deviation). Between sections, silica diver¬ 
gences are indicated by the full boxes, either top-to-bottom (single box) or surface/deep 
(double box). In the later case, the number indicates the interface (in neutral density 
units). Residuals are not given where the uncertainty exceeds 150%. 
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a sink is found in the Indian Ocean (—0.28 ± 0.24mol Si02 yr _1 m~ 2 ). This imbalance 
may be a consequence of the high Pacific-Indonesian Throughflow. (The model’s PIT 
transports 400 ± 300 kmol Si s -1 to the west, implying a source in the Pacific and a sink in 
the Indian Ocean that may be erroneous-although both are indistinguishable from zero.) 

The Southern Ocean is an uncertain net sink of silica at —480 ± 380 kmol Si s -1 
(equivalent to —0.18 ±0.15 mol Si02 yr -1 m -2 overall). This rate is greater than the 
largest burial rates measured in the high production regions by the JGOFS Antarctica 
program: 0.05 mol Si02 yr -1 m -2 (F. Sayles, personal communication, 1999). 

In most of the temperate regions, there is a depletion of silica in the upper layers 
ranging from —0.1 ± 0.05 to —1 ± 0.6mol SiC >2 yr -1 m~ 2 , with a corresponding source in 
the underlying deep layers. Although highly uncertain, this residual structure corresponds 
to what one expects from silica consumption and opal dissolution. 

The global average production rate of silica was estimated to be 0.6 mol SiC >2 yr“ 1 m -2 
by Nelson et al. (1995), with variations between 0.2 and 8.3 mol SiC >2 yr _1 m~ 2 , the highest 
production rates occurring in coastal upwelling areas. Between 10% and 50% (e.g., Honjo 
and Manganini, 1993; Nelson et al., 1995) of this is exported to the deep layers, resulting in 
an opal dissolution rate ranging 0.02 to 4 mol SiC >2 yr -1 m~ 2 . Lampitt and Antia (1997) re¬ 
ported opal rain fluxes from particle traps in the same range (0.1 to 1.2 mol SiC >2 yr _1 m~ 2 , 
excluding polar oceans, where low export production attains 0.01 mol SiC >2 yr -1 m~ 2 ). The 
magnitude of our large-scale silica sources in the deep layers in temperate regions is thus 
within the previously measured range. 

4.3.2 Nitrate and Phosphate 

Figs. 4-8 and 4-9 show the dissolved nitrate and phosphate net fluxes and residuals. As for 
silica, residuals of dissolved inorganic nutrients were integrated, either top-to-bottom (sin¬ 
gle red or blue box) or within surface and deep layers if there was a coherent sign"change. 
In contrast with silica, nitrate and phosphate residuals are estimated from anomaly equa¬ 
tions, i.e., corrected for mass imbalances. Net fluxes of nitrate and phosphate (arrows) are 
uncertain overall. There is an indication of southward nitrate flux from the Atlantic Ocean 
at 45°S (All) (i.e., a nitrate divergence of 340 ± 120 kmol Ns -1 between 48°N and 45°S) 
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as well as phosphate (16 ± 8kmolPs _1 ). 70% of the corresponding nutrient regeneration 
occurs in the South Atlantic between 30°S and 45°S. 

In the combined tropical Atlantic and “Rio” (19°S-30°S) boxes, residuals indicate a 
depletion of nutrients in the surface layers and a source in the deep layers (Fig. 4-8). 

This same feature is present in the Central and North Pacific. All residuals are highly 
uncertain but this pattern interestingly corresponds to what would be expected from a 
hypothetical particle transport that consumes nutrients at the surface and restores them 
in the deep layers through remineralization. In the rest of the Ocean, the residuals were 
too incoherent and too uncertain to make this distinction. Net sources or sinks of ni¬ 
trate and phosphate are found in several oceanic regions, noticeably in the North Atlantic 
(0.1 ± 0.1 molN m _2 yr“ 1 ). Such a source was also found by Rintoul and Wunsch (1991) 
using older data, between 24°N and 36°N (versus 24°N and 48°N here). Rintoul and Wun¬ 
sch (1991) found that the circulation derived from the older data could not be made com¬ 
patible with nitrate conservation. The apparent source led them to suggest the advection 
of significant dissolved organic matter into the region, where it subsequently remineralized. 

Nitrogen fixation is another possible explanation for the apparent sources. However, 
phosphate (Fig. 4-9) exhibits a similar source in the North Atlantic (0.02±0.007 molP m~ 2 yr _1 ), 
and no equivalent to fixation is possible. In general, phosphate residuals are similar to 
nitrate residuals (Figs. 4-8 and 4-9). A few exceptions exist (e.g., south of Africa), but 
the differences are too uncertain to draw conclusions implying different processes. Con¬ 
sequently, the nutrient divergences that I find are likely to be due to either advection of 
dissolved organic matter, or to an alias of the seasonal cycle in the hydrographic sections. 

The nitrate and phosphate residuals were integrated over latitude bands (Table 4.4). 

Both exhibit similar features, with about 100% uncertainty. High latitude oceans expe¬ 
rience nitrate and phosphate sources while the tropical band is a sink. These inferred 
sources and sinks may be the result of advection of dissolved organic matter, or temporal 
aliasing of the nutrient fluxes across the sections, but all are indistinguishable from zero. 
Dissolved nitrate and phosphate are both balanced globally. 

Net sources of dissolved nutrients in the deep layers range from 0.5±0.24 mol N m _2 yr _1 
(tropical North Atlantic) to 0.7 ± 0.2molNm~ 2 yr -1 (South Atlantic). Moored sediment 
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Figure 4-8: Same as Fig. 4-7, but for dissolved nitrate. 
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Figure 4-9: Same as Fig. 4-7, but for dissolved phosphate. 

trap measurements of particle export from various oceanic regions are of the same or¬ 
der, ranging from 0.04 to 0.4 molN m~ 2 yr _1 (Lampitt and Antias, 1997). Isotopic mea- 
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Latitude 

Nitrate (kmolNs 

Phosphate (kmolPs 

24°N-48°N 

100±110 

10±8 

30°S-24°N 

-90±200 

-13±13 

south of 30° S 

100±250 

10±15 

Global total 

40T300 

-3±20 


Table 4.4: Nitrate and phosphate residuals integrated globally over latitude bands. 

surements based on 234 Th suggest that sediment traps on sediment traps may under¬ 
estimate the actual particle flux by a factor of up to 10 (e.g., Michaels et al., 1994; 
Buesseler, 1998). Export production estimates from isotopic measurements in temperate 
and tropical oceanic regions range from 0.05-1.4 molNm~ 2 yr _1 . New production 1 esti¬ 
mates are in the same range, with about 0.5 molN m _2 yr _1 (e.g., Jenkins and Goldman, 
1985), which is in rough agreement with the nutrient sources and sinks found in the sur¬ 
face layers in the Atlantic, namely, —0.4 ± 0.3 molNm _2 yr _1 in the South Atlantic and 
—0.5 ± 0.25 mol Nm~ 2 yr _1 in the tropical North Atlantic. 

The consistency between local measurements and large scale averages may be fortu¬ 
itous. For instance, a single layer may outcrop so that both consumption and remineraliza¬ 
tion occur in the same box with the two regimes cancelling. Also, neutral surfaces do not 
correspond to the bottom of the sun-lit layer (or 150-m, the level to which sediment trap 
measurements were normalized). Nevertheless, the magnitudes of the nitrate residuals are 
in close agreement with independent measurements. Given the large uncertainties, the 
inverse model estimate provides at best an upper bound for the average remineralization 
process. Overall, nitrate residuals (including deep/surface and net in a box) are in the 
range — 0.8 ± 1 molNm -2 yr -1 to —1.8 ± 1.8 molNm -2 yr -1 . 

4.3.3 Discussion of nutrient fluxes and divergences 

Despite the large uncertainties, the magnitude of the maximum export production agrees 
closely with independent measurements from moored sediment traps and isotopic data. A 
(possible) signature of time-mean, large-scale average export production is found for silica, 

'New production is defined here as the amount of production deriving from the nutrient supply from 
the deep layers in a hypothetical one-dimensional balance. In the absence of horizontal advection, new 
production is equivalent to export production. 
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nitrate and phosphate in several temperate regions, with the exception of the Indian 
Ocean. Net sources and sinks of nutrients are found in several regions too, but all are 
indistinguishable from zero within a two-standard-deviation uncertainty. Slight differences 
between nitrate and phosphate flux divergences were noted. These are, in general, not 
significant, thus providing no evidence for nitrogen fixation or denitrification. In the Indian 
Ocean, there is an (uncertain) source of nutrients in the north and an (uncertain) sink in 
the tropical and subtropical regions. It is not clear whether this pattern is an artifact of an 
exceptionally large Pacific Indonesian Throughflow on JADE 1989 that created episodic 
nutrient imbalances, or a consequence of particle advection across the sections. 

Imbalances may also relate to alias of the seasonal or interannual variability in the 
biological activity (particle flux) and/or the advective nutrient fluxes. At the edge of the 
subtropical gyre, Williams and Follows (1998) suggested that seasonal correlations between 
nutrient concentrations and Ekman transports could be important. But their biases in the 
nitrate fluxes are of order 4kmolNs _1 (McLaren, personal communication, 1999), which 
is 10 times smaller than my uncertainties. 

Time variability in biological activity can also introduce additional biases when esti¬ 
mating a time-average because of the different time scales between production and disso¬ 
lution. For instance, Deuser et al. (1995) found a 50% change in the opal export flux to 
the deep Sargasso Sea between 1978 and 1991. Such a change should be paralleled in the 
consumption rate, but, if the opal residence time is of order 100 year (Sayles et ah, 1996), 
the change would not be seen in the dissolution on the sea floor. 
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4.4 Oxygen fluxes and divergences 

4.4.1 The oxygen cycle 

The oceanic oxygen cycle is closely related to production and to the CO 2 cycle through 
photosynthesis and remineralization. Photosynthesis produces oxygen and consumes CO 2 , 
while remineralization consumes oxygen. But in contrast with nutrients, the oceanic 
oxygen cycles include important air-sea fluxes. The amount of oxygen that is pumped 
or outgassed by the ocean depends upon the partial pressure of oxygen which, in the 
ocean, is controlled by (1) temperature and, to a minor extent, salinity (e.g., Garcia and 
Gordon, 1992); (2) biological production and respiration; and, locally (3) oceanic transport 
and upwelling. Control number (1) is called the “solubility pump” (e.g., Stephens et ah, 
1998): when water is warmed, partial pressure is increased and oxygen (and CO 2 ) are 
released to the atmosphere. The opposite happens when water is cooled. Control (2) is 
called the “biological pump”, as it consumes oxygen at depth and produces oxygen at the 
surface. 

Depending upon the oceanic region, the solubility pump, the biological pump and 
the oceanic transports have either reinforcing or opposing effects (e.g., Keeling and Peng, 
1995). For instance, in equatorial upwelling regions, where productivity is limited by 
other factors than nutrient supply, surface warming suggests (provided that the surface 
oxygen concentration is saturated) an oxygen efflux (ocean->atmosphere) while upwelling 
of low-oxygen, cold deep waters creates an influx. On the other hand, in subtropical gyres, 
warming, productivity and advection create an efflux while at high latitudes cooling and 
convection both imply an influx. 

While several studies address the seasonal cycle of oceanic oxygen, relatively few pub¬ 
lications address the average flux on longer time scales. Stephens et al. (1998) compared 
the time-average oxygen and CO 2 concentration and fluxes using ocean biological general 
circulation models, an atmospheric model and atmospheric observations from 9 stations. 
They suggested that the ocean GCMs (combined with an atmospheric tracer model) un¬ 
derestimated the southward transports of oxygen and CO 2 , yielding a large gain (air—^sea) 
of oxygen in the Southern Ocean which they suggested was not present in the observations. 
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As presented below, my analysis indeed suggests important oxygen air—»sea flux in the 
Southern Ocean. 

4.4.2 Oxygen fluxes and divergences 

Figure 4-10 shows the net fluxes and residuals of dissolved oxygen from my standard so¬ 
lution. In contrast with the nutrients, oxygen fluxes and divergences are significantly 
different from zero in many regions. There is a large influx (atmosphere to ocean) of 
oxygen in polar and temperate regions and a loss in tropical regions. As a result, oxygen 
fluxes are directed equatorward almost everywhere (except in the Indian Ocean-because 
of the oxygen transport associated with the Indonesian Throughflow). 

Analysis of the flux divergences (also called influx or efflux , in contrast with flux which 
is used for horizontal transport across sections) show that in the 24°N-48°N latitude band 
(globally) there is an oxygen influx of 1050 ± 600kmol 0 2 s -1 . North of 48°N, there is 
an additional influx of 2400 ± 650kmol 0 2 s -1 . At lower latitudes (24°N-30°S), oxygen is 
consumed or/and restored to the atmosphere at a rate of —7200±800kmol 0 2 s -1 , while in 
the Southern Ocean there is an influx of the same magnitude as the northern counterpart, 
3400 ± 700kmol 0 2 s -1 . Figure4-11 shows those sources and sinks (blue columns). The 
global oxygen budget is in balance within the uncertainty (-300 ± 800kmol 0 2 s -1 ). 

Locally, net loss of oxygen reaches values of up to 8 ± 1 mol 0 2 m -2 yr _1 in the South 
Atlantic. In this particular region the oxygen loss occurs over the first 1000 dbars, con¬ 
sistent with remineralization. (Note that if particles are created and remineralized in the 
same oceanic region, oxygen consumption in the subsurface layers would correspond, in 
the steady state, to an equivalent oxygen production and possibly outgassing in the sur¬ 
face layers through photosynthesis. By conservation of oxygen, a net oxygen utilization 
observed in a given oceanic region, even at depth, ultimately corresponds to outgassing.) 

In the North Atlantic, South Atlantic and North Pacific temperate regions, distinctive 
regimes can be identified in surface and deep layers, with an oxygen gain in the surface lay¬ 
ers, associated with atmospheric influx and, possibly, a contribution from photosynthesis. 
In deeper layers, an uncertain loss is found, probably associated with remineralization. In 
layers with oxygen gain, there is no correlation with nutrient loss and the biological pump 


170 




Figure 4-10: Dissolved oxygen fluxes and flux divergences. The caption is the same caption 
as for nutrients (Fig. 4-7), i.e., residuals were summed individually in deep/shallow layers 
when the distinction was possible, top-to-bottom otherwise. 
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Figure 4-11: Global oxygen sources and sinks for selected latitude ranges. The blue 
bars show the inverse model estimate, while the red bars correspond to the hypothet¬ 
ical contribution from the solubility pump, estimated from the energy flux divergences 
(-1500kmol 0 2 s-VPW, Holfort et al., 1999). 

is probably masked by the solubility pump. In deeper layers (e.g., in the South Atlantic, 
19°S-30°S), oxygen utilization is correlated with nutrient input. 

Table 4.5 shows a list of previous estimates of oxygen utilization rate (OUR) and pro¬ 
duction from different methods (adapted from Sarmiento et ah, 1990). The OUR estimate 
from Sarmiento et al. (1990) is 2.4-2.5mol02 m -2 yr -1 , when derived from tritium in¬ 
ventories, and 8.5±0.8molO2 m~ 2 yr _1 when derived from 228 Ra measurements. Those 
authors argue that the 228 Ra-values are more reliable than those from the tritium inven¬ 
tory. The magnitude of oxygen loss that I obtained in the South Atlantic (20°S-30°S) is of 
the same size (~ 8mol02 m~ 2 yr -1 ), although, again, I did not attempt to separate OUR 
from atmospheric exchanges in the surface layers. In the tropical Atlantic, the model’s 
OUR is lower, with 2±lmol02 m _2 yr _1 below 7 ” = 26.44. In the subtropical North 
Atlantic, OUR is not significantly different from zero beneath 7 71 = 27.16 (0.3 ±0.3 mol02 
m -2 yr -1 ) while there is a significant influx above that surface (1.1±0.2 mol 02 m -2 yr -1 ) 
(The uncertainty is small compared with the uncertainties of the oxygen transports at 
24°N and 48°N because of the correlations between errors, e.g., page 146). This contrasts 
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Location 

Method 

OUR 

New Prod. 

Reference 

Subtrop. N. Atl. 

Tritium budget 

2.4 to 3.5 


Sarmiento et al. (1990) 

Subtrop. N. Atl. 

228 Ra 

8.5±0.8 


Sarmiento et al. (1990) 

54°S-45°N 

hydrographic 
based model 

2.5±0.8 


Riley (1951) 

Sargasso Sea 

Tritium/ 3 He 

7.2 


Jenkins (1980) 

30° N, 30°W 

He age 

6.4±0.3 


Jenkins (1987) 

64°N, 28°W 

O 2 production 
in mixed layer 


5.2 

Peng et al. (1987) 

near Bermuda 

O 2 balance 


3 to 4 

Musgrave (1989) 

near Bermuda 

O 2 balance 


5.6±1.5 

Spitzer and Jenkins (1989) 


Table 4.5: Estimates of oxygen utilization rates and new production (mol yr 2 m 2 
adapted from Sarmiento et al. (1990). 


with the large average OUR of Sarmiento et al. (1990). According to GCM simulations, 
most of the oxygen input occurs in regions of large oceanic cooling such as the Gulf-Stream 
north of 36°N (e.g., Stephens et al., 1998). The oxygen influx taking place there may dom¬ 
inate the whole signal between 24° N and 48° N. Such a signal is probably underestimated 
in the Sarmiento et al. (1990) average because of limited data coverage. 

In the North Atlantic between 24°N and 36°N, Rintoul and Wunsch (1991) reported a 
net, “barely significant”, loss of oxygen of 1.3 molO 2 m _2 yr _1 . But again, the gain that 
we estimate between 24°N and 48°N may be dominated by the large oxygen influx north 
of 36°N, which is not part of the Rintoul and Wunsch (1991) estimate. 

Following the work of Watson et al. (1995), Holfort et al. (1999b) estimated the solu¬ 
bility pump as about 1500kmol O 2 s _1 per PW (Petawatt). The North Atlantic heat loss 
is, between 24°N and 48°N, 0.5 ±0.1 PW (Fig. 4-2), which would correspond to an oxygen 
influx of 750±150 kmol O 2 s -1 . The influx of oxygen diagnosed by the inverse model in the 
surface layer is in close agreement, at 620 ± 100 kmol O 2 s -1 . There is a similar agreement 
in the South Atlantic between 30°S and 45°S, with 450 ± 150 kmol O 2 s _1 predicted from 
the solubility pump versus 680±60 kmol O 2 s _1 diagnosed from the model. Globally, how¬ 
ever, the solubility effect is significantly smaller than the net observed oxygen sources and 
sinks, as shown by Fig. 4-11, suggesting that oxygen exchanges driven by biogeochemistry 
are of major importance in the ocean-roughly 50%. 

The Atlantic and Indo-Pacific fluxes are compared with other results in Figure 4-12. 
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There is an overall agreement within an uncertainty of one standard deviation, except 
for the Keeling and Peng (1995) estimate which is directed to the north at the Equator. 
The agreement in both region between the various estimates, the MIT general circulation 
off-line model, as well as the Princeton Ocean Biogeochemistry Model is remarkable (and 
perhaps fortuitous, given the uncertainties). 

In conclusion, the estimated oxygen fluxes and flux divergences reflect both solubility 
and remineralization effects. However, oxygen input through photosynthesis could not be 
distinguished from other effects in the surface layers because of the large uncertainties in 
the surface nutrient divergences. Large and significant net oxygen flux divergences are 
found, with net efflux in the tropical regions and net influx at high latitudes. The global 
oceanic oxygen budget is in balance. 
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Figure 4-12: Atlantic and Indo-Pacific dissolved oxygen transports (stars). MITgcm=MIT 
General Circulation Model (offline biogeochemical code; McKinley et al., manuscript in 
preparation, 1999); POBMl=Princeton Ocean Biogeochemistry Model Version 1 (N. Gru¬ 
ber, personal communication, 1999). 
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Chapter 5 


Conclusions 

5.1 Global summary 

The average oceanic circulation and property transports were estimated over the period 
1985-1996, based on the WOCE hydrographic program and the JADE program in the 
Pacific-Indonesian Throughflow region. Simple dynamics and budget constraints were 
used to produce an estimate that is consistent, within estimated errors, with hydrographic 
and current-meter measurements, oceanic dynamics and biogeochemistry. 

The physical oceanographic circulation was analyzed with emphasis on the meridional 
overturning. It is summarized on Figure 5-1. Large scale water mass transports were 
found to be similar to previous estimates from earlier data in the Atlantic and Indian 
Oceans. In the Pacific Ocean, differences with previous values were found: a noticeably 
weaker overturning in the North Pacific Ocean. Deep water is being formed at a rate of 
16 ±2 Sv in the North Atlantic, while 20±6 Sv of Circumpolar Deep Water is being formed 
in the Southern Ocean, the latter essentially supplied by the former. Through the use of 
anomaly equations, significant dianeutral diffusivities associated with tracer transports 
were found in the deep layers of most oceanic regions, averaging globally 3 ± 1 cm 2 s -1 in 
the 2000-4000 dbars depth range, excluding the Southern Ocean. In the Southern Ocean, 
significant diffusivities could not be determined because of the large uncertainties in layers 
that are in contact with the surface. 

The freshwater flux divergences are 100% uncertain in most places, but there is an 
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Figure 5-1: Global circulation summary. The water transports are indicated (approxi¬ 
mately) for the different density classes. The exact boundaries between surface, deep and 
bottom waters are 7 n = 27.82/28.0986 (North Atlantic); 27.72/28.11 (South Atlantic); 
27.7/28.11 (Indian); 27.8/28.10 (Pacific). The color of the upwelling arrows indicate the 
layer from which the water is coming. In the Southern Ocean, while the bottom water 
formation takes place mostly in the Weddell Sea, upwelling is uncertain and spread over 
the different regions. The South Pacific transports are given at 17°S because of the more 
complicated structure at 32°S (Fig. 3-7). 
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overall agreement with previous estimates. Energy transports evidence significant ocean- 
atmosphere heat fluxes, with a global heating of 2.3 ± 0.4 PW in the tropical band. Most 
of the cooling occurs in the Northern Hemisphere (1.7 ± 0.2 PW) while in the Southern 
Hemisphere only 0.7 ± 0.3 PW is transferred to the atmosphere. 

Silica, nitrate and phosphate fluxes and divergences were analyzed with respect to their 
respective biogeochemical cycles. Determination of nutrient divergences appears to be at 
the edge of the model and data capabilities, but despite the large uncertainties, the mag¬ 
nitude of the maximum export production agrees closely with independent measurements 
from moorings and isotopic data. The (possible) signature of a large-scale average export 
production is found for silica, nitrate and phosphate in several temperate regions. The 
nutrient divergences thus provide an upper bound on the average remineralization rate 
in several depth ranges while the lower bound is more uncertain. Net nutrient sources or 
sinks are found in several regions. Such divergences suggest advection of dissolved organic 
matter across the sections, as first pointed out by Rintoul and Wunsch (1991), unless those 
are spurious as a result of aliasing the seasonal cycle. 

Significant exchanges of oxygen with the atmosphere are found, with oxygen intake 
at high latitudes and loss at low latitudes. About half of the oxygen gain in the surface 
layer can be explained by changes in solubility through temperature changes, leaving an 
important role to the biology-driven exchanges. Indeed, significant oxygen utilization is 
found below the surface in several regions, indicating remineralization. 

5.2 Discussion and perspectives 

Exchanges of carbon dioxide between ocean and atmosphere are related to the oxygen 
exchanges. While the solubility pump produces influxes and effluxes of CO 2 that are 
similar to O 2 , the biological pump creates exchanges in the opposite direction. With 
the oxygen budget being in balance over the whole ocean, there is a suggestion that the 
solubility and biological CO 2 exchanges are in balance as well-at least within uncertainties. 
On the other hand, oceanic uptake of anthropogenic CO 2 may occur as a simple partial 
pressure equilibration, independently from the biological and solubility pumps, as has been 
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suggested by several authors (e.g., Sarmiento et al., 1992). 

Combination of ocean dynamics with geochemical cycles represents one more step 
towards producing a consistent estimate of the oceanic circulation. Use of a simple time- 
average model permits estimation of meridional property fluxes based on the highest 
resolution allowed by hydrographic data. Such a model permits estimation of the true 
uncertainties of the property fluxes. The physical circulation estimated here during the 
WOCE period provides a reference for climate studies, for instance as a way to test general 
circulation models. The estimated ocean-atmosphere heat and freshwater transfers can be 
used, similarly, as a benchmark for climatology, while average dianeutral diffusivities will 
help in parameterizing the sub-grid-scale processes in numerical models. 

I chose a standard, or preferred, solution among a large number of experiments (about 
100 in total), but there is no limit to possible configurations of the model and detailed 
regional studies will improve the model consistency. 

But a time-average budget relies upon the ergodic hypothesis, and the remaining un¬ 
certainty is dominated by oceanic variability in the density and other properties. The 
uncertainties presented here account for the errors from both the model limitations and 
the measurements. In particular, the alias errors in the mass transports resulting from 
non-synopticity of the sections and oceanic variability is taken into account by reference 
to a general circulation model simulation. 

The analysis of the present uncertainties and model limitations suggests several direc¬ 
tions for future work. The necessary next step includes the use of additional data, some 
of which may already exist (e.g., current-meters, floats, absolute altimetric data) com¬ 
bined with more complete dynamics to eventually reduce the residual uncertainty in the 
solutions. In particular, the property fluxes through the Pacific-Indonesian Throughflow 
(PIT) and other highly variable currents such as the North Brazil Current are a large 
source of uncertainty, due to the lack of measurements and the large seasonality there. 

5.2.1 New data 

Introduction of new data can be simulated in order to optimize their impact on the re¬ 
duction of the existing uncertainties. Naturally, new high-quality surveys will increase the 
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resolution of the flux divergence estimates. But if the scope is to estimate, for instance, 
air-sea heat exchanges or export production, such surveys could be designed in order to 
optimize the signal, based on existing maps of heat exchanges from climatologies or atmo¬ 
spheric models. The present surveys are mostly zonal or meridional, and several “boxes” 
contain regions of strong heating and regions or strong cooling that mostly cancel out. As 
a result, the inverse model diagnoses only the small residual air-sea exchange. 

Inclusion of the existing meridional hydrographic sections in the present model will be 
a considerable challenge. Such work is in progress in several oceanic areas. The major 
obstacle to combining zonal and meridional sections is the treatment of crossing points, 
which can lead to large mass imbalances in individual layers due to temporal changes in 
the eddy field from one section to the next (i.e., Wunsch and Grant, 1982). The a priori 
noise in the mass flux has to be modified accordingly. 

5.2.2 Improvements to the present model 

The present circulation estimate will improve with time as it benefits from suggestions 
for improvement on the different parameter choice such as the initial reference surface 
based on new current meter measurements and water mass analysis where the signals are 
important such as in the boundary current regions. The experiments that I did (e.g., in 
the Indian Ocean) suggested that large scale mass transports are relatively insensitive to 
such changes, although there is no guarantee that the same is true everywhere. 

The parameterization of dianeutral exchanges as a single diffusion and advection coef¬ 
ficients over large oceanic areas may be unrealistic. In particular, the possibility of double 
diffusion was not addressed, and use of a different coefficient of heat and tracer transfers 
may be used in future versions of the model. Dianeutral exchanges are also poorly param¬ 
eterized in regions of outcrop such as the Southern Ocean. There, an advection-diffusion 
scheme is probably unrealistic and an alternative parameterization may be to assign an 
independent transfer coefficient to each property (e.g., Sloyan and Rintoul, 1999a,b), but 
requiring extra data to resolve them separately. 

The calculation of Ekman property transports may be improved as well by, for instance, 
estimating the transport across each station pair and explicitly using the vertical property 
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distribution (e.g., Wijffels et al., 1999). 

Lastly, the a priori statistics on the model can be improved. For instance, the weight 
of anomaly equations can be directly calculated from a simulation by a general circulation 
model as I did for the mass fluxes, instead of the ad-hoc estimate that was used here. Vari¬ 
ability in the nutrient fluxes is another major possible alias that could not be quantified. 
Seasonal correlations between changes in the baroclinic transports and nutrient concentra¬ 
tion in the upper layers can create such aliases in the nutrient transport computed from a 
single hydrographic survey. Repeat survey will provide useful insights for such variability, 
would the measurement be made only over the upper layers that are the most sensitive to 
seasonal cycle. 

5.2.3 Improving the dynamics 

Including complete dynamics is the ultimate target of ocean estimation. In this model, 
variability in the density field is a major source of uncertainty and limits the impact of 
new data. Allowing a correction to the density field permits achievement of better mass 
conservation and optimizes the impact of new data. Simple steady-state (or time-average) 
models with corrections to the density field allowed through the inversion have been used 
in the past (e.g., Wunsch, 1994). But the major obstacle, apart from computational load, 
is the lack of knowledge about a priori statistics for the density variability. Global spectra 
of the density field variability may partially solve the problem (Zang, 1999), although the 
variability has a strong spatial dependence. 

Data assimilation in general circulation models is the natural next step toward ocean 
circulation estimation. Computer capabilities are currently limiting the resolution and 
prohibiting calculation of uncertainties. Assimilation of data such as hydrographic sections 
remains a considerable challenge, but the ongoing efforts, in parallel with simple statistical 
estimates such as the one presented here, will keep improving our estimate of the time- 
varying ocean circulation. 
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Appendix A 


Ageostrophy in the General 
Circulation Model 


While geostrophy is generally a reliable approximation in the ocean, the Semtner and 
Chervin model output contains substantial ageostrophic transports. Their cause and na¬ 
ture are investigated in this Appendix. 

The ageostrophic component at 36° N was computed as the difference v ageos t = v mo d e i — 
v ga , where v ga is the absolute geostrophic velocity defined by (2.8). The GCM uses a 
staggered B grid (e.g., Semtner, 1986; Killworth et al., 1991), so that v ga is staggered 
with respect to temperature, T, salinity, S and sea surface elevation, 77 (Figure A-l). In 
principle, v ga can be compared with the model velocity, v m , on the original u,v grid. 
Equivalently, the comparison was done on an intermediate grid, called the sample grid 
(superscripts s) on to which all the state variables are linearly interpolated. Namely, the 
model u,v, noted (u m ,v m ), were interpolated at the temperature points in the East-West 
direction, yielding the sample field (u s ,v s ) while T, S and 77 were interpolated in the 
North-South direction, so that all data stand on the s grid (Fig. A-la). At topography, 
the geostrophic velocity was assumed to be zero before interpolation. This corresponds to 
the no-slip boundary condition for model velocity. 

The ageostrophic velocity was then computed at each “sample” grid point 
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Figure A-l: Model B grid and sampling for a zonal section. Subscript m stands for the 
original model (u,v) grid; subscripts for the points at which the fields are interpolated 
(’’sampled”), and <7 for geostrophic. T,S,rj are interpolated in the North-South direction, 
a) top view; b) vertical and zonal section. 


(superscripts): 


Kgeost = - V U- ( A1 ) 

Locally, the ageostrophic velocities are noisy and difficult to interpret, and zonal or ver¬ 
tical integrations are helpful for determining the transports structure. Figure A-2 shows a 
snapshot of the cumulative ageostrophic transport through each vertical grid interval, inte¬ 
grated from the eastern end at 36°N, revealing strong variations of ±5 Sv (1 Sverdrup = 10 6 
m 3 s _1 ) at the western end. 

A frequency spectrum (not shown) of the net ageostrophic transport below the Ekman 
layer shows a very large amount of energy at high frequencies. It turns out that most of 
the ageostrophic transport is due to inertial and possibly internal waves. 

Because the GCM output was subsampled only every third day, inertial waves are 
aliased into frequencies that vary with latitude, as discussed in Jayne and Tokmakian (1997). 
A 15-day low pass filter, of the form —sine (yjr) (where r is the time lag) was applied 
to all the state variables of the GCM output to filter out high frequencies including the 
inertial waves. (At 36°N, the frequency of an inertial oscillation is sq = 1.176 day -1 so 
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Figure A-2: Cumulative ageostrophic transport in each layer of the model, integrated from 
the eastern end, section A36°N. 

that the 3-day alias frequency is s a u as = (6.34day) -1 .) FigureA-3 shows the integrated 
ageostrophic transport in a snapshot of the filtered field: it is now dominated by the Ek- 
man transport. In the interior and away from the western boundary current, ageostrophic 
transport accumulates up to at most a few tenths of a Sverdrup in each layer. Near the 
western boundary current (WBC) region, the integrated value increases up to ±1 Sv. 

The low-pass filtered, vertically-integrated absolute and geostrophic transports through 
each grid cell are displayed on Figure A-4a for one time period. Below the Ekman layer, 
the ageostrophy is of less than 1 Sv (Fig. A-4b, dotted lines), showing that the model 
was close to geostrophic equilibrium in the interior. However, near the Gulf-Stream the 
departure attains 5 Sv per 36-km grid cell, which is up to 50% of the total transport in 
the cell. In this ageostrophic region, the model velocity transports are weaker than the 
geostrophic ones (Fig. A-4a). 

When integrated from the East (Fig.A-4c), the cumulated ageostrophic transport 
grows up to 5Sv, and is more or less cancelled in the WBC region, to yield a total 
ageostrophic flux of 2Sv below the top 100 m. Nevertheless, this cancellation is accidental 
to this snapshot and systematic deviations of up to 5 Sv are observed persisting for a few 
weeks. There is no apparent correlation between the ageostrophy and the topographic 
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Figure A-3: Same as previous figure, but with inertial oscillations filtered out. 



features (Fig. A-4d). 

What process is responsible for a 50% deviation from geostrophy ? The GCM equation 
for u is 


du du du , dri _ 4 

m +a »i + v ^- lv = - 9 ai + ^ “■ 


At the surface, the Rossby number (based on the 36 km grid-step) is 


€ fL ° (lO- 4 x 36km) “°' 15 ' 

A direct calculation of (u|^) / (/ v) from the model velocity and horizontal shear pro¬ 
duces a similar value. The non-linear term is not able to explain the ageostrophic trans¬ 
port. If the Gulf-Stream curvature occurs on similar scales, the cyclostrophic term is of 
the same order as the non-linear terms, and is thus small in the model. 

The horizontal hyperviscosity coefficient in the model is about 10 12 m 4 s _1 (Stammer 
et al., 1996) so that the horizontal Ekman number ([viscous term]/[geostrophic term]) is 
of order: 

uV 4 u v 10 12 , 

- rv_, - — — - -———- ^ g X 10 ^ 

fv fL 4 10- 4 x (36 km ) 4 ’ 
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Distance (km) 

Figure A-4: Vertically integrated transport per grid interval from the a) model absolute 
velocity and geostrophic velocity, from 100 m to bottom; b) ageostrophic velocity and 
c) ageostrophic velocity integrated from East, d) GCM bottom topography. Inertial waves 
were removed by filtering. 
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which is too small to explain the observed ageostrophy. 

The approximations of the geostrophic calculation near a sloping bottom (which is a 
stair-case in the model) is also unlikely to explain the ageostrophy as the bottom topog¬ 
raphy did not show any correlation with the ageostrophic transport. 

An examination of the vertical structure of the velocity in the Gulf Stream shows 
that the GCM shear is close to the thermal wind, with the ageostrophic transport be¬ 
ing essentially depth-independent. (The thermal wind is referenced to the surface, and 
the geostrophic surface velocity calculated from ^ is used to produce the geostrophic 
transport.) For a 36 km grid cell, a deviation of 3cms -1 over the water column (5000 m) 
produces the 5 Sv observed deviation. However, this deviation is not justified by any of 
the terms in the momentum equations. 

The strong ageostrophy observed in GCM transport is most likely due to numerical 
problems. One possibility is the presence of an enhanced bottom friction in the barotropic 
equation of the GCM. In principle, the exact source of ageostrophy is identifiable by 
examining the size of each term in the original numerical equations. This requires a more 
complete GCM extraction than the one we have and is beyond the present work. In any 
case, the depth-independent nature of the ageostrophy makes us confident in using the 
GCM output to estimate baroclinic variability. 

An alternative explanation is purely physical: in the presence of a strong current, the 
Ekman flux is enhanced through the non-linear interaction with the wind stress. The 
sea surface elevation may be slightly shifted from the elevation due to the underlying 
geostrophic currents (e.g. Niiler, 1969). In the GCM, there is no explicit mixed layer, 
but the diffusivity is dependent on the Richardson number so that one can have a strong 
mixing is the first layers. Estimating the amplitude of this effect is beyond the scope of 
the present work. 
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Appendix B 


Model Error in a Two-dimensional 
Ocean 

A hypothetical two-layer, two-dimensional ocean is used to illustrate the effect of the 
model error on the representativeness of the heat and silica flux calculation. The “true” 
circulation is sketched in Figure B-la: the basin is closed to the north of the section so that 
it is known that net mass losses or sources do not exceed the freshwater flux, known to an 
accuracy of ±0.5 Sv. The true circulation is a 20 Sv overturning. The upper layer is 15° 
water with a silica concentration of 10/imol Si/kg while the lower layer is 2°, 30^/mol Si/kg. 
The actual heat flux divergence is given by 

H true — I pC p vO dx dz 

~ 1000 x 4000 x 20 x 10 6 x (15 - 2) kg ro~ 3 J kg- 1 K~ x m 3 s -1 K 
^1 PW; (B.l) 

and the silica divergence is: 

Si^true — / pvC dxdz 

= 1000 X 20 x 10 6 X (10 - 30) x 10~ 6 
= —400kmol Si s _1 . (B.2) 


Odb 

2000 

4000 

Figure B-l: Two-dimensional overturning example model. In the true, or time-average 
case (a) the upper layer carries 20 Sv northward of 15°C water with a silica concentration of 
1 Op,mol Si kg -1 while the lower layer carries 20 Sv of 2°C, 30pmol Si kg -1 water southward. 
Because of the baroclinic variability, a mass excess of 6 Sv is observed in the upper layer 
(b), leading to the mass-conserving circulation (c). 

(The silica flux is a unrealistic in this example as it should be small in the real ocean, 
which is not possible in this 2-D model.) The goal is to recover this solution using an 
inverse model requiring mass conservation. Because the model is two-dimensional, the 
barotropic correction to mass imbalances produces an adjustment in the property fluxes 
of Am x (Ci + C 2 )/ 2 , where Cj and C 2 are the concentrations in the upper and lower 
layers, and Am is the adjustment to the total mass flux. 

1) no measurement errors 

First, suppose that there were no errors in the measurements, (no baroclinic variability, 
no internal waves) and that we picked a reference level that is equal to 2000 dbars and 
that exactly reproduces the true circulation. Then, there is no initial mass imbalance and 
the inverse model produces the true solution. The derived heat and silica flux are the true 
ones. An a priori uncertainty in Am of 5m = ±lSv in the zero-mass-flux equation implies 
an uncertainty on the property fluxes of 

SH ±1 = P C p 5 M (T 1+ T 2 )/2 

= 1000 x 4000 x 10 6 x (15 + 2) /2 = ±0.034 PW 
5Si± 1 = 10 6 x (10 + 30) /2 x 10" 6 = ±20 kmol S* s' 1 . (B.3) 
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2) baroclinic error of 5 Sv 


Now, suppose that for various reasons the measured/calculated shear is off by 6Sv with 
respect to the time average circulation, or the “true” circulation (Figure B-lb). With the 
mass conservation constraint alone, the inverse model adds 3 Sv to each layer, producing 
the circulation on Figure B-lc. The resulting fluxes are: 

H inv = 4 x 23 x 10 12 x (15 - 2) = 1.2 PW 

Si,inv = 23 x (10 - 30) = -460kmolS z s' 1 . (B.4) 

If one would have constrained the circulation to conserve mass within ±1 Sv, the fluxes 
would be 


H inv = 1.2 ± 0.034 PW 

Si,inv — —460 ± 20kmol Si s -1, (B.5) 

which is clearly incompatible with the true solution (B.1-B.2). On the other hand, an a 
priori uncertainty of ±10 Sv in the mass equation, would yield 

H inv = 1.2 ±0.34 PW 

Si,inv = —460 ± 200kmolSiS _1 , (B.6) 

which is realistic as the model error is properly taken into account. 

If the model is not able to conserve mass properly because of other constraints so that 
the initial mass imbalance of 6 Sv remains, the solution is: 

H inv = 4 x 10 12 x (26 x 15 - 20 x 2) 

H inv = 1.4 ± 0.34 PW 
Si,inv = 10 x 26 - 30 x 20 

Si t i nv = —340 ± 200 kmol Si s -1 (B-7) 

which is consistent with the true solution (1PW, -400kmolSi s _1 ). (The uncertainties 
correspond to one standard deviation.) 
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3) anomaly formulation 

In this particular case, the anomaly formulation yields the exact answer: 

H inv = 1.4 - 4 x 10 12 x (15 x 6) 

H mv = 1±0.26 PW 
Si,inv = -360 - 6 x 10 

Si,inv = -400 ± 100 kmol S,- s _1 (B.8) 

Note the importance of using the average property concentration in the layers where the 
mass residual occurs rather than having a global average property concentration. 
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Appendix C 


A priori mass uncertainties 


The first table gives the basic a priori uncertainty contribution from each section. The 
second table gives the a priori uncertainty used in each box. (For comparison, the un¬ 
certainty including the Ekman transport uncertainty are indicated in the shaded rows. 
Because an adjustment to the Ekman transport is permitted, those are not used in the 
inverse model.) 
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Figure C-l: A priori uncertainty contribution from each section, in Sv rms. The “total 
net” (5th row) indicates the total uncertainty in the net flux. 
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Figure C-2: A priori uncertainty on each mass conservation constraint. 
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Appendix D 


Standard solution absolute 
velocities and cumulative 
transports 


For each section: 

(*) The first figure shows the absolute velocity field across the section. The initial 
reference surface is indicated by the dash-dot line on the absolute velocity plots. The 
velocities are positive northward and eastward. Negative velocity regions are shaded. 

The contour intervals are lcms -1 (thin contours, from -lOcms -1 to lOcms -1 ) and 10 
cms -1 (thick contours for stronger flows). The reference velocity (or velocity on the 
reference surface after integration) is given in the panel above the velocity contour plots, 
the shaded area giving the uncertainty (one standard deviation). 

(it) The second figure gives the transports integrated over several layers and between the 
neutral surfaces that are indicated-those roughly correspond to the water masses defined 
in the text. The thick plain line gives the absolute transport; the dashed line gives the 
contribution from the thermal wind; and the dots give the contribution from the 
reference velocity adjustment. The shaded areas give the uncertainty on the absolute 
transport (one standard deviation). The Ekman transport is not included. 
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Figure D-8: Cumulative transports, Indian Ocean, I2W (45°E) 
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Figure D-10: Cumulative transports, Indian Ocean, 12 (8°S) 
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Figure D-12: Cumulative transports, Indian Ocean, 110 (111°E) 




































-7000 - 



100 











































































donia). 


213 



















220 

21E Longitude 

























































Figure D-20: Cumulative transports, Pacific Ocean, P21E (17°S, New Caledonia to South 
America). 
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Figure D-23: Absolute velocities, Pacific Ocean, PI (47°N). 


218 




















PI, surface tor =26 





















220 


























































































-5000 


i 


I 


CO 

CO 


to 

L 

01 

o 

01 

o 

Ol 

o 

o 

o 

o 

o 

o 

o 

o 

o 

o 


I 

01 

o 

o 


cm s 


-i 


i 

01 


Ol 

















































-6000 




























J_I_I_I_I_I_I_L 

r =28.13 to bottom 

i i i ! i i i r 


_i_i_i_i_i_i_i_i_ 

30W 25W 20W 15W 10W 5W 0 5E 1 


Figure D-32: Cumulative transports, Atlantic Ocean, A7 (4.5°S). 
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Figure D-49: Absolute velocities, Southern Ocean, P12 (140°E). 
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Figure D-52: Cumulative transports, Southern Ocean, P14S (170°W). 
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Appendix E 


Standard solution property fluxes 
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Section 

Temp. 

Oxygen 

A2 

0.6 ±0.1 

-1750 ±500 

A5+F1. C. 

1.3 ±0.15 

-2070 ± 600 

FI. Current 

2.4 ± 0.07 

5260 ± 150 

A6 

0.7 ±0.3 

-1430 ±950 

A7 

1.0 ±0.55 

— 138 ±1568 

A8 

0.9 ±0.4 

-1380 ± 1250 

A9 

0.8 ±0.2 

-1590 ± 700 


0.35 ±0.15 

430 ± 680 

All 

0.66 ±0.12 

-52 ± 710 

A12 

0.75 ± 0.23 

32175 ± 1900 

A21 

1.30 ±0.14 

31385 ± 1230 

16 

0.90 ±0.6 

31555 ± 1900 

15 

-1.50 ±0.23 

-1820 ± 1090 

I9S 

1.72 ±0.3 

34248 ± 1840 

P12 

1.42 ± 0.24 

34345 ± 1800 

P14S 

1.64 ±0.26 

35551 ± 3260 

P6 

0.6 ±0.2 

4725 ±1000 

P21 

0.21 ±0.46 

3240 ± 1100 

P3 

0.53 ± 0.24 

-1410 ±680 

JADE 89-92 

-1.36 ±0.14 

-2248 ± 407 

PI 

0.0 ±0.05 

-672 ± 450 

14 

-0.77 ±0.2 

-2114 ±1126 

13+14 

-1.79 ±0.44 

-1416 ±1210 

I2W 

-0.77 ±0.8 

-3500 ± 1670 

I2W+I2+I10 

-1.57 ±0.8 

-1758 ± 1750 


Table E.l: Standard solution fluxes of temperat 
oxygen, nitrate and phosphate (kmol/s). Positi 
for Florida Current. 
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Phosphate 

Nitrate 

Silica 

1.1 ±2.5 

8 ±35 

-130 ±50 

-7.6 ±3.6 

-55 ± 55 

-220 ± 80 

14 ± 1.2 

260 ± 20 

154 ± 13 

-1.0 ±7 

-75 ± 120 

-165 ± 107 

-11 ± 10 

-210 ± 150 

-190 ± 170 

N/A 

-180 ± 140 

-230 ± 154 

N/A 

-60 ± 80 

-134 ± 107 

-4 ±5 

-70 ± 80 

-180 ±140 

-9.0 ±6 

-240 ± 90 

-107 ±190 

294 ± 17 

4304 ± 245 

9800 ± 900 

302 ± 13 

4340 ± 190 

9730 ± 880 

295 ± 27 

4095 ± 340 

9810 ± 920 

-4.3 ±9 

-85 ± 130 

250 ± 375 

320 ± 18 

4620 ± 264 

9547 ± 920 

330 ± 20 

4700 ± 266 

9520 ± 920 

326 ± 30 

4815 ± 440 

10330 ± 1300 

14 ±10 

186 ± 136 

-275 ± 380 

25 ±12 

270 ± 170 

-130 ±416 

0.8 ±8 

-1 ± 100 

46 ± 344 

-12 ±6 

-155 ±80 

-390 ± 190 

4.7 ±3 

76 ±44 

66 ± 150 

—19 ± 11 

-261 ± 150 

-311 ±311 

—19 ± 10 

-240 ± 150 

33 ± 440 

-10 ± 16 

-123 ±218 

-178 ±516 

-28 ±15 

-400 ± 210 

-260 ± 530 


5 (expressed in PW on the Celsius scale), 
is northward or eastward. FI. C. stands 
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